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all you need 


w store NITROGEN TETROXIDE 


is a rocket 


This liquid-fuel oxidizer needs no refrigeration, causes no freeze-ups 


an oxidizer for | 

Nitrogen Tetroxide has even more 

to recommend it than its high per- | 

formance (99% of theoretical I.» 

hypergolic at low altitude). The 

ease with which it can be stor 
and handled offers another major 
advantage. 
requires no refrigeration, 
no high-pressure vessels. It is non 
_ corrosive, can be stored indef 

-nitely in plain carbon steel tan 

at the launching site or right in the 

rocket itself. And it’s ready when 

it's needed—there are no freeze- 
ps in valves and motors with N.O, 
Availability: excellent. Allied i 

a major producer of Nitrogen Tet. 

roxide, can ship immediately in 

cylinders or tank car lots. Allied 
also produces ammonia, ethylene 
oxide and methanol for the rocke 
industry. Write for technical or 
other information you desire 
any of these products. 

For specifications and local offices, see 
NITROGEN DIVISION 435 and in Chemical Week 
Dept. NT 5-12.2, 40 Rector Street, New York 6, N. 


4 

; 

3 

= 

477 


JOURNAL 


A PUBLICATION OF THE AMERICAN ROCKET SOCIETY 


cormency 


EDITOR Martin Summerfield 
ASSISTANT EDITOR _ Barbara Nowak 
ART EDITOR John Culin 


ASSOCIATE EDITORS 


Ali Bulent Cambel Northwestern University 
Irvin Glassman Princeton University 
M. H. Smith Princeton University 


CONTRIBUTORS 
Marshall Fisher Princeton University 


George F. McLaughlin 


ADVERTISING AND PROMOTION MANAGER 
William Chenoweth 


ADVERTISING PRODUCTION MANAGER 
Walter Brunke 


ADVERTISING REPRESENTATIVES 


dA t James C. Galloway and Co. 

135 y N. 6535 Wilshire Blvd., Los Angeles, Calif. 
Telephone: Yukon 6-6855 Telephone: Olive 3-3223 
i dA iat R. F. Pickrell and Associates 
Waker De 318 Stephenson Bldg., Detroit, Mich. 
Telephone: Andover 3-1154 Telephone: Trinity !-0790 
Lo J. Bresnick John W. Foster 
504 Washington Ave., Chelsea 50, Mass. 239 4th Ave., Pittsburgh, Pa. 
Telephone: Chelsea 3-3335 Telephone: Atlantic |-2977 


American Rocket Society 


500 Fifth Avenue, New York 36, N. Y. 
Founded 1930 


OFFICERS 


President John P. Stapp 
Vice-President Howard S. Seifert 
Executive Secretary James J. Harford 
Treasurer Robert M. Lawrence 
Secretary and Asst. Treasurer A. C. Slade 
General Counsel Andrew G. Haley 
Director of Publications Irwin Hersey 


BOARD OF DIRECTORS 


Terms expiring on dates indicated 


James R. Dempsey 1961 Simon Ramo _ 1960 
Alfred J. Eggers Jr. 1959 H. W. Ritchey 1959 
Krafft Ehricke 1959 William L. Rogers 1959 
Samuel K. Hoffman 1960 David G. Simons 1961 
J. Preston Layton 1960 John L. Sloop 1961 
A. K. Oppenheim 1961 Martin Summerfield 1959 
William H. Pickering 1961 Wernher von Braun’ 1960 


Maurice J. Zucrow 1960 


TECHNICAL COMMITTEE CHAIRMEN 


Lawrence S. Brown, Guidance and David B. Langmuir, Ion and Plasma 
Navigation Propulsion 
Milton U. Clauser, Hydromagnetics Y. C. Lee, Liquid Rockets 
Kurt H. Debus, Logistics and Operations Max Lowy, Communications 
William H. Dorrance, Hypersonics Harold W. Norton, Test Facilities and 
Herbert Friedman, Instrumentation Support Equipment 
and Control Paul E. Sandorff, Education 


George Gerard, Structures and Materials 


Milton Greenberg, Physics of the 
Atmosphere and Space 


Stanley V. Gunn, Nuclear Propulsion 
Andrew G. Haley, Space Law and 


William Shippen, Ramijets 
John L. Sloop, Propellants and 
Combustion 


Ivan E. Tuhy, Solid Rockets 


Sociology Stanley White, Human Factors 
Samuel Herrick, Flight Mechanics George F. Wislicenus, Underwater 
Max Hunter, Missiles and Space Propulsion 

Vehicles Abe Zarem, Non-Propulsive Power 


478 
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ment of astronautics through the dissemina- 
tion of original papers disclosing new scientific 
knowledge and basic applications of suc}: 
knowledge. The sciences of astronautics are 
understood here to embrace selected aspects 
of jet and rocket propulsion, space flight 
mechanics, high-speed aerodynamics, flight 
guidance, space communications, atmospheri: 
and outer space physics, materials and struc- 
tures, human engineering, overall system 
analysis, and possibly certain other scientific 
areas. The selection of papers to be printed 
will be governed by the pertinence of the topic 
to the field of astronautics, by the current o 
probable future significance of the research 
and by the importance of distributing the in 
formation to the members of the Society and 
to the profession at large. 
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Manuscripts must be as brief as the prope 
presentation of the ideas will allow. Ex 
clusion of dispensable material and concise 
ness of expression will influence the Editors 
acceptance of a manuscript. In terms o 
standard-size double-spaced typed pages, 
typical maximum length is 22 pages of tex: 
(including equations), 1 page of references 
1 page of abstract and 12 illustrations 
Fewer illustrations permit more text, and vic: 
versa. Greater length will be acceptabl: 
only in exceptional cases. 

Short manuscripts, not more than one 
quarter of the maximum length stated for ful! 
articles, may qualify for publication as 
Technical Notes or Technical Comments 
They may be devoted to new developments 
requiring prompt disclosure or to comments 
on previously published papers. Such manu- 
scripts are usually published within two 
months of the date of receipt. 

Sponsored manuscripts are published 
occasionally as an ARS service to the indus- 
try. A manuscript that does not qualify for 
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requirements as to subject, scope or length, 
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distribution among jet propulsion engineers, 
may be printed as an extra part of the Journal 
or as a special supplement, if the author or 
his sponsor will reimburse the Society for 
actual publication costs. Estimates are 
available on request. Acknowledgment of 
such financial sponsorship appears as a 
footnote on the first page of the article. 
Publication is prompt since such papers are 
not in the ordinary backlog. 

Manuscripts must be anubie spaced on one 
side of paper only with wide margins to allow 
for instructions to printer. Include a 100 to 
200 word abstract. State the authors’ 
positions and affiliations in a footnote on the 
first page. Equations and symbols may be 
handwritten or typewritten; clarity for the 
printer is essential. Greek letters and unusual 
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lower case letters, and indicate subscripts and 
superscripts. References are to be grouped at 
the end of the manuscript and are to be given 
as follows: For journal articles: authors first, 
then title, journal, volume, year, page numbers 
for books: authors first, then title, publisher, 
city, edition and page or chapter numbers. 
Line drawings must be clear and sharp to make 
clear engravings. Use black ink on white 
paper or tracing cloth. Lettering should be 
large enough to be legible after reduction. 
Photographs should be glossy prints, not 
matte or semi-matte. Each illustration must 
have a legend; legends should be listed in 
order on a separate sheet. 

Manuscripts must be accompanied by 
written assurance as to security clearance in 
the event the subject matter lies in a classified 
area or if the paper originates under govern- 
ment sponsorship. Full responsibility rests 
with the author. 

Submit manuscripts in duplicate (original 
plus first wage? with two sets of illustra- 
tions) to the Editor, Martin Summerfield, 
Professor of Aeronautical Engineering, Prince- 
ton University, Princeton, N. J. Preprints 
of papers presented at ARS national meetings 
are automatically considered for publication. 
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ve “Operation Sky-Catch” tests, conducted jointly by 

S. Navy, Westinghouse and Lockheed, reveal effects of 

; ‘ unch stresses. To prevent damage to POLARIS test missile 
es r\ bmponents, huge overhead assembly catches it in mid-air. 

tt left: Carrying awesome destructive power, Navy’s 

u% WoLaRis is markedly smaller than other U. S. ballistic 

} \fisiles of same range. To enable one submarine to carry 

’ y missiles, Lockheed scientists achieved miniaturi- 


lion breakthroughs on all major systems of POLARIS. 
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Above: Propulsion system firings light up night sky at 
Lockheed Missiles and Space Division’s 4000-acre test 
base in the Santa Cruz mountains, Calif. A nearby vertical 
test stand permits simulated flights of complete POLARIS 
missiles, to test their propulsion and guidance systems. 


Left: Nuclear powered U. S. Navy sub (1) launches 
solid-propellant POLARIS missile which erupts from ocean 
depths (2)—then rockets its way on a ballistic missile tra- 
jectory high above the earth’s atmosphere (3) into outer 
space. Re-entering the atmosphere (4) the re-entry body 
of the PoLaRis plunges earthward to destroy its target (5). 


Launching an IRBM from submerged submarines in- 
volves technological developments that are unique in 
ballistic missile history. 


As POLARIS Missile System manager and prime con- 
tractor, Lockheed’s Missiles and Space Division is 
responsible for coordinating all phases of its design, 
development, and production. 


The PoLaris, a solid-propellant missile, will be oper- 
ational in 1960—ahead of schedule. The operational 
missiles will be delivered to U. S. Navy submarines in 
Lockheed ready-service transporters under controlled 
humidity and temperature conditions. And Lockheed’s 
ACRE~—Automatic Checkout and Readiness Equipment 
—will maintain the constant ready-to-fire capability of 
POLARIS missiles, right up to the moment of launch. 


With an ultimate range of 1500 nautical miles, the 
POLARIS will become one of the most formidable deter- 
rent weapons ever devised. For any significant military 
target on earth will be within range of its nuclear warheads. 


JET TRANSPORTS + JET FIGHTERS «+ JET TRAINERS - COMMERCIAL & MILITARY PROP-JET TRANSPORTS + ROCKETRY 
® BALLISTIC MISSILE RESEARCH & DEVELOPMENT » WEAPON SYSTEM MANAGEMENT + ANTI-SUBMARINE PATROL AIRCRAFT 
1s NUCLEAR-POWERED FLIGHT * ADVANCED ELECTRONICS « AIRBORNE EARLY-WARNING AIRCRAFT + AIRPORT MANAGEMENT 
NUCLEAR REACTOR DESIGN & DEVELOPMENT * GROUND SUPPORT EQUIPMENT +» WORLD-WIDE AIRCRAFT MAINTENANCE 
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AUTOMATIC TRANSLATION 


Computing machines aid language research at Ramo - Wooldridge 


To formulate rules for automatic language translation is 
an enormously subtle and complex project. Yet significant 
progress is being made. During the past year of research 
at Ramo-Wooldridge over 60,000 words of Russian text 
have been translated and analyzed using an electronic 
computer. From the beginning several hundred syntactic 
and semantic rules have been used to remove ambiguities 
that are otherwise present in “word for word” translation. 
Our present computer program for automatic translation 
is a considerable improvement over earlier attempts. 


Apart from the question of translation itself, electronic 
computers are invaluable for language research. The 
expansion of existing knowledge of the rules of language, 
through statistical analysis, is made practical by mecha- 
nized procedures. A clear symbiosis between linguistics 
and computer technology has emerged. 


Automatic translation research is one of many R-W 
activities addressed to problems of communication of 


scientific information. These problems are increasing at 
an accelerating pace. In this area, as in others, scientists 
and engineers find at Ramo-Wooldridge challenging 
career opportunities in fields important to the advance of 
human knowledge. The areas of activity listed below are 
those in which R-W is now engaged and in which open- 
ings also exist: 

Missile electronics systems 

Advanced radio and wireline communications 

Information processing systems 

Anti-submarine warfare 

Air navigation and traffic control 

Analog and digital computers 

Infrared systems 

Electronic reconnaissance and countermeasures 

Basic and applied physical research 


For a copy of our brochure or other information, write to 
Mr. Donald L. Pyke. 


RAMO-WOOLDRIDGE 


P. O. BOX 90534, AIRPORT STATION + LOS ANGELES 45, CALIFORNIA 


a division of Thompson Ramo Wooldridge Inc. 
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Recent Advances in 
Solid Propellant Grain Design 


JEAN A. VANDENKERCKHOVE' 


Université Libre de Bruxelles 
Bruxelles, Belgium 


The author has been Assistant at the Institute of Aeronautics of the Université Libre de Bruxelles 
since 1954. He was attaché to Professor von Karman in 1955 and is currently working with Professor 
Jaumotte, Head of the Institute and one of the Belgian representatives at the Combustion and Pro- 


pulsion Panel of AGARD (NATO). 


He received a Mechanical and Electrical Engineer degree from 


the Université Libre de Bruxelles in 1952 and an M.S. degree from the California Institute of Tech- 
nology in 1954. He has written several papers, especially in the area of solid propellant rockets, and 
is co-author of the book ‘*Rocket Propulsion,’* which, after a first French edition in 1957, is to be 
published in English and in German this year. 


WIDE variety of grain designs has been actually used or 
‘an be considered (1).2_ However all possibilities more 
or less fall into one of the three following categories.* 


End-Burning Grains 

With end-burning grains, also called restricted or cigarette- 
burning grains, the combustion is limited to the transverse 
cross section of the chamber. 

It is possible, with such rockets, to achieve near constant 


thrusts during very long burning times, up to 10 min and 


more, since the burning area remains constant, and the avail- 
able propellant thickness, which is equal to the chamber 
length, can become quite important. However, despite their 
optimum loading density, such end-burning rockets are 
mostly low performance, low thrust motors. Indeed, for long 
durations, it is difficult to protect the chamber hardware from 
heating, and the walls must then become fairly heavy to resist 
pressure. Moreover, inhibiting poses a difficult problem, and 
the center of gravity shifts during combustion. Finally, it is 
often impossible, with available burning rates, to develop 
large thrusts with reasonably limited cross sections. 


Side-Burning Grains 


With side-burning grains, the combustion takes place on a 
straight cylinder formed by the lateral surfaces of the grain, 
and the end surfaces are often protected from burning, or at 
least they do not represent a significant fraction of the burning 
area. 

In view of their small cross sections, light weights and wide 
range of thrust laws obtainable, these side-burning grains are 
used in most cases, even for relatively long burning times of 
the order of 1 min. 

Moreover, the chamber walls can often be protected from 
heating by the propellant itself, which can be free-standing or 
case-bonded. However, suitable passages or port areas must 
be provided for the gas flow within the chamber, and the 
loading density is accordingly reduced. In fact, the latter 
decreases together with burning time. From the above men- 


Received May 19, 1959. 

! Associate Member ARS. 

2 Numbers in parentheses indicate References at end of paper. 

’ Adapted by permission from material in “Rocket Propulsion’’ 
by M. Barrére, A. Jaumotte, B. Fraeys de Veubeke and J. 
\ andenkerckhove, distributed in the U. S. by Van Nostrand Co., 
and published by Elsevier, Amsterdam, 1959. 
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tioned geometrical property, it is seen that the grain design 
‘an be studied in a transverse cross section. Fig. 1 represents 
several possible cross sections; it illustrates the wide possi- 
bilities offered by this category. 

A side-burning grain is called unrestricted when the com- 
bustion takes place on almost all lateral surfaces (Fig. 1, 
vases A, B, C, L, M, N and O). The use of unrestricted 
grains results in severe wall heating and relatively large losses 
due to final grain fragmentation. Modern designers gener- 
ally try to protect the walls by the propellant itself (Fig. 1, 
cases D. E, F, G, H, I, J and K), in which case wall heating 
can be kept low, thus making for very light chambers. 

Side-burning grains are often called cylindrical grains; 
however such charges need not necessarily be cylindrical on 
their inhibited surfaces. 


Noneylindrical Grains 


This last category includes the grains which cannot enter 
one of the two previously described classes. It is character- 
ized by the fact that the generatrices of the burning area 
are not parallel straight lines any more. Fig. 2 illustrates 
three designs representative of this category, the so-called 
slotted. evlindroconical and cruciform grains. The compu- 
tation of the characteristics of such charges is often difficult, 
because it cannot be reduced to a two-dimensional problem. 
However, chamber slenderness is often quite high, and it is 
then possible to simplify the problem by considering, as a first 
approximation, that the grain is composed of several portions 
of straight side-burning surfaces (slightly tapered port), or 
that the grain is a straight side-burning surface on which 
small modifications are added (influence of end surfaces or 
other peculiarities, such as smal] holes drilled for eliminating 
irregular burning). 

It must be noted that the grain characteristics sometimes 
heavily depend on the inhibiting pattern. For instance, the 
cruciform grain with no inhibition at all, is strongly regres- 
sive, whereas it is progressive if all arms are completely in- 
hibited. It has neutral characteristics only when 50 per cent 
of the supporting surfaces are restricted, in which case it is 
rigorously a side-burning grain only at the beginning of the 
combustion, since noneylindrical areas develop afterward near 
the edges of the inhibitor strips. (The arrangement of inhibi- 
tor strips shown in Fig. 2 has been selected to avoid irregular 
burning.) 

On the other hand, the amount of final losses or slivers also 
depends heavily upon grain geometry. Near the end of the 
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Fig. 1 Wide variety of side-burning grain cross sections 


combustion, the burning area is often strongly regressive. 
Indeed, when the major part of the burning surface dis- 
appears, near burnout, a fraction of the propellant remains 
attached to the walls or to the inhibitor, and continues to burn 
at low pressure without producing an appreciable thrust. 
That fraction, which should be kept to 5 per cent or less, must 
be considered as a loss whose importance can be computed by 
geometry (“theoretical unburnt residue’). This phenom- 
enon is known as the tail-off. Since it is difficult to tell how 
much of the slivers will actually burn, variability is thus 
added to the thrust law. 

On the other hand, with some designs, the combustion takes 
place on both sides of the propellant thickness (Fig. 1, cases 


Fig. 2 A: slotted grain, B: cylindroconical grain, C: cruciform 
grain with its classical arrangement of inhibitor strips 


A, B, C, K, L, M, N and O), and final breakup takes place in 
an erratic manner. A loss of about 5 per cent of the web 
thickness is often assumed to account for final fragmentation. 
The actual value, however, depends upon grain geometry, 
the higher the grain compactness, the lower the loss. The 
free slivers are then discharged through the nozzle. The char- 
acteristics of many simple grain designs are discussed with 
some detail in the literature (2 through 8). 


Internal Ballistics 


The proper grain design must give the requested thrust 
program /'(t) which can be neutral, progressive, regressive or 
assume any other pattern, a dual thrust law, for instance. 

Due to the slight dependence of thrust coefficient on expan- 
sion ratio, the thrust is roughly proportional to chamber pres- 
sure which, certain second order effects being neglected, is 
given by the well-known relation 


| 
= (ppaKe*)! (l—n [1] 
with 
K = Aj/A_ 
if the burning rate is expressed by 
r = dy/dt = ap.” {2} 
where 
Pp» = specific mass of the propellant 
c* = its characteristic velocity 
a = constant which depends only on grain temperature 
n = constant called the combustion index which is always 
smaller than 1 
K = ratio of burning area A, to nozzle throat area A, 


The burning rate r is, indeed, the derivative with respect to 
time of the burnt distance y which varies from 0 to w, w being 
the web thickness. 

In grain design, it is desirable to use the burnt distance y as 
the independent variable. This can be easily done by inte- 
gration of Equation [2] in which the pressure is known as a 
function of time p,(t). We thus obtain, successively, the 
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functions y(t), p-(y), and finally, through Equation [1], 
K(y). The problem can thus be reduced to a geometrical 
one. For instance, a progressive thrust law is obtained when 
A, steadily increases with y (9). 

Similarly, a dual thrust law in which the high thrust level 
is developed during a short booster phase and is followed by a 
long sustainer period at moderate thrust, can be obtained with 
a single grain which, at the beginning, provides a large burning 
area, Which after a while, suddenly drops to a much lower 
value (see the wagon-wheel design description in section on 
design charts). 

Equation [1] can be viewed as the transfer function be- 
tween grain geometry and chamber pressure, the term 
1/ 1 — n) which only depends on the propellant nature, acting 
as the amplification factor. 

{t is seen that for large values of the combustion index (n = 
0.7 for instance) small variations in burning area result in very 
important thrust changes, a property which could, conceiv- 
abiy, be used to advantage, although in most cases small 
values of n (n < 0.4) are to be favored for the purpose of re- 
ducing dispersion. 

his brings us to study briefly the problem of propellant 
selection in relation to grain design. For this purpose we 
shall assume that the chamber diameter is known (see sec- 
tion on overall grain dimensions) and that the thrust law is 
approximately neutral. The selection of propellant, chamber 
pressure and grain design is more or less arbitrary. It de- 
pends primarily upon the web thickness w needed to obtain 
the design burning time ¢, through the relation 7 4, = w, in 
which 7 is the mean burning rate. 3 

With side-burning grains which are used in most cases, de- 
sign studies indicate that it is difficult to get web thicknesses 
much over one-half the chamber radius (Fig. 1, cases F and 
G). However, small web thicknesses can be easily ob- 
tained in a wide variety of shapes (Fig. 1, cases H, I, J, K, 
M, N and O), but a decrease in web thickness generally results 
in increased production complexity and a higher sliver frac- 
tion, especially when the combustion takes place on both 
sides of the propellant thickness. It thus appears that, what- 
ever the burning time, it is favorable to select a design which 
protects the walls and whose web is, very roughly, equal to 
one-half the chamber radius, the burning rate being chosen so 
as to correspond to optimum pressure for which the chamber 
weight can be minimized. Two cases will be considered. 


Long Burning Times 


When the preceding “rule of thumb” cannot be readily sat- 
isfied with a suitable propellant, the web thickness can be 
slightly increased. However it seems difficult to exceed w = 
0.7 Re. 

The chamber pressure can also be somewhat lowered at a 
slight cost in specific impulse. If this is not sufficient, the 
burning rate can, in some cases, be reduced without much 
lowering the specific impulse, but often a less energetic formu- 
lation must be used. 

In extreme cases, the chamber radius must be increased, or 
an end-burning geometry must be chosen. In practice, how- 
ever, burning times in excess of 30 sec can be obtained with 
side-burning grains of reasonably low diameter. 


Short Burning Times 


Obtaining short burning times is much simpler, since the 
web can be reduced easily at the cost of increased production 
complexity and higher sliver fraction. 

It must be emphasized, however, that throughout the pre- 
ceding discussion, the terms “long” or “‘short’’ burning times 
do not have an absolute meaning, since the scale effeet must 
he accounted for. Indeed, it is much easier to obtain a long 
burning time with a large motor than with a small one. 
‘herefore, high burning rates must be favored either when the 
hurning time is short or when the rocket is large. 
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Fig. 3 Grain design using two propellants and its thrust law. 

During the first phase, a wagon-wheel design provides the initial 

boost; it is followed by a sustainer period at the end of which a 

layer of fast burning propellant is reached and provides a final 
impulse 


If a complicated thrust law is to be obtained, as with a dual- 
thrust motor, it is often difficult to vary, with a single grain, 
the burning surface within sufficiently wide limits, and two 
propellants having different burning rates can be used (10, 
11, 12). Such a ease is represented on Fig. 3 adapted from 
(10). 

Before leaving the subject of internal ballistics, two impor- 
tant phenomena must be mentioned. 

The first one is the pressure drop along the grain length 
which, for a constant port area, results from the mass addi- 
tion and the corresponding flow acceleration. The maximum 
pressure is reached at the chamber front end, and its value 
‘an be approximately computed by Equation [1] in which 
ratio K is replaced by an effective value 


Kerr = K{l + (T'2/2)J?] 


with 
= A;/Ap 
where 
I. = function of the specific heat ratio alone (for y = 1.25, 
we have’ = 0.658) 
A, = port area at the moment under consideration 


In high performance designs, the loading density is quite high 
together with ratio J (J > 0.65). 

A sizable pressure peak thus takes place right at the begin- 
ning of the burning time when J still has its maximum value. 

This phenomenon can be compensated by reducing slightly 
the initial burning area so as to keep A,,;, constant during 
combustion. 

The second phenomenon consists in an increase in burning 
rate under influence of the flow velocity in the port. It is 
called erosive burning. 

These two phenomena result in a variation in burning rate 
along the grain length which increases tail-off, and, in severe 
cases, may lead to burn-through due to premature contact of 
the flame front with a portion of the chamber walls. 

By tapering the port area so as to increase the passage area 
in proportion to the mass addition in each section, the mean 
flow velocity and the pressure can be kept constant, together 
with the burning rate, along most of the grain length (see 
Fig. 4). 

The loading density can thus be improved significantly with 
short burning times, and a constant web thickness can be 
used. 


Similarity Rules 


In designing a solid propellant rocket, the thrust law (thrust 
vs. time) is generally imposed, but the choice of motor dimen- 
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Fig. 4 Longitudinal section of a case-bonded rocket with tapered port 


sions, chamber pressure and propellant nature is left, within 
limits, to the designer. 

Tn certain cases, however, the motor length, its diameter or 
even the chamber wall thickness can be imposed by ballistical 
or structural considerations. 

The choice of one of the powerplant characteristics, its 
length for instance, influences almost all of the other param- 
eters. 

The similarity rules of solid propellant motors provide the 
means of easily arriving at a first approximation of the char- 
acteristics of a new design by extrapolating the dimensions of 
existing powerplants (13). 

Let us consider a family of geometrically similar solid pro- 
pellant rockets (same propellant and same grain design). 
Each member of the family is characterized by a length / (the 
chamber diameter for instance), and it can be readily shown 
that: 

Burning time t is proportional to /. 

Thrust is proportional to 72. 

Total impulse J is proportional to /*. 

Chamber pressure p, (or its variation as a function of a di- 
mensionless time scale) is independent of 1. 

Loading density of the chamber and the ratio J of the throat 
area to port area. which characterize internal ballistics, are 
independent of 1. 

These similarity rules often allow for the reduction of the 
number of full-scale tests, since much data can be obtained 
with scaled-down models. 


Overall Grain Dimensions 


In most cases, the overall grain dimensions, i.e., diameter 
and length, are approximately equal to the overall chamber 
dimensions, allowance being made for the inhibitors, liners, 
supporting elements and constructional features, such as the 
end closures. 

The grain design is thus directly affected by the motor di- 
mensions. When the latter are imposed, the choice of a suit- 
able grain configuration can be carried out almost straight- 
forwardly. 

However, for system analysis where optimum design has to 
be determined, it has been pointed out by Thackwell among 
others, that the motor dimensions should be kept as flexible 
as possible (14). 

Most of the published works have taken advantage of this 
concept, either by direct scaling of known powerplants by 
means of the above mentioned similarity rules, or in a some- 
what less restrictive manner, by assuming a priori some kind 
of geometrical condition and by allowing the other parameters, 
the chamber pressure or the grain slenderness, to vary within 
limits (14 through 17). 
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A widely used condition has been that the grain web thick- 
ness must always be kept equal to one-half the chamber 
radius. The loading density of the chamber can also be 
chosen a priori. 

Such a view, although it has provided very useful results, is 
somewhat restrictive, because it requires the preliminary 
choice of the propellant kinetic characteristics (burning rat 
vs. pressure). From a theoretical point of view, such an ap- 
proach does not provide a truly optimum design, since th: 
chamber dimensions are determined a posteriori from the 
grain design choice. 

A much more general approach is obtained by minimizing 
the total hardware weight while limiting the assumptions made 
on grain design and propellant performance to a minimum 
(18, 19). 

In fact, this method requires only the choice of the propel- 
lant energetic properties (its characteristic velocity c*) and 
the assumption that a side-burning grain is used (possibly 
with a tapered port area). 

The propellant volume required for obtaining the specified 
total impulse can thus be computed readily. It decreases 
slightly when chamber pressure increases due to the favorable 
influence of expansion ratio on thrust coefficient. 

Similarly, for a given thrust level, the nozzle throat area, 
varies approximately as the reciprocal of chamber pressure. 
It must remain proportional to the initial port area, since the 
value of the parameter J = A,/Ay, which must always be 
smaller or at most equal to unity, can be chosen a priori, 
taking into account the magnitude of the difficulties that could 
be encountered during the development. The ratio J is rep- 
resentative of the empty space in the chamber. In the case of 
a tapered port area, its mean value along the grain length 
must be used. 

For each pressure, the propellant volume and grain port 
area are known, and the inert hardware weight can be given 
as a function of chamber radius alone. Thus an optimum 
chamber radius exists for each pressure. 

When chamber pressure increases, the minimum dead 
weight steadily increases, while the propellant weight de- 
creases as shown by Fig. 5. 

An optimum pressure therefore exists for which the total 
rocket motor weight passes through a minimum. It can be 
shown that this true optimum does not correspond to the 
maximum loading density or to the maximum value of the 
ratio A of the propellant weight to total motor weight. 

This method provides true optimum dimensions which can 
often be adhered to. In some cases, however, this may not be 
actually possible, and it then provides the tool for evaluating 
the penalty resulting from the specifications. For instance, 
an overlong burning time at moderate thrust may lead to 
excessive chamber diameter, whereas a large impulse produced 
in a very short time may result in a poor loading density. 
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Fig. 5 Inert motor weight, propellant weight and minimum 

total motor weight as functions of chamber pressure. The locus 

of the minimum inert weight is the envelope of the curves of the 

inert weights obtained by taking arbitrarily different chamber 

radii. The motor is due to give a 5000 kg thrust during 4 sec. 
Its ratio / is taken equal to 0.5 


Fortunately, it is seen that, in most cases, the optimum is 
not very pronounced, especially when construction is light 
and burning time is long. 

Finally, it must be emphasized that the preceding method 
maximizes the ratio of total impulse to total motor weight. 
In many cases, however, the burnout velocity of the missile 
is to be optimized, and the method must be slightly modified, 
taking into account the fact that the dead weight must still 
be minimized. 


Fundamentals of Grain Design 


A solid propellant grain burns in such a manner that any 
point on its surface moves perpendicularly to the tangent 
plane at that point. Therefore, in a cross section, a cusp on 
the perimeter, convex toward the gas phase, remains a cusp 
as burning progresses, while a cusp concave toward the gas 
phase becomes an arc of a circle, with its center at the original 
cusp. As combustion proceeds, new convex cusps can be 
formed, as the result of the interaction between two smooth 
curves (see Fig. 6, case A); less frequently, convex cusps also 
disappear. 

In most cases, as pointed out earlier, the grain characteristics 
can be studied in a cross section perpendicular to the grain 
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axis. In such a two-dimensional study, the most interesting 
parameters are: 

The length S of the burning perimeter, which is a function 
of the burnt distance y. (For a side-burning design, we have 
SL = Ap, L being the grain length.) 

The port area A, from which the ratio J can be readily de- 
duced. 

The geometrical amount of slivers, which can be defined 
as the final grain cross-section area A, when the burning area 
begins to drop sharply, near burnout. 

Epstein has established two useful properties (20): 

1 The length S,,) of a portion of perimeter AB, free of 
convex cusps, can be expressed as a function of its initial length 
Si =o, of the burnt distance y, and of the angles and gz 
made by the end tangents with a direction of reference (see 
Fig. 6 case A), by the relationship 


= Siy=0 + — GB) 


If the perimeter is a closed curved, we thus have 


¢s = + 20 


and the rate of increase in Jength of a closed perimeter free of 
convex cusps does not depend on its length or form. 

As an illustration of the usefulness of this property, it can 
readily be shown that, during combustion, the increase in 
burning perimeter of the multiperforated grain, represented 
by case K of Fig. 1, is equal to 142w, w being the web thick- 
ness. 

It appears, therefore, that an internal burning grain without 
cusps is progressive, the larger the number of perforations, 
the more pronounced the progressivity. Likewise, an un- 
inhibited and cusp-free outer surface has a regressive influ- 
ence; hence a tubular grain is neutral. 

2 The length S,,) of a portion of perimeter AB containing 
one convex cusp, can be expressed as a function of the 


py 


Fig. 6 Evolution of a portion of burning perimeter as combustion 


proceeds 
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parameters S,,- 0), ¥, ¢4, and of the cusp apex angle @ 
(see Fig. 6, case B), by the relationship 


Sy = 8 ( ) a 
w = Sy =o) + — Ga) + 2 2 y 


As pointed out by Epstein, it is thus possible to compute a 
suitable design from the a priori choice of a few parameters 
(number of axis of symmetry of the perimeter). 

For instance, neutral characteristics can be obtained with 
an internal burning design having n identical convex cusps or 
star points (Fig. 1, cases F and G), if, during combustion, 
can be kept constant and equal to the value 6 given by 


3 

Nim 

a 


+ cot , 


This relationship yields: 


6/2 (deg) 


28.21 
31.12 
33.53 
35.55 
37.30 
38 .83 
40.20 
42.52 


3 


— 


Design Charts 


In practice, the grain designer has at his disposal a series of 
charts on which the characteristics of numerous relatively 
simple designs are represented in a wide range of values of the 
different parameters. Very often, one of these basic geome- 
tries, generally built with segments of straight lines and with 
radial ares, is well suited and can be readily selected, with 
possibly minor changes. 


The computation of such charts is a necessary, straight- 
forward but tedious geometrical job, the description of which 
is outside the scope of this review. Almost every group work- 
ing on solid propellants has developed its own satisfactory 
method around the same basic equations. Some methods, 
however, are believed to be more convenient than others for 
the purpose of providing the maximum amount of information 
in the minimum number of charts covering the widest ranye 
of the design parameters. A dimensionless presentation s, 
of course, to be favored for this purpose. 

With such charts, the grain designer can investigate rapid'y 
and systematically the complete range of possible shapes, aid 
is thus in a position to optimize the design. 

The most important case of internal burning, single-per- 
forated grains (internal burning star and wagon-wheel co:- 
figurations) has been extensively studied in (5 and 21 through 
25). 

As an illustration, Fig. 7 represents the charts giving tle 
characteristics S/l, A,/l? and A ,/I? of the eight-pointed inte ’- 
nal burning star, as functions of the reduced burnt distance 
(f+ y)/L. 

Radius / of the basic perforation without fillet has been 
chosen as unity, and the other geometrical parameters n, 
f, 9 and w are defined on Fig. 7. 

Depending upon the value of angle 6, the combustion cai 
be progressive, neutral or regressive during the first phasc, 
during which the straight star points have not yet disappeared. 

During the second phase, the star points have disappeared, 
and after a while the grain becomes always progressive. 

As an example of the information provided by the charts, 
Fig. 8 represents the characteristics of one particular eight- 
pointed wagon-wheel design which gives a dual-thrust law. 
The geometrical parameters are defined on the left of the fig- 
ure. 

It must be noted that the second phases of the wagon-whee! 
and internal burning star are identical, and that the second 
phase progressivity can be compensated by allowing the grain 
ends to burn unrestricted, and by tapering the grain length 


fty 
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Fig.7 Geometrical definitions of the internal burning star design and charts representing the characteristics of the eight-pointed star in 
a dimensionless form. Burning perimeter during the first and second phases, port area during the first phase and sliver area during the 
second phase, as functions of the dimensionless burnt distance. From (24) 
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at the outer diameters (Fig. 4)... Both methods result in a de- 
creasing grain length, as combustion proceeds. 

It must be emphasized that, in addition to suitable geo- 
metrical characteristics, the choice of a design must also take 
into account several other important factors: 

| In many applications, the reproducibility of the thrust 
law is of prime importance, and the tail-off period must be re- 
duced toa minimum. This could be achieved with an inter- 
nul burning star by filling with an inert material the dark zone 
representing A, on the left of Fig. 7. Other more elaborate 
designs can conceivably be developed, for instance by means 
o! grains made with two propellants having slightly different 
burning rates. 

In the case of ballistic missiles, however, thrust termina- 
tion methods must be resorted to. 

2 Some designs are more subject than others to irregular 
or resonance burning as discussed for instance by Green (26). 

3 Finally, the actual realization of cores or mandrels repre- 
sinting the chamber free-volume is sometimes difficult and 
limits the choice of grain shapes. 


Mechanical Behavior of Grains 


In many cases, due to the mechanical and therinal stresses, 
the grain resistance is critical and can be improved through 
proper design. The physical properties of the propellant 
niust thus be carefully taken into account in the choice of the 
grain geometry. 

As a rule, hard propellants are used for unrestricted and 
free-standing grains which are locally supported. In these 
cases, the allowable grain deformation under acceleration and 
internal pressure difference has to remain sufficiently sinall. 
lor instance, Wimpress has shown that in the simple case of a 
tubular unrestricted grain, the maximum possible value of 
ratio J increases together with the modulus of elasticity of the 
propellant (3). 

On the other hand, soft propellants are needed for case- 
honded grains for which a good elongation is imperative in 
order to accommodate the chamber expansion. The strains 
are then imposed, and the stresses must be kept within the 
possibilities of the propellant which depend, among other fac- 
tors, on temperature and rate of the strain. 

The discussion of the rheological properties of the propel- 
lants is outside the scope of this survey, and we shall only sum- 
marize here the few fundamental papers which have been pub- 
lished on this subject and which are of direct interest for grain 
design. 

It must be noted, however, that these published works are 
more or less oversimplified, since they assume a perfectly elas- 
tic behavior whereas the propellant can exhibit viscoelastic 
properties. Even so they are limited, by necessity, to simple 
two-dimensional models. However, these works provide the 
grain designer with very useful qualitative and even quanti- 
tative indications for stress analysis. 


Thermal Stresses 


An important environmental effect on a rocket is the ther- 
mal stress produced in the propellant by the thermal cycling 
during storage. This problem has been excellently studied 
by Geckler and his collaborators, in the case of an infinitely 
long tubular grain (27, 28, 29). 

These authors have shown that the most severe stresses oc- 
cur with case-bonded grains after a sufficiently long time, 
when thermal equilibrium is achieved. In this particular 
case, the outer tangential strain of the propellant must be 
equal to the differential expansion between the case and the 
propellant. Therefore, the pressure applied by the chamber 
on the propellant can be computed, and the stresses can be 
expressed by means of Lamé’s equations. The maximum 
tensile stress is seen to be the tangential stress a» at the inner 
radius after a cooling period. It can be expressed in a dimen- 
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Fig. 8 Geometrical definitions of the wagon-wheel design an 
characteristics of one particular design of this family 
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Fig. 9 Maximum thermal stress at the inner radius of a case- 


bonded tubular grain for several values of the ratio Ri/R, 
(Poisson’s ratio v = 0.4) 


sionless form by the relationship 


Ea,(T.— T) — (Ri/R.)? 
20 = ») 
where 

v = Poisson's ratio 
E = modulus of elasticity of the propellant 
a, = coefficient of linear thermal expansion of the propellant 
a. = coefficient of linear thermal expansion of the chamber 


To = stress-free temperature 
T. = final temperature 

R; = inner radius of the grain 
R,. = outer radius of the grain 


Fig. 9 represents the variation of 


ad 
Ea,(T. To) 


as a function of a,./a, for different values of R;/R, in the case 
y = 0.4. 
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For our purpose, the important conclusion is that the maxi- 
mum tensile stress increases together with the ratio R,/R; 
which is representative of the web thickness. 

The same conclusion holds true with free-standing grains, 
but the stresses are then considerably less severe, and their 
maximum occurs early during the heating period. 
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Fig. 10 Maximum tangential stress at the inner radius as a 
function of the ratio R;/R.. The curves corresponding to case- 
bonded grains have been computed with » = 0.4 
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WEB THICKNESS 
FILLET RADIUS 
Fig. 11 Design curves showing the influence of the web thick- 
ness and fillet radius on the stress concentration factor at the 
inner surface, for six- and eight-slotted configurations. Adapted 
from (30) 
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Mechanical Stresses 


During combustion, stresses are also induced, even if the 
pressure difference between the front and nozzle ends and the 
axial acceleration are neglected. 

Lamé’s equations can again be used for computing the maxi- 
mum tensile stress, which is the tangential stress at the inner 
radius, except in the simple hydrostatic case of a free-standiny 
grain whose inner and outer surfaces are subjected to the samo. 
pressure (pe = pi) (2). 

For an internal burning free-standing grain whose outer 
annulus between the propellant and the wall is not pressur- 
ized by the combustion gases (p, = 0), the maximum tangen- 
tial stress varies from p; to infinity when the web increases. 

Such a high stressing cannot be tolerated, and it is thus cur- 
rent practice to allow for gas passage into the outer annulu~ 
in order to provide external support for the grain and preven! 
its failure. 

A more interesting fact, for our purpose, is that for case- 
bonded grains, tensile stresses only occur when the web thick- 
ness becomes larger than approximately one-half the chambei 
radius. Above that value, the maximum tensile stress in- 
creases fairly rapidly, and elasticity of the chamber wall: 
tends to increase its value. Fig. 10 represents the dimension- 
less maximum tangential stress o¢/p; as a function of the ratio 
R;/R, in the cases discussed above (p; is the effective pressure 
in the perforation). 


Stress Concentrations 


The preceding discussion has been limited to hollow, circu- 
lar cylinders. Actual grains, however, exhibit much more 
complex shapes, which are associated with stress concentra- 
tions. 

For instance, the internal burning grains with their complex 
perforations are subject to the inherent weakness that stress 
concentrations occurring around the small angles of the star 
points tend to aggravate any susceptibility of the propellant 
to mechanical failure. The actual stress can easily exceed 
two or three times that which occurs in a simple thick walled 
cylinder. The grain must then be designed with a configura- 
tion that maintains the stress concentration at a mini- 
mum possible level. Using photoelastic methods, Ordahl 
and Williams have investigated the quantitative relationships 
between the various parameters of the grain geometry and the 
elastic stress distribution (30). Figure 11 represents some of 
their results. 

On this figure, the stress concentration factor K is defined 
as 

— oR 
where og and oz are the actual tangential and radial stresses, 
and o»° and o,° are the stresses which would occur in the 
equivalent tubular grain. Within the range of validity of the 
experimental method, these authors have shown that: 

1 The stress concentration increases as the ratio of the 
web thickness to fillet radius increases. (It is thus impera- 
tive to use sufficiently large fillets.) 

2 The stress concentration decreases as the web fraction 
increases. 

3 Fora given web thickness and fillet radius, the stress con- 
centration decreases as the number of star points is increased. 

4 The stress concentration decreases as the star point 
angle increases. 

5 The stress concentration decreases as the width of the 
slot increases 

Reference (30) only considers steady-state stresses due to 
pressure or differential thermal expansion of case-bonded 
grains. Ina recent paper, Williams has extended his inves- 
tigation to the case of transient thermal stresses for which 
two concentrations must be considered: A geometric concen- 
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tration factor K, representing the difference in planar rigidity 
of the star from a thick walled cylinder of the same web 
fraction, and a temperature factor K, due to the concentration 
of heat flux at the star points (31). 

The actual concentration factor K is the product K,A,;. The 
geometric concentration factor K, is found to be very similar 
to the one measured with pressure loading, while the tempera- 
ture factor K;, is always smaller than unity. 

his favorable influence is more pronounced for small web 
fractions and small fillet radii. It decreases slightly when the 
number of star points increases, and appears to be the result 
of the relaxation (negative concentration) of heat flux away 
from the star points toward the regions of higher thermal iner- 
tia. 
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Explorer Rocket Research Program 


G. ROBILLARD' 


California Institute of Technology 
Pasadena, Calif. 


Nine Jupiter-C missiles have been launched at Cape Canaveral since September 1956, the first 
three testing the propulsion system, airframe, guidance components and re-entry conditions of 
the nose cone, the remaining six were used as carriers of the Explorer satellites I through VI. This 
was a joint program of Jet Propulsion Laboratory and Army Ballistic Missile Agency with consider- 
able assistance from the University of Iowa, Air Force Cambridge Research Center, the Vanguard 
tracking teams and many others. The Jupiter-C is a four-stage vehicle burning solid propellant. 
The first stage is a modified Redstone missile, developing a thrust of 83,000 Ib. The second stage 
is powered by an annular ring of 11 solid propellant motors, patterned after the Sergeant 
missile. The third stage consists of three motors, identical to the second, fitted inside the ring 
of second-stage motors. The fourth stage is a single motor, identical to the motors of the other 
high speed stages. The satellite payloads of the Explorer series were contained in a shell 80-in. 
long by 6 in. in diameter, which, together with its package of instruments for Explorer I, weighed 
30.8 lb. The experiments were designed especially to study and measure cosmic radiation in great 
detail. Explorer IV, the latest one to attain a successful orbit, achieved a perigee altitude of 160 


miles and an apogee altitude of 1390 miles. 


INCE September of 1956, nine Jupiter-C missiles have 

been launched from the firing pad at Cape Canaveral. 
The first Jupiter-C firing tested the propulsion system, air- 
frame and guidance components of the missile, and the second 
and third firings tested a model of the Jupiter nose cone under 
realistic re-entry conditions. The remaining six Jupiter-C 
missiles were used as the launching vehicles for Explorer 
satellites I through VI (Fig. 1). Of the six satellite firings, 
Explorers I, III and IV achieved satisfactory orbits. 

The Jupiter-C missile was designed and developed as a 
joint program under the technical direction of the Jet Pro- 
pulsion Laboratory and the Army Ballistic Missile Agency. 
The Jet Propulsion Laboratory developed the three high 
speed stages, and the Army Ballistic Missile Agency handled 
the development, construction and operation of the first stage 
booster rocket and the guidance system. Many other 
organizations have contributed to the success of the Explorer 
satellite program, most notably the State University of 
Iowa, the Air Force Cambridge Research Center, and the 
satellite tracking teams of the Vanguard organization. 


Design of the Jupiter-C 


The Jupiter-C missile is a four-stage vehicle. The first 
stage consists of a modified Redstone rocket. The three 
high speed stages are comprised of clusters of solid propellant 
rockets placed in a spinning tub which is mounted in the nose 
of the first stage. 

The standard Redstone missile delivers a thrust of 75,000 
lb and employs alcohol as the fuel and liquid oxygen as the 
oxidizer. For the Explorer missions, the propellant tanks 
were enlarged in order to permit an increase in burning time 
to 155 sec, and the alcohol was replaced by another fuel, 
called hydyne. This latter change resulted in an increase 
in thrust of about 8000 lb, giving a total thrust of 83,000 lb. 
The fuel change was effected without major modification of 
the engine. 

The last three stages of the missile are spun in order to 
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minimize dispersion of the payload due to any thrust  mis- 
alignment of the unguided solid propellant rocket clusters. 
The tub is spun by two battery-powered electric motors 
mounted in the top of the instrument compartment, which is 
emplaced between the tanks of the first-stage rocket and the 
spinning tub and contains the guidance and control equip- 
ment for the missile. This compartment also holds a spatial 
attitude control system for horizontal alignment of the high 
speed stages before the ignition of the second stage 

The second stage of the Jupiter-C is composed of an 
annular ring of 11 solid propellant motors (Fig. 2), each of 
which is a 6-in. OD scale model of an early version of the 
powerplant of the Sergeant missile. All of the solid pro- 
pellant motors of the three high speed stages are identical in 
configuration. This motor was chosen for this task because 
it was a reasonably high performance rocket and was known 
to be extremely reliable, having undergone over 300 static 
tests. The 11-motor ring of the second stage is mounted in 
bulkheads inside the spinning tub and is held accurately in 
position by launching lugs. Extremely accurate positioning 
and balance of the three upper stages and the payload are 
required in order to prevent vibration caused by the spinning 
tub and to minimize dispersion of the stages in launching 
and firing. Because of the centrifugal forces on the spinning 
second-stage motors, it was considered necessary to test 
these motors under simulated flight conditions on the ground. 
The radial load amounts to 180 g at the full spin rate of 700 
rpm. Pairs of motors were spun and fired in a static-test 
jig. The motors operated well, but it was necessary to make 
a small change in the nozzle design in order to prevent the 
cantilevered nozzles from deflecting outward under the spin 
loads. 

The third stage (Fig. 3) consists of three motors, identical 
to those in the second stage, which are connected by bulkheads 
in a fashion similar to that used for the second-stage instal- 
Jation. The three-motor cluster fits inside of the ring of 
motors comprising the second stage and is held in place by 
launching lugs. 

A conical structure attached to the top of the third stage 
contains the timer and batteries for firing the third and fourth 
stages and serves to support the single fourth-stage motor. 
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Fig. 1 Cutaway views: Explorer I and Explorer III 


This motor is identical in configuration to those of the other 
high speed stages, except in the manner in which it is sup- 
ported and in its nozzle exit diameter which is somewhat 
smaller to allow adequate launching clearance. 

By the time permission was obtained to fire Explorer I, 
sufficient tests had been made of a propellant having improved 
performance over that used in Sergeant to permit its use in 
the fourth-stage motor. The fourth-stage rocket places the 
lowest requirement on the physical properties of propellant, 
in that the rocket is spun on its own axis; hence the radial 
loads are a minimum. Previous to the firing of Explorers 
IV and V, this propellant was qualified by extensive static 
spin-testing for use in the third stage, and an even more 
energetic propellant had qualified for use in the fourth stage. 
These propellant changes, together with a decrease in the 
weight of some of the inert parts by design refinements, ac- 
count for the increased payload weights of the latter satellites. 
Some potential for increased payload using this missile con- 
figuration still exists, but the return in terms of payload 
weight is rapidly diminishing with respect to the amount of 
development effort required. 

The satellite payload (Fig. 1) is rigidly attached to the 
head end of the fourth-stage rocket motor through an in- 
sulating ring which forms the antenna gap for one of the two 
payload transmitters. The second transmitter operates 
across a similar gap between the cylindrical portion of the 
payload and its nose cone. 

The entire high speed assembly, including the payload, is 
very carefully aligned and balanced in its tub before the as- 
sembly is mounted on the Redstone booster. 

The design of the Jupiter-C was carried out to afford a 
vehicle which would provide maximum reliability with mini- 
mum development effort. For this reason, proved compo- 
nents were used wherever possible, and in order to assure reli- 
ability, redundant systems were used in many places to pre- 
vent system failure caused by malfunction. The ignition 
system for the high speed stages is one example of such re- 
jundaney. The problem of reliable ignition at high altitude 
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was resolved in three ways: 

1 An igniter was designed which would fire the motors 
reliably in a laboratory vacuum. 

2 The rocket motors were sealed so that they would hold 
atmospheric pressure. 

3 The igniter was sealed in an inner container capable of 
holding atmospheric pressure if the motor seal should leak. 


Firing Sequence 


Prior to takeoff, the tub is spun up to 550 rpm and held at 
this speed. About 70 sec after takeoff, a governor controlled 
by the missile programmer increases the rotational speed 
gradually up to 750 rpm at a time about 20 sec before first- 
stage burnout. This procedure eliminates resonance be- 
tween the spin frequency and the natural bending frequen- 
cies of the missile. These bending frequencies increase as 
the propellants of the first stage are consumed; thus the spin 
rate is increased gradually, so that the spin frequency always 
remains between those of two of the natural modes of vibra- 
tion of the missile. 

The missile takes off vertically under its thrust of 83,000 
lb. During the 155-sec burning time of the first stage, the 
vehicle is tilted into a trajectory which is approximately 40 


Fig. 2 Second-stage cluster of solid propellant motors 


Fig. 3 Third-stage motor cluster 
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Fig. 4 Explorer launching 


deg inclined to the horizon at burnout. A depletion cutoff 
technique is used to maximize the impulse delivered by the 
Redstone motor. Contacts are placed in the discharge lines 
from the fuel and oxygen pumps; whichever of these contacts 
first senses a pressure drop triggers a relay which closes both 
propellant main valves to the combustion chamber. Five 
sec after burnout, a timer activates the six explosive bolts 
which secure the guidance compartment to the tank section of 
the first-stage motor. Coil springs wrapped around these 
bolts exert a gentle push on the guidance compartment and 
separate it cleanly from the booster. The failure of Explorer 
V to achieve orbit is attributed to a malfunction at this point 
in the launching sequence. In this one instance, the Red- 
stone booster was still exerting some thrust after separation; 
hence the separated booster section was able to overtake and 
ram the guidance compartment. This impact caused the 
gyros in the guidance compartment to tumble; consequently, 
the reference planes were misoriented during the ensuing 
portions of the launching sequence. 

From the point of separation, the two portions of the missile 
coast through a vacuum trajectory until the apex is attained 
after about 404 sec after takeoff. During the free-coasting 
period between separation and apex, the spatial attitude 
control system in the guidance compartment aligns the high 
speed stages into an exactly parallel position with respect to 
the local horizon. This reorientation uses the gyros of the 
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guidance compartment as reference and is accomplished by 
means of small compressed air nozzles located at the base 
of the guidance compartment. 

Because of the relatively crude cutoff technique based on 
propellant depletion used in the first stage, it is impossible to 
predict the exact time when apex will be attained. Only by 
igniting the high speed stages at the moment of apex and by 
accurate alignment of the stages in the horizontal direction 
can orbital flight be assured. Three independent methods 
are employed to determine the time of apex. First, a radar 
tracking plot is used to determine the time at which apex 
will be reached; second, an accelerometer in the missile relays 
by telemetry the buildup of velocity in the first stage; an: 
third, standard Doppler tracking furnishes the same inform: - 
tion. The data obtained from these sources are fed into 
small computer that provides an average estimate of apex 
time, but allows for the respective quality of the three types 
of data. This value is then used to set a timing device whic): 
sends a radio signal to the missile to fire the second stag: 
The computation is accomplished in the 4-min interval be- 
tween separation and apex. 

Each of the three high speed stages has a burning time o/ 
about 6 sec, and 2 sec of coasting time are allowed between th: 
firing of successive stages to permit complete thrust decay. 
The second-stage firing signal is received by the command 
transmitter in the instrument compartment, and the stag: 
is fired by an electrical signal from batteries in the guidance 
compartment. The pressure rise in the second-stage motors 
signals the start of a timer housed in the conical support for 
the fourth stage, which, 8 see later, fires the third stage and, s 
sec after that, the fourth stage. As each stage fires, it 
breaks shear pins holding it to the preceding assembly and 
leaves that stage behind to fall back to Earth. 

The failure of Explorer II to achieve an orbit has been 
traced to a structural failure in the igniter support for the 
fourth stage. The light plastic cone which supported the 
igniter from the nozzle end of the rocket motor apparently 
failed under the combined spin, inertial and vibration loads. 
thereby allowing the igniter to slip out through the nozzlc 
into the conical structure supporting the fourth stage. 
Strengthening the igniter support prevented this malfunc- 
tion in subsequent rounds. 

The payload of Explorer VI contained a fifth small solid 
propellant rocket motor stage oriented in such a manner that 
its nozzle exhausted through the nose cone of the payload. 
This motor was designed to be fired by a timer at apogee after 
the first half-revolution of the satellite around the Earth. 
At this time the motor would have been oriented in such a 
way that it could impart additional velocity to the payload 
along its path, thereby raising the perigee height well above 
its initial value. This ‘“‘kick-in-the-apogee” technique is an 
economical way of achieving a high orbit, because it eliminates 
the necessity of raising all of the missile hardware to a high 
altitude. 


Payloads and Orbits 


The entire orbiting assembly in Explorer I, that is, the 
empty rocket motor shell from the fourth stage plus the in- 
strument compartment, weighs 30.8 lb. The instrument 
compartment alone weighs 18.0 lb. The total assembly is 80- 
in. long and 6 in. in diameter. This satellite achieved an 
initial orbit having a perigee altitude of 225 miles and an 
apogee altitude of 1594 miles, corresponding to a period of 
revolution of 114.78 min. Tracking data showed that the 
satellite was injected into orbit at an angle only about 0.8 
deg from the horizontal, although a suitable orbit would have 
been achieved even if this error were as great as 4 deg. 

The shell of the instrument compartment was constructed 
of 410 stainless steel, 0.025-in. thick. The outside surface 
was sandblasted and then coated in stripes with aluminum 
oxide (Rokide A). The region of the stagnation point of the 
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nose cone was completely covered to guard against aero- 
dynamic heating during launch. The payload surface was 
treated in this manner so that the heat absorbed by radiation 
from the sun and by reflected radiation from the Earth plus 
the heat generated internal to the payload would balance the 
heat radiated by the payload to space. By this means, it 
was possible to maintain the electronic components within 
u temperature range in which they could function properly. 
Because this technique could not be expected to maintain a 
constant temperature of the satellite skin as the satellite 
moved into and out of the Earth’s shadow, the sensitive 
components were thermally insulated from the skin to 
smooth out the temperature fluctuations. The success of 
this technique is indicated by the fact that the maximum and 
minimum temperatures recorded within the instrument 
compartment were 104 and 43 F, respectively. 

Explorer I carried two minimal weight transistorized radio 
transmitters. The low power transmitter radiated 10 mw of 
energy, Whereas the second transmitter radiated 60 mw. 
Because of its higher power requirement, the 60-mw unit 
transmitted for only about 2 weeks; however, the low power 
transmitter was tracked for several months after launching. 

The primary mission of the transmitters was to provide a 
signal for tracking the satellite. The high power trans- 
niitter could be received by sensitive amateur receivers, but 
the low power transmitter could be tracked only by the more 
sophisticated narrow-band microlock and minitrack re- 
ceivers. In addition to providing a tracking signal, the 
transmitters were used to telemeter scientific data from the 
satellite; these data included measurements of temperature, 
micrometeorite impacts and cosmic-ray intensity. 

Temperatures were measured at four points: One directly 
under the stagnation point of the nose cone, two others at 
two locations on the outer skin, and another directly within 
the electronic package itself. 

Micrometeorite impacts were counted by two devices. A 
sensitive microphone attached to the inside of the payload 
shell, by means of an amplifier and a scale-of-two circuit, 
caused a frequency shift of a subcarrier oscillator each time 
an impact occurred anywhere on the satellite. In addition 
to the microphone, a micrometeorite erosion gage consisting 
of 11 grids of fine brittle wire was mounted just forward of the 
nozzle of the last stage rocket motor. These grids were wired 
in parallel; each time a grid was broken by micrometeorite 
impact the resistance of the circuit changed, and the change 
was reflected in the frequency of a subcarrier oscillator. 

The intensity of cosmic radiation was detected by a Geiger- 
Miieller tube. This instrument was able to count only the 
number of impingements within its sensitivity range and was 
unable to differentiate between particles of various energy 
levels. Inasmuch as the cosmic-ray intensity had been pre- 
dicted to be of the order of 30 to 40 counts per sec, the range 
of the instrumentation was set to bracket this range. The 
field of cosmic radiation encountered at altitudes greater 
than 600 miles, however, was so intense as to saturate the 
instrument completely. The only conclusion which could 
be reached was that the intensity was at least 1000 times 
greater than expected. 

With Explorer I, telemetering data could be gathered only 
when the satellite was over one of the ground receiving sta- 
tions. These stations were located in Florida, Nigeria, Singa- 
pore and California, and at a number of places along the 65 
deg west meridian. Because of the relatively small number 
of receiving stations, much of the orbit was not observed, 
and for this reason much of the telemetered data was lost. 
In order to eliminate this difficulty, Explorer III carried a 
miniaturized tape recorder which stored the information 
acquired throughout one complete orbit and reported it on 
command when the satellite passed over a suitably equipped 
receiving station. Playback was accomplished at high 
speed, so that approximately 2 hr of data could be trans- 
mitted in 5 sec. Between interrogations, the high power 
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transmitter was turned off to conserve power. Only the 
cosmic-ray information was transmitted in this fashion; the 
other experiments, similar to those in Explorer I, were tele- 
metered continuously by the low power transmitter. 

Explorer III attained an orbit having a perigee altitude of 
120 miles and an apogee altitude of 1740 miles, with a period 
of 115.8 min. The orbit was inclined about 333 deg to the 
equator. The greater ellipticity of this orbit was caused 
by a significant elevation error at injection. In retrospect, 
this represents a particularly favorable malfunction, since the 
high apogee orbit has made it possible to obtain exceedingly 
interesting data from the cosmic-ray experiment. 

The instrumentation in Explorer IV was designed specifi- 
cally to study the cosmic radiation in greater detail. Im- 
provements in the performance of the high speed stages for 
this firing permitted the payload weight to be increased by 7 
lb. By eliminating all of the other experiments, including 
the tape recorder, it was possible to package four instruments 
for measuring cosmic-ray data in this satellite. Two scin- 
tillation counters were carried in addition to two Geiger- 
Miieller tubes similar to those carried in the earlier satellites. 
One of the Geiger tubes is shielded with lead to eliminate 
data below fixed energy levels. The other tube is set to 
respond to any particle which penetrates the shell of the 
satellite and the casing of the Geiger tube. One of the scin- 
tillation counters records a pulse each time a particle passes 
through it. The other can measure the energy level of the 
particles which it encounters. Both scintillation counters 
are shielded by lead, except for one window in each where 
the shielding has been made as light as possible and a hole 
has been cut in the shell of the satellite. 

Explorer IV was fired on a much more northerly orbit than 
its predecessors. It achieved a perigee altitude of 160 miles 
and an apogee altitude of 1390 miles; the period of the satel- 
lite was 110.3 min. Its orbit inclined 50.4 deg to the equator. 
By this means, it was able to investigate a far greater volume of 
space with respect both to latitude and to altitude than were 
the earlier vehicles. 


Scientific Data 


The justification for a project of this magnitude, of course, 
is the scientific data which are obtained. The Explorer 
project has, in fact, acquired a great deal of very interesting 
data; so much so that it may be several years before this 
extensive accumulation is completely analyzed. However, 
sufficient preliminary information is now available to sketch 
a rough picture of space in the vicinity of the Earth. 

Although the temperature measurements taken inside the 
satellite were not intended to measure any fundamental 
quality of the space environment, they served to prove the 
success of the scheme which had been employed to control 
the satellite’s temperature. We now have demonstrated 
that, at least insofar as temperature is concerned, we shall 
be able to provide a habitable environment for man in space. 

The rate of decay of the orbit, particularly that of Explorer 
III, which lasted in orbit for about three months, has pro- 
vided a measure of the density of the upper regions of the 
atmosphere. Explorer III completed approximately 1260 
revolutions before its disappearance. Its orbit ranged from 
a perigee of 116.86 miles and an apogee of 1740.77 miles on its 
first pass to 107.2 miles and 650.2 miles, respectively, about 
two weeks before it vanished. The rate of decay of apogee 
was 11.36 miles per day between March 26 (the date it was 
launched) and April 11, 11.95 miles per day between April 11 
and 25, 14.99 miles per day between April 25 and May 31, 
14.16 miles per day from then until June 10, and 16.9 miles 
per day thereafter. The exact significance of this somewhat 
erratic path in terms of atmospheric density is still being 
analyzed. 

The measurements of micrometeorite density have sub- 
stantially confirmed previous estimates based on such in- 
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direct. evidence as the data from high altitude sounding 
rockets, and the measurement of ocean floor sediment. The 
results of a partial analysis of the data from Explorer I by the 
Air Force Cambridge Research Center, which is responsible 
for this experiment, show for example that the average flow 
of particles in the 4-4 range is 30 particles per yd? per hr, 
whereas for a larger particle of about 10u in diameter, the 
rate is three particles per yd? per hr. For the spaceship de- 
signer, these numbers offer no problem. However, they take 
on an interesting meaning when the surface considered is the 
surface of the Earth. In these terms, we find that cosmic 
dust accounts for an increase in the mass of the Earth of 2000 
tons per year. 

It is interesting to note that the Russians carried on a 
similar experiment aboard the Sputniks. Surprisingly, they 
have reported a density of particles which is approximately 
1000 times as great as that measured by Explorer. This 
result is difficult to understand not only because it disagrees 
with the previous estimates, but also because the increase in 
mass of the Earth at this high rate would be expected to have 
had a noticeable effect on its rotation speed over the centuries. 

In contrast to the micrometeorite measurements which 
appear to verify preliminary estimates, the measurements of 
cosmic-ray intensity gave some surprising new _ results. 
The first two Explorers indicated that the level of cosmic radi- 
ation was about as expected up to 600 miles altitude in the 
region near the equator. From this altitude up to the highest 
altitude reached by the satellites, the cosmic-ray intensity 
increased at a much higher rate than had been expected. 
At the farthest reaches of the orbits, the count was more 
than 1000 times the expected value, and the rate of increase 
showed no tendency to level off. 

Preliminary analysis of the data from Explorer IV shows 
that the number of cosmic-ray particles passing through the 
unshielded Geiger tube is of the order of 36,000 per cm? per 
sec, from all directions. At least 60 per cent of these particles 
have enough energy to penetrate the lead shield around the 
other Geiger tube. In addition to these very energetic par- 
ticles, there is an even greater density of particles of low 
energy. Atthe northern and southern limits of the trajectory, 
these low energy particles predominate. Here, the flux of 
particles of low energy at an altitude of about 1400 miles is 
about one million per cm? per sec, from all directions. As 
the satellite rotates around its axis, the density measure- 
ments change. On the basis of this result, it appears that 
more than 80 per cent of these particles come from one 
favored direction. 

On the basis of the preliminary estimates, it is not possible 
to separate altitude effects from latitude effects, and, because 
of the magnetic field, both of these effects should cause changes 
in cosmic-ray intensity. The data of the first two weeks 
indicated an increasing intensity with altitude at all latitudes. 
The radiation has been found to vary both in quantity and 
in quality with latitude and altitude. A low energy com- 


ponent similar to that suggested on the basis of the results of 
the first two Explorer satellites is usually present, whereas 
high energy particles become more and more predominant as 
the satellite approaches the magnetic equator. Thus radi- 
ation intensities in excess of 100 roentgens per hr have bee 
observed by the satellites. 

In regard to cosmic radiation, it is interesting to note that, 
although the Russian satellites also carried instruments for 
measuring this effect, the preliminary announcements 0° 
scientific results made by the Russian scientists did no: 
indicate that they had discovered this high intensity region. 
Their preliminary information was confined to altitude. 
below 600 miles. At a recent meeting sponsored by the IG) 
committee held in Moscow, the Russians entered into a mor 
thorough description of the equipment onboard the Sputnik: 
They revealed that the Sputniks carried tape recorders whic! 
stored scientific information around the orbit and played i: 
back as the Sputnik passed over Russia. Also, the Sputnik: 
carried continuous transmitters which sent the informatio) 
to Earth steadily for any receiver which the satellite happene: 
to be passing over. The cosmic-ray measurements wer 
being transmitted on the continuous transmitter in the sam: 
way that information is telemetered from Explorers I and IV 
but the cosmic-ray information was evidently not put on th: 
tape recorders as was done with Explorer III. Thus th: 
Russians could observe cosmic-ray data only when the satellit: 
was near its lowest point. When the satellite was passing 
through the high altitude region, and presumably detecting 
the intense cosmic-ray field, the satellite was not over Russia. 
Instead, when the satellite was passing through this interest- 
ing region, it was over some part of the world near or below 
the equator. 

It is true that many radio receiving stations over these por- 
tions of the world were recording the signals sent by the 
Russian Sputniks. However, these signals were meaningless 
without the telemetering code, and this code was not pub- 
lished by the Russians. Furthermore, in spite of the request 
by this country and others, the Russians would not enter into 
any arrangements to exchange preliminary satellite data. 
Thus, although the Sputnik was observing this surprising 
high intensity cosmic-ray field, the Soviet government, by 
its own decision, refused to accept this data or make it pos- 
sible for anyone else to interpret it. As a result, although 
the Soviet scientists were the first to launch an artificial Earth 
satellite, they were the last to learn about this vital dis- 
covery of the high intensity cosmic-radiation field. Instead, 
it remained for the American scientists to discover and an- 
nounce this completely new geophysical phenomenon. 

The Jupiter-C research vehicle then, however unsophisti- 
cated it may seem as compared with Vanguard, and however 
diminutive it may be as compared with Sputnik, has been 
to date exceedingly successful in accomplishing its objective: 
The exploration of space and the acquisition of fundamental 
knowledge of our planet’s environment in the universe. 
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General Solution for Optimization of 
Staging of Multistaged Boost 
Vehicles 


CARL H. BUILDER’ 


The Marquardt Corporation 
Van Nuys, Calif. 


A general solution for the optimization of staging of multistaged boost vehicles based upon 


weight-dependent criteria is presented. 


Optimization may be in terms of any criterion that can 


be expressed as proportional to stage or step weight, such as cost, volume, empty weight, propellant 
weight, as well as the usual criterion of gross weight. The solution is applicable to a vehicle of any 
number of stages having different performances and characteristics. The possibility of stages with 
variable performance, such as an air-breathing engine stage, is included. Analytic procedures 
for the solution of multistage vehicles are outlined, and an example of the optimization of staging 


for minimum cost is presented. 


\N THE past few months, the staging of step rockets for 

minimum gross weight has received considerable atten- 
tion (1 to 4).2 However, these solutions are applicable to 
only one special case of the general problem of optimum 
staging: When minimum gross weight is the optimization 
criterion. 

The optimization of rocket staging for minimum gross 
weight has been treated by a number of authors during the 
past two years. Malina and Summerfield, in an early paper 
(5), optimized staging for maximum burnout velocity with 
homogeneous stages. Vertregt (6) was perhaps the first 
to outline the general solution for minimum gross weight 
staging; however, explicit solution was limited to stages hav- 
ing equal specific impulses. Goldsmith (7) offered a solution 
for two-stage rockets when the structural weights are pro- 
portional to the fuel weight and the engine weights propor- 
tional to the stage gross weight. Schurmann (8) considered 
burnout accelerations as a variable in the optimization for 
minimum gross weight when all stages have equal specific 
impulses. Weisbord (1), Ten Dyke (2), Subotowicz (3) and 
Hall and Zambelli (4) have all presented, in the last year, 
general solutions for minimum gross weight with nonhomo- 
geneous stages. Thus, it would appear that if gross weight 
were the criterion for optimum staging, solution would now be 
amenable through any one of nearly a half-dozen papers. 

However, minimum gross weight is not likely to be the only 
criterion of interest in the optimization of staging. Under 
certain circumstances other criteria, such as cost, volume, 
empty weight, etc., may be equally important. In the de- 
sign of boost vehicles for any ambitious space program in- 
volving large payloads in orbit or in extraterrestrial trajec- 
tories, the cost of the boost vehicle is likely to be one of the 
most important criteria for staging. While it is recognized 
that optimized staging does not necessarily result in the opti- 
mum system, the ability to optimize the staging for various 
criteria provides a set of sub-optimum systems which can aid in 
the final system selection. The ability to define optimum 
staging for various criteria is the objective of this paper. 

A general solution for the optimization of staging for any 
criterion proportional to stage weight is derived. The solu- 
tion is applicable to any number of dissimilar (performance- 
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wise) stages, and takes into account the possibility of a stage 
or stages having variable specific impulse, such as an air- 
breathing engine stage. Since a large number of pertinent 
parameters, such as cost, empty weight, propellant tank 
weight, engine weight, volume, etc., can be approximated 
over a reasonable range as being proportional to stage 
weight, this solution provides a means of optimizing staging 
for many criteria of interest, including the criterion of mini- 
mum gross weight. 


General Solution for Optimum Staging 


The solution is taken to be that set of stage velocity incre- 
ments AV; which corresponds to a stationary value of the 
objective function K subject to the constraints imposed by 
the stage parameters, J;, 8;, »; and the system final velocity 
V.. The objective coefficients 7; may be specific volume, 
specifie cost or any other parameter per unit weight of stage 
i. So long as these coefficients are independent of the ab- 
solute magnitude of the stage weights 


n 
K= >» Hiwi {1] 
i=1 


If the system is defined for a payload weight of unity, then 
from momentum considerations, for the ith stage 


n 
j= du j=it 
—— = loge | —— = [2] 
Tig 
AVi 1+ Byw; + Wj 
j=l j=itl 
AVj 
j=1 du = [3] 
Tq 
AVi 
j= 
If the function z; is defined as 
vi one 
then Equations [2, 3 and 4] may be combined to yield 
n 
w;, = — 1) II [5] 
j=itl 
497 
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Substituting Equation [5] in Equation [1] 


II 
i=t+l 


[6] 


[n(x — 1) x] 
i= 


The objective function K is a stationary value when, for i = 


oK/ oz; = 0 [7] 
subject to the constraint 
n 
Vi= > AV; [8] 
i=1 


However, since 2; is a function only of 7;(8; is a constant), 
Equation [7] is equivalent to 


oK/dy: = 0 [7a] 


The method of LaGrange’s multipliers provides a solution for 
constrained extrema in terms of n simultaneous equations of 
the form 


oK op 
oy: [9] 
plus the constraining equation 
{10} 
In this case 
n 
¢=)> AV; [10a] 
t=1 
If, by definition 
=1/(1 — (11] 
then Equation [6] may be differentiated to 
oK n i-1 t-1 
uf II +> Ee [12] 
k=t i=1 k=i+1 
where the value of the product 
i-1 
I 
k=i+1 
is taken to be unity forj + 1> 7-1. 
Differentiating Equation [10a] 
op oAV; 
U3} 


The evaluation of this total differential requires examination 
of Equation [3]. If the impulse of the stage J; is constant, as 
is frequently assumed for the rocket motor, then Equation 
[3] may be integrated to 


AV; _ 
Lig [3a] 
and the differential becomes 
op oAV; 
oY: [18a] 


However, if the impulse varies with velocity, as is character- 
istic of air-breathing engines, evaluation of the Equation [3] 
involves differentiation under the integral sign 


i i 

avs Av; 
dAV; Ji-1 oAV; oAV; lig 


i 
> avi 


i=1 


> avi 
j=1 


Since the impulse J; is a function of the velocity u and not 
the stage velocity increment AV;, the partial differentia] 


under the integral is zero. Also 
t 
AV 
(15) 
i-1 
and 2 >> AV;,/dAV; = 0 16] 
j=l 
dy; 1 
Thus = {— «= AVi 
us aAV, i [1 ta] 
For a stage with variable impulse, Equation [13] becomes 
do 
— = — = AVi {1 sb] 
oY: ov: 


Equation [13b] is equivalent to [13a] if the impulse J, is 
taken at the burnout velocity of the ith stage. Substituting 
Equations [12 and 18a] in Equation [9] results in n equations 
of the form 


plus the constraining Equation [8] 


i—1 


Alig = 0 


k=j+l 


ni + 1) 


j=1 


[17] 


The LaGrangian multiplier \ may be eliminated from these 
n + 1 equations by substitution, leaving n — 1 equations 
of the form 


(18} 


i-1 
m j=1Lm k=j+l1 


plus the nth equation, Equation [8]. Evaluation of Equation 
[18] for each of the n stages is relatively simple, since the 
right-hand side of the equation contains no terms higher 
thani — 1. Thus, if AV; is assumed, y2 may be calculated, 
and subsequently ys may be calculated, and so on until yn. 
The resulting values of AV; may be substituted in Equation 
[8] to calculate V, and, thus, check the validity of the as- 
sumed value for AV;. Successive approximations lead to 
that value for AV,, which, through Equation [18], leads to 
satisfaction of Equation [8]. 


Discussion 


The procedures for solution of a specific boost system by 
means of Equations [8 and 18] have been outlined. Solution 
by successive approximation techniques may be carried out 
by automatic computers or by hand. Simple iteration meth- 
ods do not generally converge and must be supplemented by 
damping the indicated change in AV, (in computer calcula- 
tions), or by judgment (in hand calculations). For hand 
calculations, tables or graphs of the functions x; and y; as 
functions of y;, with 8; as a parameter, are very useful. It is 


(7) 
[14] 


u= > avi 


i=1 


AV; 

j=l 
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Table 1 Assumed parameters for a three-stage boost 
vehicle 
Stage, 7 1 2 3 

specific impulse, J; 285 298 250 

residual weight, 8; 0.08 0.11 0.09 
objective coefficient, 7; 

(stage cost in $/Ib) 10 17 vf 


not feasible to present graphs of the desired precision to 
accompany this paper. Such graphs are easily prepared 
using Equations [4 and 11]. 

Note that the objective coefficients ; appear in Equation 
[1] of the solution in ratio form. Thus, the absolute magni- 
tule of the objective coefficient is never required for estab- 
lishing the optimum staging. This is useful, for example, 
when agreement on absolute costs is not attainable. Often, 
agreement may be confined to, “This costs 20 per cent more 
thin that.” Under such conditions, optimum staging can 
be established, even if the total system cost cannot. 

Ten Dyke (2) and Hall and Zambelli (4) have proposed 
siniple procedures for taking into account the real effects of 
drag and gravity when using the ideal rocket equations, such 
as Equation [2]. These procedures involve the use of an 
eff-ctive terminal velocity which exceeds the real terminal 
velocity by an amount that corresponds to the energy con- 
suined by drag, gravity and turning. The incorporation of 
such effective velocities in the optimization techniques pre- 
sented here is acceptable provided that the relationship be- 
tween the effective and real velocities is accounted for in the 
determination of the specific impulse of variable-impulse 
stages. This may require that the relation between effective 
and real velocities be established over the boost vehicle 
spectrum. This, of course, is tantamount to assignment 
of drag and gravity losses over the entire trajectory. 

As an example of the application of optimum staging 
for minimum vehicle cost, consider a three-stage orbital 
boost vehicle with the characteristics listed in Table 1. The 
effective velocity of the system is taken to be 32,000 fps which 
includes a reasonable allowance for drag and gravity effects 
along the trajectory. 

If the example vehicle is staged for minimum gross weight. 
by setting 

= 1 [19] 


the resulting mass ratio (gross weight to payload weight) is 
79.0. If the assigned cost figures are applied to this staging, 
the total cost of the booster vehicle is $850 per Ib of payload. 

When the vehicle is staged for minimum cost, by setting the 
4; values equal to the stage costs per pound in Table 1, the 
mass ratio is increased to 79.6, but the total cost is reduced 
to $842 per lb of payload. However, the differences in mass 
ratio and cost are sufficiently small, in this case, so that choice 
between the two sub-optimum systems is probably arbitrary. 
This is because the assigned costs for each stage are of the 
same order of magnitude. 


Now, however, let us assume that the relative stage costs 
change markedly by recovering and reusing the first or bottom 
stage of the booster system. If we assume that recovery 
incurs no penalty in stage performance and that fuel costs are 
negligible when compared to hardware costs, then 


m = (20) 
m 
where m is the number of flights made by the first stage. 
The differences in mass ratio and cost between vehicles staged 
for minimum gross weight and for minimum cost increase 
steadily with m, as shown in Table 2. 

Although the example illustrates optimum staging for 
minimum cost, optimization could as easily have been for 
minimum empty weight, minimum volume, minimum fuel 
weight or any other parameter which can be expressed per 
pound of stage weight. 


Nomenclature 

g = acceleration due to gravity 

I; = propellant specific impulse of the 7th stage 

K ~~ = objective function of the booster system per unit weight 
of payload 

n = total number of stages in the booster system 

u = velocity 

AV; = velocity increment due to the ith stage 

V, = final velocity attained by the payload 

w; = net weight of the 7th stage, including propellant, but ex- 
cluding payload and all stages above the 7th stage. 
Expressed as weight per unit weight of payload 

zi; = function of y; and 6; 

yi = function of y; and 6; 

8; = fractional weight of non-propellants of the 7th stage 
which is jettisoned at the burnout of the 7th stage 

yi = integral function of u and /; 

ni = objective coefficient of the ith stage. Expressed as 
volume, cost, etc., per unit weight of the 7th stage 

d = LaGrangian multiplier 

o = constraining function = 0 


The mathematical symbols d, 0, 2, I, e, ete., are used in their 
usual contexts. 
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Table 2 Comparison of weight-optimized and cost-optimized systems 
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Number of flights for first stage, m 1 2 3 5 10 20 
System optimized for gross weight: 
mass ratio, lb/lb of payload 79.0 79.0 79.0 79.0 79.0 79.0 
cost, $/lb of payload 850 525 416 330 265 232 
System optimized for cost: 
mass ratio, lb/lb of payload 79.6 82.7 85.9 91.6 104.0 124.0 
cost, $/lb of payload 842 494 369 261 172 118 
System comparison: 
difference in mass ratio, lb/lb of payload 0.6 3.7 6.9 12.6 25.0 45.0 
difference in cost, $/lb of payload 8 31 47 69 93 114 
499 
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The Useful Life of Solid Propellants 


at Very High Temperatures—Part 1 


M. VISNOV! 


Frankford Arsenal 
Philadelphia, Pa. 


The advent of the heating problem associated with supersonic speeds directs attention to the ther- 
mal stability of propellants under unusual conditions. It requires a sharp definition of the point 
in time-temperature coordinates where the propellant ceases to function correctly. The auto- 
ignition point is not sufficient to meet the designer’s needs, as there are two successive zones of 
danger in the decomposition of propellants. The first zone is where propellant function is faulty; 
the second is where the propellant becomes a hazard to adjacent equipment. Knowledge of these 
zones is vital to the design of propellant powered devices. Examples are given of the data which the 
problem demands. Typical data obtained to date are presented, as well as a brief description of 
automatic equipment employed in the current project for obtaining the data. 


HE HEAT that accompanies the operational conditions 

of current and planned supersonic aircraft and missiles 
has added stil] another problem area to the list of those facing 
technical personnel in these fields. Survival and reliable 
functioning under such extremes are now required of the com- 
pact source of work we know as the propellant powered or 
cartridge actuated device (CAD). For those that deal in 
structures, the emphasis has been on the search for materials 
which maintain form and strength over the increased tempera- 
ture range. The field of propellant powered devices faces an 
even more pressing problem. Here, relatively sensitive chem- 
ical compositions which are inherently unstable to various 
degrees must not only maintain form and strength at these 
high temperatures, but, in addition, when called upon must 
decompose into specified volumes of gas at controlled rates. 
When it is realized that there exist today high temperature 
requirements for propellant actuated devices severe enough to 
pose the hazard of spontaneous ignition, then some apprecia- 
tion of the problem facing the propellants field can be ob- 
tained. 

A study of the high temperature problem as applied to pro- 
pellants for cartridge actuated devices must deal with certain 
conditions peculiar to these systems which tax the stability of 
a propellant to an unusual degree. The severity of these con- 
ditions provides a virtual proving ground for propellant life 
under extremes of exposure. These factors include: 

Location of cartridge actuated devices. 

Extension of high temperature range of performance. 

Repetition of high temperature exposures (recycling). 

The sealed cartridge case. 

Limitations for protection of propellant powered devices. 

Some idea of how each of these factors contributes to the 
unique demands on propellants for CADs will be discussed. 


Location of Cartridge Actuated Devices 


Where possible, personnel ejection devices are placed in 
areas served by forced cooling and hence have some means of 
protection against heat. However, certain cartridge powered 
initiators and thrusters for canopy removal and related func- 
tions are placed by necessity close to the canopy itself. 
Other propellant systems perform multiple functions simul- 


Presented at the ASME-ARS Aviation Conference, Dallas, 
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taneously or in series at divergent points. It is unavoidable 
that parts of this network of devices must lie in uncooled con- 
partments. Here, the rise in temperature which the propel- 
lant must withstand depends upon its distance from the skin. 
For those CADs which must be mounted a few inches from 
the skin, the problem is, of course, aggravated. A related 
problem is the potential mounting of propellant devices to 
perform certain functions close to heat producing mechanisms, 
such as engines or afterburners. More severe temperature 
extremes are encountered by stores or payload ejection de- 
vices. Such devices are often mounted with a minimum of 
shielding close to the skin directly exposed to air flow. Speci- 
fications exist today for such devices to withstand tempera- 
tures higher than 300 F. There are requirements for certain 
CADs now under development for exposure to 500 F for 
periods in excess of 4 hr. 

Requirements for propellant powered devices for missiles 
vary with use. Gas generators for fuel pressurization are usu- 
ally positioned within the interior of the body and require 
little protection. Certain designs, however, must be shielded 
from internal heat effects. Nose cone and warhead detaching 
devices, if skin mounted, require protection from aerodynamic 
heating. A current specification for a cartridge powered de- 
vice that manipulates missile airfoils bears a high temperature 
requirement of 350 F for 1 hr. A potential application for 
propellant devices is in the functioning of components to alter 
flight during atmospheric re-entry. These devices could be 
actuated by temperature sensing mechanisms. For such a 
system, it follows that a greater resistance to high tempera- 
ture by the propellant would allow a wider latitude in inere- 
ments of time and, consequently, distance, within which the 
necessary flight adjustment could be made. 


Extension of the High Temperature Operation 
Range 


For most rocket or gun applications, ordnance require- 
ments stipulate satisfactory performance over a temperature 
span of —65 to 160 F. Inaddition, it is desired that tempera- 
ture coefficients of performance be ata minimum. However, it 
is not uncommon for performance at either extreme to be mar- 
ginal. Propellant devices exposed to much higher tempera- 
tures than described may also be required to function at the 
peak of heat exposure. Propellants for this application face a 
temperature range of performance hitherto unknown to the 
field. 
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Recycling 


Propellant powered devices which are an integral part of the 
operation or mission of supersonic vehicles are subjected to 
only one thermal cycle before being fired. But other CADs, 
such as personnel ejection or certain stores ejection devices 
serve an emergency purpose. As such, they are required to 
undergo an unknown number of missions and attendant high 
temperature cycles. If this is coupled with potential ground 
soik exposures between missions, then the heat problem as 
applied to CAD propellants assumes its real, albeit disturbing, 
proportions. Although this possibility of continuous cycling 
and its additive destructive effects on propellants could be 
kept under control by a frequent replacement program, a 
hizh temperature resistant propellant could keep such a pro- 
grim within reasonable bounds or else do away with it alto- 
gether. In reference to ground soak, studies show that in 
parked aircraft, temperatures well over 200 F have been ob- 
served beneath the closed canopy at pilot heat level. In one 
day, under these conditions, data show temperatures exceed- 
ing 160 F for 6 hr and exceeding 200 F for 3 hr. The number 
of times that this cycle may be repeated cannot be determined, 
since it depends on the operational status of the aircraft and 
ground handling procedures, among other factors. The fact 
remains that, whether the aircraft is flown or not, the CAD 
propellant, unlike weapons propellants, is a fixed aircraft 
component subject to whatever extremes occur. 


Sealed Cartridge Case 


A specification requirement for CADs is that the case con- 
taining the propellant be sealed to withstand 14 psi pressure 
differential. This is to forestall moisture condensation within 
the cartridge case due to extremes of humidity and tempera- 
ture, between ground level and flight altitude. Sealing, how- 
ever, works to disadvantage at elevated temperatures, since 
this accelerates the decomposition rate of many propellants. 


Limitations for Protection of Propellant 
Powered Devices 


The possibility of protection against heat by insulation has 
been considered and in certain instances is feasible. How- 
ever, it is difficult to extend insulation over a broad range of 
application, since CADs are designed under sharp limits of 
allotted space, usually within fractions of an inch. Dimen- 
sional increases due to the addition of thermal barriers must 
be brought into line with with the ever-present drive toward 
miniaturization. 

The Research and Development Group of Frankford 
Arsenal is charged with the mission of research and develop- 
ment in the field of cartridge actuated devices. The need for 
propellant-initiator systems capable of operation at tempera- 
tures far above the presently accepted maximum of 160 F has 
been recognized for some time. As a result, a program is cur- 
rently under way for the Air Force whose object is the develop- 
ment of propellant-initiator systems for use in environments 
up to 400 F. This program is divided into five phases which 
can be summed up as follows: 

1 Determine the useful life at high temperatures and 
auto-ignition hazard time-temperature relationship of cur- 
rently standard propellants for CADs. 

2 Determine the same limits for other existing propellant 
types. 

3 Establish the feasibility of substituting more thermally 
stable existing propellants for the present ones in CADs. 

4 Develop new and radical propellant-initiator systems 
for high temperature resistant CADs. 


* Fendick, R. B. and Stuck, D. H., ‘Phase V High Temper- 
ature Climatic Test of an F101A Aircraft,’’ Directorate of Flight 
and All-Weather Testing, Wright Air Development Center, Tech- 
nical note no. WADC-TN-56-404, Nov. 1956. 
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5 Establish the safe operating time-temperature limits of 
standard and new systems in propellant powered devices. 

This paper describes some of the data required, the ap- 
proach taken to obtain this data, and some of the equipment 
developed for that purpose. 


Discussion 


Nitrate Ester Propellants 


The bulk of high temperature work on propellants in the 
past has usually taken either of two directions—research on 
storage stability and surveillance procedures or direct work 
on the phenomena of auto-ignition or cookoff. 

For storage stability, the maximum temperature of work is 
in the area of the 160 F requirement imposed by the armed 
services. In this work as in stability tests conducted above 
this temperature, heat is employed as a means of further ac- 
celerating those propellant decomposition processes which 
would normally take long periods of time to bring about. 
Norms are set in terms of these accelerated results, and esti- 
mates of propellant life at ambient temperatures are made by 
correlating the results obtained under these extreme test con- 
ditions. For work at temperatures higher than the 160 F 
range, investigators over the years have resorted to rigorously 
defined empirical test procedures to obtain reproducible re- 
sults at these accelerated conditions. These include the well- 
known methy!] violet heat stability test at 274 F (134.5 C) 
for nitrocellulose single-base propellants, and at 248 F 
(120 C) for double-base propellants, and the gas evolution 
tests, such as the Taliani test at 230 F (110 C) or the vacuum 
stability test at 212 F (100 C) for single-base and 194 F (90 C) 
for double-base propellants. 

The literature shows lesser emphasis on high temperature 
work dealing with auto-ignition of propellants. The hazard of 
spontaneous ignition during propellant manufacture has pro- 
vided the impetus for most of this work. Some research of a 
more basic nature deals not with auto-ignition as an end in it- 
self, but, under controlled conditions, as a means of deter- 
mining physico-chemical constants of propellant compositions. 
Typical is the establishment of the apparent activation energy 
for ignition by applying the Arrhenius equation for the in- 
fluence of temperature on rate of reaction, in this case reduced 
to 


1 
log = + constant 


where 
t = time to ignition 
Eq = activation energy 
R = gas constant 
T = absolute temperature 


A plot of log t vs. the reciprocal of the absolute temperature 
yields a straight line from whose slope the activation energy 
E, can be calculated. 


Composite Propellants 


The variety of fuels developed for composite propellants 
since World War II pose an extremely difficult problem to 
those concerned with correlating propellant stability under 
various conditions. The complicating factor is that the ac- 
cepted chemical stability tests of today, either the color indica- 
tor or gas evolution types, are applicable only to the nitrate 
ester family of propellants. Unfortunately, the chemical basis 
for applying these tests and norms to the composite pro- 
pellants does not appear to exist, since in many cases the latter 
do not evolve gaseous nitrogenous decomposition products. 
In the absence of a standard yardstick, other means of com- 
parison must be employed. For the program described in 
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this paper, progressive degradation of different propellant 
types under severe heating is being measured by infrared 
spectrophotometry and by mass spectrometry techniques. 


Factors Affecting CAD Propellant Life 
at High Temperature 


For our purpose the title “propellant life at high tempera- 
tures” may be subdivided. We can describe the profile of 
high temperature life as in Figs. la and Ib. If we raise the 
temperature surrounding the device above 200 F, there is 
first an interval where, outwardly, nothing takes place. But 
it is during this period that present CAD propellants begin to 
change. So long as these changes are not significant enough 
to affect ballistics when fired, the propellant can be said to 
still retain ‘‘useful life.’’ Useful life is terminated when this 
decomposition has increased to the point where ballistics fall 
outside of minimum or maximum tolerances. Beyond this 
point, as exposure to heat continues, the rate of decomposition 
increases sharply. Though seemingly intact, the propellant 
device is now useless. Ultimately, of course, if the time- 
temperature conditions are severe enough, the zone of auto- 
ignition hazard is entered. 

The earlier heating stage where useful life may be measured 
in minutes or hours, depending on heat input history, is an 
unexplored area. This useful life can fall well within one or 
more flight times of supersonic vehicles. The changes in the 
propellant in these exposures short of auto-ignition include 
chemical degradation, migration of ingredients, or plastic de- 
formation of calculated grain geometries. Any one can render 
the charge incapable of functioning as originally designed. 

Various factors are being investigated for their effects on 
the time-temperature coordinates of the propellant life zones 
pictured above. These factors are: 

Rate of temperature rise. 

Rate of rise plus exposure time at a series of temperature 

levels. 

Thermal lag through the case wall. 

Composition of the propellant. 

Propellant dimensions and mass 

What is unique is that their influence is being exerted in this 
system under constant-volume conditions within the limits of 
the CAD seal. For nitrate ester propellants, at least, it is 
known that if the products of decomposition remain in contact 
with the propellant under pressure, the decomposition rate 
will be accelerated. 


400° F. — 
MAX. - 7 Rote Of Rise To 
‘onst. Temp. 
MIN. 
PERFORMANCE 
SPAN 
(units optionel) 
MAX. 
ACCEPTABLE PERFORMANCE 
— 
USEFUL LIFE TERIORA TIONS 


Time 


Fig. la Propellant life at high constant-temperature levels, one 
system 


Experimental Methods 


Procedure 


The hazard zone of potential auto-ignition is determined 
under actual sealed CAD cartridge case conditions. The 
auto-ignition time-temperature relationship is evaluated 
from heat runs at various rates of temperature rise and ra‘e 
of rise coupled with sustained constant temperatures. The 
latter are carried out at regular intervals beginning at 200 F 
and extending to auto-ignition maximum. The next step s 
the determination of the zone of useful life. This is obtained 
from simultaneous exposures of numbers of loaded CAD ca- 
tridges. The specific time-temperature schedules are in tle 
safe zone determined from the foregoing auto-ignition dete :- 
minations, but are arranged to become increasingly sever, 
ending just short of the indicated auto-ignition conditions. 
From each exposure, certain of the cartridges are then tested 
ballistically for change in performance, while others are exan:- 
ined chemically and microscopically for extent of decompos'- 
tion and physical deformation of propellant. Concurrentl:., 
study is made of the heat transfer characteristics of the systein 
and effect on primer and igniter. 


Materials and Equipment 


Propellants 


The propellant types currently under study include: 

Nitrate ester types, single-base compositions (nitrocellu- 
lose) and double-base compositions (including both high and 
low nitrogen nitrocellulose and a range of nitroglycerin con- 
tents and additives). 

Composite propellant types (including perchlorate and ni- 
trate oxidizers and rubber, plastisol, resinous, nitropolymer, 
and newer, more exotic fuels). 


Thermal Programmer 


Valid and reproducible data for this program require the 
closely controlled application of heat. To accomplish this, 
an electronically controlled thermal programmer was designed. 
The unit is capable of providing linear rates of rise from 0.5 to 
30 F per min. It can be preset to level off at any temperature 
up to 750 F and maintain that temperature for any time inter- 
val from 1minto12hr. Time-temperature history and auto- 
ignition point are recorded on a strip chart. Accessory de- 
vices include a delayed start timer for lengthy overnight runs 
and a signal actuated at shutoff to operate safety and cooling 
devices. 


Auto-Ignition Apparatus 


This device is operated by the thermal programmer. In it, 


= CAD System 
= Propeliont 


~ °F fmin Rote of Rise 
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Fig. 1b Propellant life at constant-temperature levels, one 
system 


ARS JouRNAL 


3 
| sence 


nined 

The 
uated 
rate 

The 
200 F 
fep 
ained 
) 
n the 
lete- 
ver, 
tions. 
ested 
Xan i- 
i- 
ntl, 
‘stein 


cellu- 
and 
con- 


d ni- 
rmer, 


» the 
this, 
med. 
).5 to 
ture 
nter- 
wuto- 
y de- 
runs 
oling 


it, 


NAL 


Fig. 2 Gas release chamber 


a loaded cartridge can be safely heated until auto-ignition 
occurs. The heated vessel is a solid block of commercially 
pure copper containing a deep 1-in. diameter well. Within the 
well, the cartridge is supported by a retainer which permits 
the case wall surrounding the propellant to be suspended in 
air. The control thermocouple from the programmer is in 
contact with the cartridge wall. Heating is accomplished by 
a dual system, one system comprised of helically wound ele- 
ments in a heating jacket, the other, a set of six cartridge heat- 
ers buried in the copper block radially surrounding the well. 
The standard M38 cartridge case serves as a test container. 
For this purpose, it is modified to permit the gases from the 
ignited propellant to rupture the sealing disk. At auto-igni- 
tion, the blast of gases actuates a trigger, which, through a 
microswitch, instantly shuts off the entire heating, controlling 
and recording mechanisms. The operation is automatic and is 
carried out in an armored cubicle. 


Apparatus for Estimation of Propellant Useful Life 


For the determination of the zone of useful life of propel- 
lants, fully loaded cartridges are used. Safety considerations 
dictate that personnel cannot be present during these tests. 
For this purpose, a specially designed oven, located in an 
armored room, is operated by remote control from the pro- 
grammer. The severe heating schedules are relatively short, 
and must be terminated abruptly. When a particular sched- 
ule ends, the programmer shuts off the oven, signals a device 
which opens the oven door, and a blower purges the oven of 
heated air. Since the exposed cartridges are sealed, they still 
retain gaseous decomposition products. To transfer these 
products to the mass spectrometer for analysis, the cartridge 
is opened in a sealed-off chamber designed for this purpose 
(lig. 2). The sealed chamber is first flushed with helium. 
A sharp pin driven by a solenoid positioned above the cartridge 
riiptures the cartridge sealing disk, and the gases are released. 
‘lhe chamber is then opened to an evacuated line leading to 
tle mass spectrometer, and the gases are drawn into the in- 
sirument for analysis. 
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Results 


The first phase of the program is approaching completion 
with the establishment of time vs. temperature auto-ignition 
hazard data for CAD cartridges loaded with the most widely 
used solid propellant types. Over eight hundred runs have 
been made covering various rates of temperature rise and 
elevated constant temperature “soak.” Due to the precise 
programming of heat input by the apparatus described above, 
reproducible data have been obtained. Typical individual 
run data are shown in Table 1 for a CAD cartridge loaded with 
various propellants. 

The effect of rate of temperature rise on auto-ignition tem- 
perature for nitrate ester propellants is shown in Fig. 3. In 
general, a lower rate will result in a lower auto-ignition tem- 


@---- Single bose propellant, #10 


Double bose propeliont, M2 


AUTOIGNITION TEMPERATURE, ° F 


250—— 


10 


4 


LINEAR RATE OF RISE, ° E/MINUTE 


Fig. 3 Auto-ignition temperature vs. linear rate of rise, sealed 
CAD cartridge 


Table 1 Auto-ignition temperature individual run data, 
sealed CAD cartridge 


Auto-ignition temperature, F 
—Linear rate of rise, F/min— 

Propellant type 1 5 6.10 20 30 
double base 284 325 354 414 441 
284 325 «4350 411 442 

283 324 350 412 440 

284 325 355 410 441 

282 325 350 412 445 


ammonium perchlorate/. 460 514 554 602 637 
polymer 460 516 552 598 637 
460 515 553 601 647 

457 513 555 612 637 

463 516 559 596 652 


ammonium nitrate / 382 402 437 497 527 
polymer 382 401 437 493 542 
387 400 437 496 532 

379 4388 497 532 

387 399 488 485 527 


black powder 529. 554 572 594 616 
522 553 575 598 617 
522 550 573 598 615 
522 554 573 594 609 
520 550 571 590 611 


350 
300 
© = 
| 
a 
#3 
one 
ae 


a. Perchlorate/rubber: 2 hours at 350° F 


b. Perchlorate/resin: 12 hours at 400° F 


c. New composition: 20 hours at 400° F 


Fig. 4 Composite propellants before and after high tempera- 

ture exposure: a. perchorate/rubber—2 hr at 350 F. b. per- 

chlorate/resin—12 hr at 400 F. c. new composition—20 hr at 
400 F 


Table 2 Auto-ignition temperature, sealed CAD cartridge 
case 
Auto-ignition 
temperature, F 
Grain Linear rate of rise® 
Propellant dimensions, in. F/min 
type OD Length 5 10 20 30 
M2, double 0.040 0.046 321 347 403 436 
base 0.073 0.243 322 351 399 430 
0.203 0.4382 329 349 403 449 
0.245 0.478 326 350 396 447 
0.495 0.600 319 352 403 442 
M10, single 0.049 0.135 336 364 410 446 
base 0.202 0.475 337 369 418 455 
* Values represent averages of five runs. 


Table 3 Auto-ignition time, min, at constant temperature, 
sealed CAD cartridge 


Constant temperature level, F7———\ 
Propellant 225 250 275 300 325 350 


double base 720+ 286 30 1l 6 4 
singlebuse 720+ 720+ 480° 19 8 4 


2 Values represent averages of five runs. 

> Average of seven cartridges. Eight cartridges sur- 
vived 720 min. Fifteen tests made due to this anomalous 
behavior. 
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perature. At these low rates of rise, the temperature gradient 
induced by the cartridge case is minimized, and the observed 
auto-ignition temperature is more a characteristic of the pro- 
pellant itself. At high rates of rise, the indicated auto-igni- 
tion temperature reflects the heat transfer rate of the cartridge 
system. (In this connection, it should be borne in mind that 
the work described herein is necessarily the study of a car- 
tridge exposed to a heated environment. As a result, the 
values obtained are not ascribed to propellant alone, but rather 
to the propellant-cartridge system.) Asshownin Table 2, tlie 
grain size of the nitrate ester propellant described here had no 
significant effect upon auto-ignition temperatures. It shou!d 
he noted, however, that the grain sizes are relatively small. 

Other data obtained show that nitrate ester propellan's 
now standard for CAD will not meet 400 F exposure. Sever:: 
series of propellant loaded sealed cartridges were heated at 01 ¢ 
rate of rise at 30 F per min up to level-off at predetermined 
constant temperatures. Isothermal exposures were term'- 
nated at 12 hr if auto-ignition did not occur sooner. The 
sharp drop in total life at constant temperature exposurcs 
above the 250 F level is shown in Table 3. Note the sligl\t 
advantage of single-base over double-base propellants below 
the 300 F level. 

When tested under similar conditions, the composite pro- 
pellants as a class show greater tolerance for elevated temper::- 
tures prior to auto-ignition. Many can be held at 300 F or 
somewhat higher for considerable periods. (Exact data and 
compositions are, of course, classified.) Some of these, how- 
ever, do not survive these exposures in usable condition. U)- 
on examination after exposures up to 400 F, they have shown 
slumping, melting, exudation of ingredients or such gener 
degradation as to be unrecognizable as propellant grains. 
Some are pictured in Fig. 4 before and after exposure. New 
composite propellants specifically tailored for high tempera- 
ture resistance are now being developed. Some are being 
evaluated in this program and to date are a marked improve- 
ment over earlier compositions. One is pictured in Fig. 4 
for comparison. 

An interesting sidelight is the resistance to auto-ignition 
shown by black powder, various grades of which are used as 
igniter charges in CAD. A sample of the data is shown in 
Table 1. Ona temperature rate of rise basis, there was little 
difference between the various granulations. 


Summary 


The first phase of this program can be summed up asa search 
for propellant systems capable of simple survival at high tem- 
peratures. By-products of this evaluation are the data 
necessary for establishment of auto-ignition hazard zones. 
With the completion of this phase, some observations can be 
made. Present nitrate esters fall short of the new high tem- 
perature requirements for propellant actuated devices. Cer- 
tain composites show promise up to 400 F exposures for mod- 
erate durations. Even here the composite propellants face 
certain problems at high temperatures. Oxidizers must sur- 
vive these temperatures without phase change or other 
phenomena. Fuels or binders will require changes in their 
formulations. Propellant curing cycles may need adjustment. 
(For this program, one of the earlier composites which showed 
exudation of one of the fuel ingredients was reformulated and 
now shows promise of increased high temperature life.) 
For CAD application, composite propellant grains in very 
small sizes and burning web thicknesses of close tolerance will 
have to maintain these dimensions without swelling or slump 
over the increased temperature range. There is some diffi- 
culty in obtaining composite propellants in small geometries, 
but some new extrudable high temperature resistant formula- 
tions have been obtained. Similar difficulty exists in obtain- 
ing burning rate data of composite propellants at the high 
pressures characteristic of some cartridge actuated devices. 
A very difficult problem is an acceptable temperature coeffi- 
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cient of performance over a temperature spread of approxi- 
mately 460 F. 

The second phase of this work is the delineation of zones of 
substandard ballistic functioning at high temperature expo- 
sure levels short of auto-ignition. Concurrently, these heated 
propellants are being studied by instrumental techniques to 
determine the mechanism of ballistic failure. This requires 
matching the ballistic performance of nitrate esters with the 
new propellants (phase 3 of the program). Some of this has 
already been accomplished, and firings of CADs exposed to 
high heat schedules are now in progress. It is hoped that a 


report on the completion of this phase will be made in the near 
future. 

In this problem of propellant actuated devices exposed to 
severe heat, a great deal of data must be obtained and exa- 
mined. Elevated temperature conditions are now real, and 
these data must come from a considerable amount of work in 
areas previously subjects for theoretical study. Stability and 
performance problems in these areas must be overcome if we 
are to make use of the merits of propellant power packages at 
high temperatures. The approach described here encompas- 
ses such a study with a view toward end item development. 


Vaporization Rate Limited Combustion 


in Bipropellant Rocket Chambers 


E. MAYER! 


Rocketdyne Division, 
North American Aviation, Inc. 
Canoga Park, Calif. 


A one-dimensional model of the bipropellant rocket combustion chamber is analyzed on the as- 


sumption that vaporization of liquid droplets is the rate controlling step in the combustion process. 
The dependence of the gas flow variables on position coordinate in the chamber is determined in 
elementary analytic form subject to the assumptions that the vaporized propellants are instan- 
taneously mixed and reacted in stoichiometric ratio, and that droplet motion is governed by full 
entrainment in the product gas flow; i.e., droplet velocity is equal to local gas velocity. The formal 
solutions obtained permit comparison and contrast with the more fundamental case of monopro- 
pellant combustion. It is concluded that in the bipropellant system gas velocity distribution is de- 
termined primarily by the droplet spectrum and volatility of the more slowly vaporizing species, 
and that the temperature distribution is strongly dependent on the droplet size spectrum and 
volatility of the more rapidly vaporizing species. Depression of the temperature levels associated 
with the presence of excess unreacted propellant gas near the injector implies the existence of flame- 


holding stability limits due to excessive mismatch in the local vaporization rates. The possible 
scope of applications and limitations of the vaporization rate controlled combustion model is dis- 
cussed on the basis of relevant data available in the literature on combustion in rockets. 


URING the past decade a great deal of experimental and 

theoretical work has been accomplished in basic studies 
of droplet and spray combustion (1)? with the objective of 
elucidating a variety of combustion problems in jet propul- 
sion. However, applications of these basic studies to com- 
bustion in liquid propellant rocket chambers have been lag- 
ging because of complexities arising from the interactions 
among the atomization, fluid dynamic and chemical effects. 
In recent investigations Priem (2,3) obtained successful cor- 
relations of combustion efficiency with chamber length based 
on a theoretical combustion chamber model in which vapor- 
ization of hydrocarbon droplets was the rate controlling 
process. It is the purpose of this paper to present a general 
description of the vaporizaton rate limited bipropellant 
combustion chamber by analysis of a representative model 


Presented at the ARS 13th Annual Meeting, New York, 
N. Y., Nov. 17-21, 1958. 

1 Engineering Specialist, Advance Design. Member ARS. 

? Numbers in parentheses indicate References at end of paper. 
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which retains the essential features of the bipropellant process 
without undue mathematical complications. Elaborate treat- 
ment does not appear warranted at this time, owing to un- 
certainties in droplet size and velocity distributions obtained 
upon injection of liquid propellant streams into the prevalent 
gas environment in the rocket chamber. Nevertheless, by 
the use of simplified representations of droplet distributions 
and vaporization laws, it is possible to ascertain the effects 
of bipropellant volatilities on the gas velocity, efficiency and 
temperature steady-state distributions in rocket combustion 
governed by vaporization rates. The simplified treatment 
herein presented permits the comparison and contrast of the 
bipropellant process with the more elementary case of 
vaporization rate controlled monopropellant combustion. 
The applicability of the results obtained is limited to high 
temperature reaction under relatively rapid mixing condi- 
tions. Consideration of published data on the internal 
chemical and flow conditions in experimental chambers 
(4,5) reveals the importance of mixing and chemical time 
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lags which, to some extent, mask the role of vaporization rates 
in the dependence of combustion efficiency on combustor 
length. Nevertheless, as shown by Priem (2,3), useful tech- 
nical correlations may be deduced by treatment of a vapori- 
zation rate limited combustor model. More generally, it is 
reasonable to expect that clarification of bipropellant com- 
bustion governed by vaporization rates will simplify the inter- 
pretation of combustion phenomena arising from mixing and 
chemical time lags present in the rocket combustion chamber. 


Model of the Vaporization Rate 
Limited Combustion Chamber 


In the combustion process limited by vaporization rates, 
the continuous, steady-state conversion of liquid propellants 
to reacted gases is visualized to occur in the following se- 
quence. Upon injection of the liquid jets at the upstream 
end of the chamber, a droplet spray is formed here with initial 
droplet size and velocity distributions dependent on both jet 
impingement configurations and prevalent “steady-state” 
gasdynamic conditions. In the subsequent history of the 
droplets, evaporation is governed by heat transfer from the 
surrounding gas stream, and velocity modifications ensue by 
entrainment in the product gas flow. The vaporized pro- 
pellants are rapidly mixed and promptly reacted (in stoichio- 
metric ratio) wherever they are formed, thus providing for 
a continuous increase of the gaseous mass flow rate with in- 
creasing distance along the flow direction. The combustion 
is completed downstream of the injector where all liquid 
droplets capable of providing reactants have been vapor- 
ized; i.e., in a bipropellant, with nonstoichiometric overall 
oxidant-to-fuel ratio, combustion is complete when the de- 
ficient propellant has been vaporized. 

An idealized model of the combustion chamber representing 
the features described above is analyzed in this paper on the 
basis of the following assumptions: 

(a) One-dimensional, constant-pressure gas flow. 

(b) Ideal gas law. 

(ec) Instantaneous mixing and (stoichiometric) chemical 
reaction of the vaporized propellants within the local product 
gas flow. 

(d) Negligible axial heat conduction and radiation effects. 

(e) Vaporization rate independent of instantaneous loca- 
tion of droplet, i.e., dependent primarily on the time elapsed 
after injection. 

(f) Droplet motion governed by full entrainment in the 
main gas flow, i.e., droplet velocity equal to the local gas 
velocity. 

The idealizations listed in these assumptions appear to be 
more applicable to large-scale, high performance rocket en- 
gines than to small-scale model combustion chambers, where 
the nonrepetitive pattern of injected jets and the influence 
of relatively nearby chamber walls impair the validity of the 
one-dimensional treatment. To be sure, even in the large- 
scale engine with planar injector face, the presence of discrete 
liquid jets creates complex multidimensional flow conditions 
near the injector face. In particular, the continuous effective 
ignition of the vaporized propellants in the large-scale engine 
must be accomplished by anchoring the reaction process in the 
low velocity recirculation zones created in the neighborhood 
of the solid portions of the injector face. Although details 
of the recirculation process are not accounted for in assump- 
tions (a), (b), (c) and (d), the one-dimensional treatment 
does permit an appraisal of the relative importance of bipro- 
pellant vaporization rates on flameholding against the in- 
jector face, which is a necessary condition for combustion 
stability in the large-scale, high performance engine (see dis- 
cussion of flameholding stability). 

In regard to assumption (e), basic investigations of both 
combustion (1,6) and vaporization of droplets (1,7) indi- 
cate that the diameter vs. time relation of typical propellant 


droplets of initial diameter D; can be approximately repre- 
sented by the parabolic form 


Dt = D2 — Kt 


where K is a rate parameter dependent on physico-chemic:] 
properties of the droplet and of its ambient atmosphere. In 
application to rocket combustion, use of K data pertaining 
to vaporization of droplets in turbulent surroundings is more 
appropriate than use of rate parameter data determine 
from droplet combustion in quiescent atmospheres (see next 
section). 

Finally, in regard to assumption (f), it is clear that rapi| 
entrainment of the droplets in the main gas can be expecte:| 
only for the small droplet component of the spray forme | 
upon injection. In high performance engines, this spra, 
component appears to be of particular importance in achiev- 
ing the high gas generation rates and heat release rates pres- 
ent in the zone adjacent to the injector. Full entrainmen' 
in one-dimensional gas flow is an essential assumption fo 
avoiding mathematical complications which arise if dra; 
effects due to relative velocities between droplets and ga 
flow are taken into account. A more realistic treatment o 
the entrainment problem appears to be further complicate: 
by the unknown initial droplet velocity distribution withii 
the spray formed upon impingement. 


Representative Vaporization Rate Data 


In rocket engineering practice, a basic datum in the design 
of a combustion chamber adequate for efficient combustion 
of given propellants is the characteristic length defined by 


L* = V./A, 
where 


V, = chamber volume (including the converging section of the 
nozzle) 
A, = nozzle throat area 


Empirical values of L* are listed in Table 1 reproduced from 
(8). 


Table 1 Empirical values of L* at P, = 300 psia 
L*, in. 
Propellants (Approximate values) 
LOX-aleohol 50-100 
nitric acid-aniline 40-175 
nitric acid-JP; 50-60 
nitromethane >190 
LOX-JP, 35-55 


The extended ranges of L* values indicated in the table 
arise from such complicating hydraulic factors as initial axial 
velocities, droplet size distribution and mixing pattern of the 
injected propellant present under the experimental conditions 
yielding the listed data. Further complications arise from 
chemical lags incurred in the conversion of gas phase pro- 
pellants to product gases. In spite of these complications, 
the influence of vaporization rates on the reported L* values 
becomes apparent in view of representative vaporization rate 
data as given, for example, by Penner (7), shown in Fig. 1. 
The solid curves in Fig. 1 describe the dependence of droplet 
radius on time calculated in (7) on the assumption that heat 
transfer from the product gases to the vaporizing droplet 
occurs across a film thickness approximately equal to the 
instantaneous droplet radius. The dotted curve in the figure 
is obtained by application of Spalding’s transfer-number 
method (9), whereby droplet lifetimes under comparable 
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Table 2 Lifetimes of fuel droplets of initial diameter, D; = 200 u 
Lifetime, K = D;?/ty, 


Fuel Process ty, sec 
n-heptane Quiescent combustion in O, at 0.008 
ethanol Py = 300 psi 0.009 
| 
n-octane Vaporization in product gas at » 0.003 

ethanol T, = 3006 K, Py = 300 psia 0.004 


| 
| 
3 


em?/sec Source 
0.050  Goldsmith’s (10) experimental data corrected for 
0.045 pressure level by 0.25 pressure exponent (11) 


0.133 Ref. (7), see also Fig. 1 
0.100 


product gas conditions are interrelated by the transfer 
numbers’ B dependent on the vaporizing species. 

Comparison of the empirical L* data with the calculated 
droplet lifetimes indicates that both quantities increase in the 
propellant sequence, 1 hydrocarbon, 2 alcohol, 3 aniline, sug- 
g sting the importance of vaporization as a rate controlling 
step in the overall combustion process for these high per- 
formance propellants. 

The dominant role of the vaporization process is also re- 
vealed in the relatively long L* of monopropellants as nitro- 
methane and hydrazine (L* > 200 in.). However, low per- 
formance monopropellants of this type are known to undergo 
slow gas phase reactions which preclude the applicability 
0: single rate controlling step (see discussion). 

It may be noted here that the droplet vaporization rates cal- 
culated in (7) are considerably greater than those obtained 
when vaporization occurs under quiescent conditions of drop- 
let burning governed by molecular diffusion of vapor from the 
droplet surface to a surrounding reaction front, established 
across an intervening vapor film of some 3 to 10 radii in thick- 
ness (1,6). Although an abundance of quiescent droplet 
burning data has become available during the past decade 
(1,6), the applicability of such data to the interpretation of 
rocket combustion is questionable because of turbulent con- 
ditions in the rocket chamber. These conditions preclude 
the existence of an extended molecular diffusion field surround- 
ing the droplet. It is likely that the high temperature gaseous 
products penetrate within a radius or less of the droplet sur- 
face, thus increasing heat and mass transfer rates severalfold 
as compared with droplet burning. 

For purposes of orientation, lifetime data on hydrocarbon 
and alcohol droplets of initial diameter D; = 200 yw, vapor- 
izing under quiescent and turbulent conditions, are com- 
pared in Table 2. The rate constants K listed in this table 
are related to the lifetime t; by the parabolic law Equation 
[1], ic, by K = D,?/t; which, on theoretical grounds, is ap- 
plicable to droplet combustion. For the vaporization rates 
calculated in (7), use of the parabolic law is only approxi- 
mately valid as indicated in the nondimensional diameter 
vs. time curves in Fig. 2. The agreement shown in this curve 
is sufficiently close to warrant the application of the parabolic 
form for vaporization rates in the illustrative calculations of 
the next two sections. 

An order-of-magnitude estimate of injected droplet diam- 
eters based on D; = +~/Kt; may be made from droplet life- 
time data inferred from photographic studies of combustion 
in transparent experimental rocket chambers. The growth 
of gas velocity in the axial direction, associated with the 
progress of combustion, is conveniently determined by fast 
photography of moving luminous pockets, as described by 
Berman and Logan (12). Gas velocity data obtained under 
severe atomization conditions in liquid oxygen/alcohol com- 


3 The transfer numbers Boctane = 5.02, Bethanot = 2.95 ((9) p. 
130) inserted into 


lethanol = log (1 Bostane) 


loctane log (1 + Bethano1) 


aterrelate the times required for a given decrease in the droplet 
radius in curves (c) and (d) of Fig. 1. 
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bustion at P, = 300 psi (12) indicate that combustion is 
accomplished within 3 in. of the injector with an approximate 
droplet lifetime of t; = 0.002 in transit through the combus- 
tion zone. With the vaporization parameter K = 0.1 cm? per 
sec (see Table 2), the estimated maximum droplet diameter 
is D;~ 1.4 X 107? cm. 


Combustion of Uniform Droplets 


Many of the essential features of vaporization rate limited 
combustion are readily deduced by treatment of the elemen- 
tary case of uniform initial droplet size D; for each propellant 
type. The influence of droplet size distribution is conven- 
iently discussed in terms of the elementary results derived 
in this section. 


Monopropellants 


Since the heat released per unit mass of the vaporized (and 
reacted) monopropellant is constant, the gas temperature 7’, 
and density p, are also constant under the assumptions (a) 
and (b). If N liquid droplets of diameter D; are injected 
per unit area per unit time, the mass flux m, consisting of liquid 
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Fig. 1 Radius vs. time data for vaporizing fuel droplets (7, = 
3000 K, P, = 300 psia) 
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Fig. 2 Nondimensional diameter vs. time dependence of vaporiz- 
ing fuel droplets based on data in (7) 
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and gas phase fluxes m, and mz,, is 
m =m + m, = (7/6)D3p.N [2] 


where p; is the liquid density. The gaseous mass flux at any 
station along the combustion chamber 


mM, = — Mi = [3] 


is readily calculated from given vaporization rate laws (e) 
under the assumption of full entrainment (f) of liquid droplets 
in the product gas velocity uy. 

Full entrainment implies that the droplet flux N is constant 
at any section of the one-dimensional chamber, hence m; for 
the species is 


D 3 
= 5 = ti (3) [4] 
where D is the local droplet diameter. 
It is convenient in the calculations to employ the time 
elapsed after injection ¢ as a parameter which connects u, and 
the coordinate x downstream of the injector 


x -f; at [5] 


Thus, if the droplet history is given by the vaporization law 
D = D(t), Equation [3] yields immediately 
6) 
m Dé uy 


where uy is the final gas velocity attained after an interval 
ty corresponding to the droplet lifetime defined by D(t;) = 0. 
In accordance with the assumed ideal gas law (b), u; is de- 
termined from 

m mRT, 


Po P, 


Fig. 3 Nondimensional gas velocity vs. coordinate in monopro- 
pellant combustion chamber with uniform injected droplet size 
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Fig. 4 Sketch of incremental relations for elementary zone of 
combustion chamber model 
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The velocity u, in Equation [6] is attained at the coordinate 
determined from Equation [5] 


r= Uy t- Ds 


For the special form D(é) in Equation [1], we have in non- 
dimensional form with r = t/t; = Kt/D;? asa parameter 


1 (1-7) (9 
us 
x 5 2 
ac (1 T) } 
where 
3 
zy = 5 {1 1] 


is the distance traversed by the droplet in its lifetime 7 = |, 
i.e., the combustion length. These results are illustrated in 
Fig. 3 showing the dependence of the velocity on coordinate 
along the monopropellant combustion chamber. With suii- 
able modifications, the prototype parametric results in Equ:- 
tions [9 and 10] prove to be useful in calculating velocity an] 
coordinate results for bipropellant systems with inclusion of 
the droplet size distribution effects, as will be shown. 


Bipropellants 


In this case we distinguish between the histories of fuel and 
oxidant droplets Dr(t) and Do(t) as they may be character- 
ized by different vaporization parameters Kr and Ko, respec- 
tively. Since the chemical reaction is assumed to occur in « 
stoichiometric mixture ratio, M,, = (O/F)., the presence o/ 
excess unreacted propellant sharing in the local heat release 
may, in general, lead to spatial dependence of the steady-state 
temperature 7’,. For the sake of simplicity, we shall assume 
further that the local vaporization rates always provide excess 
oxidant, so that the heat release rate is expressible in pro- 
portion to the local fuel vaporization rate. 

Consider now the steady-state gas mass balance and rela- 
tions for the elementary zone between the planes X and X + 
dX, as sketched in Fig. 4. The incremental mass flux, con- 
sisting of both reacted products and unreacted excess gas 
(oxidant in this instance), is 


dm, = + diigo = [12] 


The incremental enthalpy flux is referred to the heat value of 
the fuel Hy, under standard conditions at temperature JT = 
Te 


= (Hr, » — AHp) dtitg, — AHodtiy.g [13] 


where the enthalpy corrections AHr and AH are so deter- 
mined as to allow for possible nonstandard initial conditions 
of the fuel and oxidant, respectively. (These corrections are 
small compared with Hy,z for most propellant combinations.) 

Integration of Equations [12 and 13] yields the continuity 
equation 


M, = Mg, r + = [14] 
and the energy equation 
— = (Hr — AHr)my,, AHom,,o [15] 


wherein the integration constant cp,s T's has been introduced 
to refer the enthalpy of the product gas to the standard tem- 
perature 7's. With p, eliminated by means of the equation 
of state 


Pp, = P,/RT, [16] 
Equations [14 and 15] yield the following results for gas 
ARS JouRNAL 


of 
: 
T | 
| ast t Pg 4 
) 4 4 j 
° 02 04 08 
4 
| | 
: 


ordinate 


[8] 


in non- 
er 


[9] 


[19] 


(11) 


r= 
uted in 
‘dinate 
h suit- 
Equ:- 
ty 
sion of 


el and 
‘acter- 
espec- 
rina 
nee of 
elease 
-state 
sume 
2XCeSS 
pro- 


rela- 
X + 
con- 
gas 


eter- 
bions 
3 are 
) 


uity 


[14] 


[15] 
iced 
em- 
tion 
[16] 
gas 


VAL 


velocity and temperature 


U, = F — GoM, 0 [17] 
where 
R 
(Hr — AHr + Ts) [19] 
Jo = — T's) [20] 


The changing chemical composition within the bipropellant 
chamber produces changes in cp, R, gr, go, Which, on the basis 
of simplified reaction kinetics, can be expressed as functions 
of the auxiliary variables m,, » and ti, 9 (see Appendix). 
T iese variables, in turn, are implicitly dependent on the time 
i lapsed after injection, in analogy with the monopropellant 
case 


Dr (t 

tity, = the (: [21] 

Tig, 0 = tho (: ~ [22] 


where mp and mo designate total mass flux of the injected 
species. Thus, with ¢ as a parameter Equations [17, 18] and 


X =f, U,dt [23] 


constitute the parametric relations connecting U, and T, 
with the coordinate X downstream of the injector. 

In order to express the salient features of the bipropellant 
combustion field in simple analytic form, we now assume 
constant gas properties c, and R. The coefficients gr, go are 
then also constant, in view of Equations [19 and 20]. Intro- 
ducing the terminal values attained upon complete vapori- 
zation of the fuel at t = tr > to 


Uy = — Jomo 
[25] 
tho + the 
> 
X; = f, — o)dt [26] 


and employing the dimensionless notation 


= 
~ 


mo/Mr Q = qo/qFr 
nr(l) = my, Noll) = o/tmo 


we arrive at the parametric relations for assumed constant 
gas properties 


Udt) ne — QMno 

1-QM 
Tot) _ ( 1+M 129) 
Ty nr +Mno/\ 1 QM 


t 

Jo. QM np) a 
t 

Xr (ne at 


[30] 


For liquid oxygen/hydrocarbon systems, the overall mixture 
ratio M ranges from 2 to 3, whereas the enthalpy correction 
ratio Q is of the order of 10~?, tending to make QMno < nr, 
even though, due to the greater volatility of oxygen, no > nr 
in the early stages of droplet vaporization. With QM <1, 
Equations [28-30] indicate that the principal effect of early 
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Fig. 5 Influence of relative vaporization rates on temperature 
and velocity distributions in bipropellant combustion chamber 
(mixture ratio J = 2.5) 


vaporization of the oxidant is on the temperature distribu- 
tion T, vs. X, whereas U, vs. X is determined primarily by 
the fraction of fuel vaporized nr. It is to be noted that 
parametric relations of the type Equations [28 and 29] ex- 
pressing U, and 7, in terms of the vaporized propellant frac- 
tions nx and no may be derived, in general, from aerothermo- 
dynamic principles. However, the connection between X and 
the vaporized fractions can be determined only by taking into 
account the kinetics of gas generation, i.e., the motion of 
vaporizing droplets. The assumption of equal droplet and 
gas velocities (f), whereby nr and yo are interrelated by a 
single parameter t, yields the particularly simple expression 
for X in Equation [30]. The present treatment can be 
readily extended to account for the influence of propellant 
droplet velocity lags under representative frictional entrain- 
ment laws. 

A graphical illustration of the gas flow conditions de- 
scribed by Equations [28 to 30] is illustrated in Fig. 5 based 
on the vaporization law 


Dy? = Dp, — Krt [31] 
= Do Kot [32] 


with the assumption that the oxidant droplet lifetime to = 
Dé,i/Ko is not greater than tr = D7, i/Kr. The influence 
of relative vaporization rates is illustrated for three different 
values of the droplet lifetime ratio, to/tp = 4, % and 1, ata 
given overall mixture ratio M = 2.5and Q = 0.01. 

It is seen from Fig. 5 that the changes in to/tr have a far 
greater influence on the temperature distribution (dotted 
curves) than on the velocity development (solid curves) in 
the combustion chamber. This is readily understood on the 
following physical grounds. Reduction of the ratio to/tr 
increases the oxidant-to-fuel vaporization rates, and, since 
the heat release depends primarily on the vaporized fuel, 
less heat will be released per unit mass of the vaporized 
mixture. The consequent reduction of temperature at con- 
stant chamber pressure results in increased gas density. 
Thus, the greater mass flux due to increased oxidant vapor- 
ization with lower temperature and greater density in the 
chamber can be accommodated without large changes in the 
gas velocity. 

It is to be noted that the monotonic increase of 7, vs. X 
shown in Fig. 5 is predicated on the assumption that, in the 
oxidant-rich gas stream considered here, tr > to. If tr < to, 
a temperature peak in excess of 7’; (see Eq. [25]) is attained 
at the end of fuel vaporization. More generally, it is readily 
seen from Equation [29] that temperature maxima and 
minima along the vaporization rate limited bipropellant 
combustion chamber can occur if 

aT, dnr dno [33] 


dt 
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~ 

= 

= 

~ 

[12] 
[13] 


Finally, we note that at injector, X = 0, nr ~ 0, no > 0, 
the temperature ratio in Equation [29] is calculated by a 
limiting process: t > 0, no/nr — tr/to, which yields 


_ ( 1+M QM site) (34) 
\14+M te/to 1 — QM 


The Influence of Droplet Size Distributions 


If the number of droplets injected per unit area per unit 
time N contains a spectrum of initial droplet sizes specified 
by the distribution function n(D;), we have 


aN = n(D;) dD; [35] 


as the component of the spectrum in the size range from D,; to 
D;+ dD;. Let the time history of droplets in this size range 
be designated by D(i, t) with lifetime t,;,; determined from 
D(i, t;,:) = 0. Let, further, the identity of droplets be fixed 
upon injection in the flux component dN defined by Equation 
[35]. As the droplets in this group move along the chamber 
with the velocity of the product gases, the droplet flux re- 
mains constant at the value dN, while the associated elemen- 
tary liquid mass flux is 


dm, = dN [36] 
Hence, for the model under consideration with total mass flux 
Di = Dn 
f Dt n(D,) dD; [37] 
6 Di =0 


the liquid mass flux surviving up to a time ¢ after injection is 
given by 


Dm 


m Dm 
n(D;)dD; 
0 


[38} 
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Fig. 6 Comparison of droplet size spectra for constant mass 
flux m and maximum diameter D,. [m/Dm = (3/6) p, Dm?) 
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Fig. 7 Velocity vs. time elapsed after injection for various in- 
jected droplet spectra with maximum droplet diameter D,, 
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Table 3 Specification of droplet size distribution functions 


pt 
Spectrum a b 6 m CarDm 
I 1 4 630 
II 1 1 30 
III 4 1 72 


where D,, is the maximum diameter in the injected spectrun:. 
In exhibiting the role of droplet size distribution, the surviving 
liquid fraction m,/m in Equation [38] replaces the fraction 
D*(t)/D;° appearing in the previously derived expressions fcr 
uniform injected droplet size (Eqs. [6, 8, 21, 22]). 

The droplet size distribution produced by the atomization 
process is generally described by empirically determine | 
distribution functions sufficient to specify such basic featur: s 
of the spray as mean droplet size, skewness and dispersio: 
about the mean. In the illustrative calculations below w: 
shall employ for each propellant species a polynomial repr« - 
sentation of the injected distribution function in the form 


n(D;) = CyDi*(Dm — Di)? [34 
where a, b are integers, and C,» is a constant related to the 


injected mass flux of the species by 


Dm 
0 


3 
ll 


Equation [40] is adequate for approximate representation 
of some actual distribution functions, such as, for example 
the Nukiyama-Tanasawa equation ((1), pp. 4 to 53). More- 
over, it is of great convenience in treating a broad class of 
possible distribution functions without entailing undue 
mathematical complexities. 

The influence of the injected droplet spectrum is illustrated 
for three distribution functions (Spectra I, II and III in Table 
3) having peak droplet populations respectively at 4, } and 
% of a reference maximum diameter Dm. Graphical com- 
parison of these spectra referred to fixed values of m and D,, 
is shown in Fig. 6. 


Monopropellants 


Applying the distribution functions to a monopropellant 
obeying the vaporization law (Eq. [1]) 


Di, t) = D2? — Kt [41] 


we obtain from Equation [38] 


1 


m [ere — dé 
0 
where 
= Kt/D,,? 
= D,;/D,, 


The integrals required here are readily evaluated from stand- 
ard formulas (13) for 
bed 
where p is an integer. The results showing nondimensional 
gas velocity vs. 7’ 
Uys m m 


[43] 


are given in Fig. 7. A second integration with respect to 
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time indicated in Equation [5] leads to 


Uy, dr’ 
= = Jy [44] 


Ty; 1 
[mar 
0 


in which the required integrals are also readily evaluated by 
standard methods. 

Sased on Equations [43 and 44] with 7’ as a parameter, 
the calculated results for the velocity dependence in a mon- 
propellant combustion chamber are shown in Fig. 8. This 
figure, in effect, represents the efficiency vs. length for the 
combustor model considered here. In the high efficiency 
recime, u,/uy ~ 0.9 significant reductions in length are indi- 
cated by improvement of atomization. However, since 
complete combustion depends on vaporization of the largest 
inected droplet (Dm is the same for the distributions con- 
si cred here), failure to reduce the max mum droplet size 
under enhanced atomization results in increasing the com- 
bustion length required for 100 per cent efficiency. 


Bipropellants 

In the bipropellant combustion field given by Equations 
[25 to 30], the role of the injected droplet size spectra is 
implicit in the vaporized fractions nr and no defined by 
Equation [27]. The latter quantities can be alternatively 
defined in terms of the surviving liquid fractions 


9 = 1 — (mi/m) [45] 


and, hence, can be calculated by application of Equation [38] 
to each species. Insertion of the results for nr and no 
based on Equations [38 and 45] into Equations [28 to 30] 
yiclds the parametric equations connecting U,, 7, and X corre- 
sponding to droplet size distributions nr(Dr, :), No(Do, i) and 
vaporization laws Dp(i, t), Do(t, t). 

In view of the role of nr and no in the parametric equa- 
tions, the influence of bipropellant droplet size distributions 
may be summarized as described below. The conclusions 
in the next paragraph refer to the assumed oxidant-rich con- 
ditions in the derivation of the prototype solutions in 
Equations [28 to 30].4 

Since for typical bipropellant combinations @ <1,MQ <1, 
the effect of the droplet spectrum on U,/U, vs. X/X + appears 
primarily through the modified nr. However, the modified 
no has a strong influence on the 7',/7'; vs. X/X,; distribution. 

The terminal values U; and 7, are determined independ- 
ently of vaporization by given propellant injection rates, 
heat released per unit mass of the overall mixture, and the 
pressure level in the combustion chamber. The combustion 
length X, depends primarily on 7p, i.e., on the vaporization 
rates and droplet size spectrum of the fuel. Graphical illus- 
trations of these results are shown in Fig. 9 for oxidant and 
fuel droplet spectra corresponding to the droplet spectra I 
and III, in Table 3, with assumed lifetime ratio tg/tr = 4 and 
mixture ratio M = 2.5. 


Flameholding Stability 


The previously derived results were obtained on the 
assumption that vaporization is followed by prompt local 
mixing and chemical reaction. The latter can occur only 
if the locally vaporized propellants mixture ratio is within 
sufficiently close limits about the stoichiometric ratio and the 
local temperature is sufficiently high, dependent on the pro- 
pellant combination. For oxygen-inert-hydrocarbon systems 
in particular, fundamental experiments indicate that flame 
propagation can be achieved only if the resultant flame 


* Analogous results are readily obtained within the scope of 
the present analysis if the local vaporization rates.provide excess 
fuel everywhere, or if transitions occur from oxidizer-rich to fuel- 
rich vaporization rates along the combustion chamber. 
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Fig. 8 Velocity vs. coordinate in chamber for various injected 
droplet spectra with maximum droplet diameter D,, 
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Fig. 9 Influence of inj2cted droplet size spectrum on tempera- 
ture and velocity distributions in bipropellant combustion cham- 
ber (lifetime ratio (9 /tr = 1%; mixture ratio M = 2.5) 


temperature is above ~1200 K (14). If we introduce such 
a minimum temperature as a threshold temperature at which 
reaction is achieved in the bipropellant combustion chamber, 
we arrive at a stability condition necessary for continuous 
ignition of the propellants (i.e., flameholding) in the low ve- 
locity region near the injector face 


T(0) => Tmin [46] 
In view of Equation [29] the flameholding stability condi- 
tion can be expressed as 
T, limit 
t—>o (nr + Mno) 


(nr — QM 
QM) = T min [47] 


where the 7’; is the terminal temperature, dependent on the 
overall mixture ratio M = mo/m,y. For the particular case 
of uniform oxidant and fuel droplet size with lifetimes to and 
tr, the stability condition yields 


(1+M) (1 — QMtr/to) 
(1+ Mip/to) (1 — QM) 


In application to oxygen/hydrocarbon systems with assumed 
Tmin = 1200 K and typical data M ~ 2.5, T, ~ 3500 K, 
Equation [48] implies that the stability limit is reached if 
tr/ to ~ 3. 


> Twin [48] 


Discussion 


The idealized conditions of the vaporization rate controlled 
combustor model are approached when mixing and chemical 
time lags, present in actual combustion chambers, are at a 
minimum. The limitations and possible scope of application 
of the theoretical results are briefly discussed on the basis of 
relevant experimental data. 
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Monopropellants 

By injection of fuel and oxidant in chemically homo- 
geneous or physically premixed liquid jets, as in the case of 
monopropellants, the mixing lags are minimized; high flame 
temperature (~3000 K) and associated reduction of chemical 
lags are additional prerequisites for the emergence of vaporiz- 
ation as a rate controlling step. A survey of available mono- 
propellant combustion data (15) reveals that the latter pre- 
requisites are not satisfied, inasmuch as currently familiar 
low performance monopropellants generally undergo slow 
gas phase chemical reactions. It would be interesting to as- 
certain the applicability of the vaporization model to re- 
cently developed high performance monopropellants (16). 
In particular, the vaporization model implies that the product 
gas in monopropellant combustion is at a constant temper- 
ature (flame temperature); this result could be most simply 
ascertained by noting to what extent the exhaust gas tem- 
perature is independent of combustion efficiency determined 
under experimental conditions. 


Bipropellants 

In regard to bipropellant combustion, if vaporization rate 
is the slow step in the combustion process, unreacted gaseous 
propellants of both species cannot coexist in the product gas 
flow. However, the latter may be diluted by the presence 
of an excess, unreacted gas phase propellant. As the re- 
sults summarized in Equations [28 to 30] and illustrated in 
Figs. 5, 8 and 9 indicate, the velocity or efficiency distribution 
is primarily determined by the droplet spectrum and vapor- 
ization rate of the slowly vaporizing species. The principal 
effect of varying the droplet size spectrum of the more rapidly 
vaporizing species appears in the temperature distribution 
which, in general, is depressed by the presence of excess un- 
reacted propellant gas. 

In Priem’s (2,3) combustion model it had been assumed 
that hydrocarbon vaporization is rate controlling while the 
product gas remains constant at the flame temperature 7’. 
Although the gas temperature variations influence the vapor- 
ization rates to some extent, in view of the results summar- 
ized in the preceding paragraph, Priem’s assumptions appear 
satisfactory for the calculation of gas velocity or combustion 
efficiency vs. combustor length. 

Experimental data on temperature and chemical distribu- 
tions have been obtained by sampling of the gas flow with 
pneumatic probes in liquid oxygen-ammonia (4) and nitric 
acid-hydrocarbon (5) combustion chambers. The investiga- 
tions reported were done under comparable pressure (P, ~300 
psia) and thrust (ca 1000 Ib) conditions with a variety of 
hydraulic injection patterns. Evidence of mixing delay is 
shown in both investigations. The temperatures reported 
appear to be generally well above the ignition temperature 
and increase progressively in the axial direction from the in- 
jector toward the exhaust nozzle, but radial variations are 
observed associated with the hydraulic injection patterns. 
In the oxygen-ammonia combustion investigated by Baker 
(4), coexistence of unreacted fuel and oxidizer did not occur 
except in some low temperature gas pockets associated with 
poor mixing. On the other hand, the nitric acid-hydrocarbon 
combustion data obtained by Trent (5) shows the existence 
of nonequilibrium concentrations of nitric oxide at high 
temperatures, suggesting, as a rate controlling step, the 
chemical process whereby the intermediate nitric oxide is 
converted to product gases. Similar investigations of oxygen- 
hydrocarbon combustion might be expected to reveal the role 
of carbon monoxide conversion as a rate controlling step 
under sufficiently rapid mixing conditions. 

From the preceding considerations it is concluded that a 
fundamentally adequate account of rocket combustion phe- 
nomena cannot be based solely on the vaporization rate con- 
trolled model. However, the latter model provides a ra- 
tional approach to the determination of technically useful 
correlations between rocket combustion and vaporization 
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data, particularly in high performance combustors. High 
performance, characterized by high volumetric reaction rates 
in the combustion chamber, is achieved in practice with 
highly reactive chemical propellants under adequately rapid 
mixing conditions favoring the emergence of vaporization 
as a rate controlling process. 


Appendix 


The influence of variable gas properties is treated by con- 
sidering the presence of product gases and excess gas (oxi- 
dant in the example treated) with molecular weight oN), and 
Mo, respectively. With reaction rate occurring at assumed 
constant stoichiometric ratio M,,, the species referred to have 
the flow rates 


Mopr = Mg, + Ms) [A-1] 
= Mo,0 — Mg, FMst [A-2} 


Expressing the additivity laws for c, and R in terms of specie: 
flow rates (to which the densities are proportional by virtu: 
of the common transport velocity U,), we have 


oO , Cp pr 
= + —— Mor A-3 
a 9,0 P [ 
mo, Oo Mor 
mk = Ral — + — A-4! 
( Mo 
where 
C, = molal specific heat at constant pressure 


Rox = universal gas constant 


Elimination of mp, and mm, 9 from Equations [A-3 and A-4] 
by substitution from Equations [A-1 and A-2] yields the ex- 
pressions for c, and R in terms of the auxiliary variables m,, r 
my, 0, and hence in terms of the time parameter ¢ implicit 
iN tg, Fr, ™g,o (see Eqs. [21 and 22]). In this manner, 
variable gas properties c, and R can be treated within the 
system of parametric Equations [17, 18 and 23] which con- 
nect U, and T, with the coordinate X. 

In particular, we note that the preceding modifications lead 
to a simple expression for the enthalpy ratio (see Eq. [29]) 


eA(T,)T, ( 1+M [A-5] 
nr + Mno/\ 1— Q;M 


from which the temperature 7, may be evaluated by suitable 
iterative procedures if temperature dependence of the molal 
specific heats C,(7,) is taken into account in the calcu- 
lation of C,(7,). Application of Equation [A-5] to the 
flameholding stability limit considered in the example shows 
that, because of lowered gas specific heats at Tmin = 1200 C, 
the stability limit is reached if the fuel-to-oxidant lifetime 
ratio is increased to tr/to ~ 3.5. 


Nomenclature 


a,b = exponents in polynomial representation of dimen- 
sionless droplet size distribution, 
gas specific heat at constant pressure, cal/gm K 


Cc => 

4 = injected propellant mass flux, gm/cm? see 

Mg = gas phase mass flux, gm/cm? sec 

m = liquid phase mass flux, gm/cm? sec 

mp = total fuel mass flux, gm/cm? sec 

mo = total oxidant mass flux, gm/cm? sec 

n = droplet size distribution function, no./em? sec 
q = velocity coefficient, cm*/gm 

t = time elapsed after injection (droplet age), sec 
ty = droplet lifetime, sec 

lp = maximum fuel droplet lifetime, sec 

lo = maximum oxidant droplet lifetime, sec 

us = terminal gas velocity (monopropellant), cm/sec 
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Ug = velocity in chamber (monopropellant), cm/sec pr = product gas 
x = coordinate in chamber (monopropellant), em B = refers to standard conditions 
zy = combustion length (monopropellant), cm F = fuel 
Ca = rate constant in polynomial representation of in- O = oxidant 
jected droplet size distribution, no./sec em!¢+6+3) 
Cy = molal specific heat at constant pressure, cal/mole K Po 
D = instantaneous droplet diameter, em 
D; oe injected droplet diameter, cm 1 Battelle Memorial Institute Personnel, ‘‘Injection and Combustion of 
De = maximum injected droplet diameter, cm Liquid 
. ais ¢ 2 riem, R. J., ‘‘Propellant Vaporization as a Criterion for Rocket- 
Hr, — yor of fuel referred to standard temperature Engine Design; Calculations of Chamber Length to Vaporize a Single n- 
a, cal/gm : Heptane Drop,” NACA TN 3985, 1957. 
AH», AHo = enthalpy corrections for fuel and oxidant, cal/gm 3 Priem, R. J., ‘‘Propellant Vaporization as a Criterion for Rocket 
K = vaporization rate constant, cm?/sec Engine Design; Calculations Using Various Log-Probability Distributions of 
L* characteristic combustor length, Temperature Surveys of Ammonia- 
M = mo/mp = overall mixture ratio, dimensionless Oxygen Rocket Motor Combustion Chambers,’’ JetT Proputsion, vol. 25, 
AI = molecular weight, gm/mole 1955, p. 217. 
N = inie 3 a 2 5 Trent, C. H., ‘Investigation of Combustion in Rocket Thrust 
5 injected droplet flux, no./cm? sec : Chambers,’’ Ind. Engng. Chem., vol. 48, 1956, p. 721. 
Ig = combustion chamber pressure, gm/cm’ : 6 Bolt, J. A. and Saad, M. A., ‘‘Combustion Rates of Freely Falling 
Q = go/qr = ratio of velocity coefficients, dimensionless Fuel Droplets in a Hot Atmosphere,” in ‘‘Sixth Symposium on Combus- 
Ty = temperature at the end of combustion zone, K tion,” Reinhold, 1957, p. 717; see also Kumagai, 8. and Isoda, H., ‘‘Combus- 
= local gas temperature K of Chamber,” in “‘Sixth Symposium on Com- 
ustion,” Reinhold, 1957, p. 726. 
Uy = terminal velocity (bipropellant), cm/sec 7 Penner, S. S. ‘‘On Maximum Evaporation Rates of Liquid Droplets 
Uy = velocity in chamber (bipropellant), cm/sec in Rocket Motors,’’ JouRNAL Or THE AMERICAN Rocket Society, vol. 23, 
x coordinate = chamber (bipropellant), — 8 Besserer, C. W., ‘‘Missile Engineering Handbook,"’ van Meetwnnd, 
Xs = combustion length (bipropellant), em Princeton, 1958. 
r = m,/m = vaporized fraction of propellant species, 9 Spalding, D. B., ‘‘Some Fundamentals of Combustion,’’ Academic 
dimensionless Press, New York, 1955. 
= D,/Dmn = li . 1 initial d let diamete 10 Goldsmith, M., ‘‘Experiments on the Burning of Single Drops of 
= i/'m ae nondimensional initia. ropiet diameter, Fuel,” Jer Propuusion, vol. 26, 1956, p. 172. 
dimensionless 11 Hall, A. R. and Diederichsen, J., ‘‘An Experimental Study of the 
Pg = gas density, gm/em* Burning of Single Drops of Fuel in Air Pressures up to Twenty Atmos- 
pl liquid density, gm/cm* on Combustion,’ Williams and Wilkins, 
tT = Kt/D? nondimensional time (uniform initial drop- 12 Berman, K. and Logan, 8. E., ‘‘Combustion Studies with a Rocket 
let size), dimensionless Motor Having a Full-Length Observation Window,’’ JouRNAL OF THE AMERI- 
h = Kt/D,? nondimensional time (initial droplet spec- can Rocket Society, vol. 22, 1952, p. 78. 
( P 13. Dwight, H. B., ‘‘Tables of Integrals and Other Mathematical Data,’’ 
trum with maximum diameter D,,), dimension- 
less 14 Egerton, A. C., ‘Limits of Inflammability,”’ in ‘‘Fourth Symposium 
on Combustion,”’ Williams and Wilkins, Baltimore, 1953, p. 4. 
. ‘ 15 Altman, D. and Penner, S. S., ‘‘Combustion Processes,’’ Vol. 2, in 
Subscripts ‘‘High Speed Aerodynamics and Jet Propulsion,’’ Princeton University 
Press, 1956, p. 489. 
Mis final, compl tely vaporized 16 Clark, J. D., ‘‘Dark Horse in the Propellant Race,’’ ASTRONAUTICS, 
g = gas vol. 3, no. 9, 1958, p. 34. 
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Correlation and Prediction of Flame 
Properties With Special Reference to 


MELVIN GERSTEIN! 


Liquid Hydrazine 


NASA, 
Lewis Flight Propulsion Laboratory 
Cleveland, Ohio 


A large background of experimental data, analytical studies and correlations exist for hydrocarbon 
air flames. Much less work has been done with propellants of interest in the rocket field. It is 
possible, however, that some of the experimental and analytical techniques developed for flames of 
fuels burning in air will be useful when applied to the fuels and oxidants used in rockets. This 
paper presents a summary of some of the predictions and correlations which have been reported. A 
brief discussion of the application of some of these correlations to the liquid hydrazine decomposi- 
tion flame is also included. Such a flame, using a propellant with relatively simple flame kinetics 
and which can support a flame when introduced either as a liquid or a gas, may well serve to bridge 


the gap between gas phase combustion studies and combustion of c P 


and solid monopropellants. 


Flame Velocity 


| Ppa INAR flame velocities have been calculated for several 
systems in which the chemical reactions can be repre- 
sented by relatively simple kinetics. Typical of these are 
the hydrazine and ozone decomposition flames and the hydro- 
gen-bromine flame. It is particularly fortunate that chem- 
istry has, in this case, favored rocket combustion interests, 
hydrazine and ozone being of obvious importance, whereas the 
hydrogen-bromine flame may well be a better introduction 
to the hydrogen-fluorine flame than to the much more com- 
plicated hydrogen-air flame. 

The equations used in predicting flame velocity involve 
straightforward use of the equations of conservation for a 
chemically reacting system. An example of the results ob- 
tained from such calculations is compared with experimental 
data in Table 1. 


Table 1 Flame velocity, cm/sec 


System Calculated Experimental 
H,-Brz (40 per cent Bro, 34 20-32 
60 per cent He) (1) 
O, (25 per cent Os, 47 55 
75 per cent O, (2) 
N2H, (3) 184 185 


Considering the approximations required to obtain solu- 
tions and the lack of good high temperature kinetics and gas 
property data, the agreement is good. 

Attempts to calculate the laminar flame velocity of systems 
having more complicated reactions have not been very suc- 
cessful, although it has been possible to predict changes in 


Presented at the ARS 13th Annual Meeting, New York, N. Y., 
Nov. 17-21, 1958. 
1 Assistant Chief, Propulsion Chemistry Division. Member 
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d d ph » such as liquid 


flame velocity produced by changes in initial mixture tem- 
perature and to some extent fuel-air ratio and pressure. 


Quenching Distance 


Recently, quenching distances have been calculated for the 
same flames (Table 1) which have been the subject of flame 
velocity calculations. Although quenching distances may 
be of limited interest in most practical combustion systems, 
it will be shown later that the quenching distance is a rela- 
tively easily measured property, which can be used in con- 
junction with the flame velocity to obtain information related 
to the chemical reaction rates in flames. A simplified model 
of flame quenching suggests that the flame is quenched when 
sufficient energy has been removed, so that the chemical 
reactions are no longer self-sustaining. The heat loss can be 
calculated using conventional heat transfer methods, and the 
rate of chemical reaction can be obtained when the kinetics are 
known. A further criterion is needed, the quantity of energy 
which must be abstracted before the flame reactions are no 
longer self-sustaining. The criterion used in (4)? was the 
temperature of the inflection point in the temperature-dis- 
tance curve through a plane flame. The results obtained 
in this way are shown in Table 2. 

For the H.-Br, flame, the agreement is excellent, not only 
for the calculated quenching distance at 1 atm shown in the 


Table 2 Quenching distance, cm 


System Calculated Experimental 
H.-Brz (55 per cent He, 0.66 0.6-0.7 
45 per cent Bre) 
O; (25 per cent Os, 0.35 0.25 
75 per cent Oz) 
0.05 0.12 


2 Numbers in parentheses indicate References at end of paper. 
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table, but also the pressure dependence is calculated to be 
Q.D. « P-** 
which can be compared to the experimental value 


« 


The agreement for ozone and hydrazine is not as good. 
‘The experimental data for ozone are open to some questions, 
so that the validity of the comparison is not certain. The 
only experimental quenching distance found for hydrazine in- 
volved the limiting tube diameter for flame propagation 
through a column of liquid hydrazine. The experimental 
\alue is 2.4 times the calculated value. The reaction rate 
used in the calculation involved gaseous hydrazine. If the 
-ombined evaporation and reaction required for the liquid 
column is a slower process, then the deviation is in the 
proper direction. The liquid hydrazine flame will be dis- 
cussed in more detail later. 

Those who have attempted to calculate quenching dis- 
iances for flame reactions more complex than those involved 
in Table 2 have been troubled by the same problem mentioned 
in connection with the calculation of flame velocity, namely, 
the lack of suitable knowledge about the chemical kinetics. 
The use of a simplified model has been reasonably successful 
in predicting the effects of pressure, temperature and tube 
veometry on quenching distance. The relationship between 
reaction order and pressure dependence of quenching dis- 
tance has been useful in obtaining approximate reaction orders 
for the overall process occurring in a flame. The failure of 
these reaction orders to predict the effects of changes in fuel- 
air ratio (5) indicates that a single overall reaction is not ade- 
quate for the description of hydrocarbon-air flame properties. 


Relation Between Flame Velocity and Quenching 
Distance 


Dimensionally, flame velocity may be considered propor- 
tional to the square root of the ratio of a transport property a 
and a reaction time 7; the quenching distance may be con- 
sidered proportional to the square root of the product of a 
and 7, where 7 is related to the reaction rate W and initial fuel 
concentration Ao by the relation 


W= Ao/T 


The ratio of flame velocity U, to quenching distance D 
should then be given by 


D Ao 


and the product U;D « a. 

Reaction rates calculated from the first expression are found 
to give reasonable reaction orders and rates comparable with 
those deduced from other flame properties (6). 

It is shown in (7) that the second expression is in reasonable 
agreement with experiment. Large deviations were found 
for hydrocarbon-nitric oxide flames and the ethylene oxide 
decomposition flame. Whether these represent failures of 
the correlation or experimental errors is not certain. 


Liquid Hydrazine Flames 


The rate of burning of a liquid column of hydrazine varies 
lirectly with pressure and the limiting tube diameter in- 
versely with pressure for sub-atmospheric pressures (8). The 
ratio U/D thus varies with the pressure squared at low pres- 
sures, indicating that a second order reaction is involved. 
Since the rate of the gas phase decomposition of hydrazine is 
second order (9), it is suggested that the gas phase reaction 
rate may be a controlling factor at low pressures, a conclusion 
lrawn in (8) from the burning rate study alone. The inverse 
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pressure dependence of the limiting tube diameter is also 
consistent with a second order reaction. 

Using the relationship between U/D and W developed in 
(7), and including a factor for the ratio of moles of product 
to moles of reactant, one can estimate a reaction rate of 
0.036 mole/cm*-sec, which is three-fifths of the value neces- 
sary to predict the quenching distance in Table 2. Since 
the quenching distance depends on the square root of the 
reaction rate, the actual error in predicting quenching dis- 
tance would be about 30 per cent. Note that while the 
pressure dependence of both burning rate and limiting tube 
diameter suggests that the gas phase reaction rate may con- 
trol, the estimated value of the reaction rate for liquid hydra- 
zine is an order of magnitude slower than that calculated for 
the gas phase reaction (4). This suggests that an additional 
process depending only slightly on pressure may be involved. 

It has been suggested that the burning of solid propellants 
is a two-step process (10) involving the gasification of the 
propellant followed by a gas phase reaction. If the same con- 
cept is applied to the liquid hydrazine flame, the additional 
process required, evaporation of the liquid, is only slightly 
dependent on pressure (11). 

A very simple model produces some interesting results. 
Consider that the rate of disappearance of the liquid U, 
which is also the burning rate, is determined by the rate of 
heat transfer from the gas phase to the liquid surface 


poAH, AX 
where 
rN = mean thermal conductivity 
Po = liquid density 
AH, = heat absorbed by the liquid 


AT/AX = assumed linear temperature gradient 


Since the flame temperature is relatively independent of 
pressure, the change in gradient required to make U vary 
linearly with the pressure is that AX vary inversely with 
pressure, where AX is now the thickness of the reaction zone. 
For a second order reaction, the thickness of a gas phase 
flame does vary approximately as 1/P (12), which is then 
consistent with the low pressure experimental data. At 
higher pressures (1 to 8 atm) the burning rate of liquid hydra- 
zine varies roughly as P'’, while above 8 atm the burning 
rate remains constant or may even drop slightly (13). Since 
the gas phase reaction rate varies as P? and the evaporation 
rate is relatively independent of pressure, it is not unreason- 
able that at low pressures the rate is controlled by the gas 
phase reaction and at high pressures by evaporation. Ina 
simple way, one can assume that the flame thickness is made 
up of two factors, AX = ho + hi/P. where ho is controlled 
by heat and mass transfer and is relatively independent of 
pressure, while h, is thickness of the gas phase reaction zone 
at 1 atm and for a second order reaction varies as 1/P. The 
data for the burning rate of liquid hydrazine from (8 and 13) 
are plotted in Fig. 1. The curve drawn through the data was 
calculated from the relation 


0.0017 
ho + hi/P 


The values of Ao and h; were obtained empirically from ex- 
perimental values of U at 3 and at 16 atm. It is seen that 
the simple model produces a curve of the proper shape. 

The addition of small amounts of water should not change 
the evaporation characteristics of the liquid mixture but 
should affect the rate of the gas phase reaction. Examina- 
tion of Table 3 shows that there is no systematic variation 
of burning rate with water content at high pressure, whereas 
at 1 atm the burning rate decreases as the water content in- 
creases. 


515 


rbon 
It is 
es of 
This 
lL A 
etics 
idge 
quid : 
the 
ame 
nay 
‘ms, 
ela- 
cal 
ica] 
be 
are 
rgy 
no 
the 
lis- 
ned i 
nly 
the 
| 
AL 


O S m 
Ve 
4 
4 a 
3 
© 
006) 
c 
5 
a2 
ee 2 4 681 4 6 610 20 40 60 80100 


2 
Pressure, atmospheres 


Fig. 1 Variation of the burning rate of liquid hydrazine with 
pressure for 98.9 per cent N.H;. Ordinate—burning rate U, 
cm/sec. Abscissa—pressure, atm 


Table 3 Effect of N2H, concentration on burning rate 


Average burning rate (13) 


NeH, concentration 1 atm Plateau 
98.9 0.028 0.077 
97 0.028 0.092 
04-95 0.024 0.098 
94.5 0.022 0.069 


The calculated curve is compared with the experimental 
data of (13) in Fig. 2, assuming that ho is independent of pres- 
sure and obtaining h; from the 1-atm burning rate for a 94 to 
95 per cent N.-H;-H.O mixture. 

The agreement is encouraging. Certainly, the model is 
oversimplified, but it suggests that the liquid hydrazine flame 
cannot be treated in terms of the gas phase reaction alone 
or as a simple evaporation problem neglecting chemical reac- 
tion rate over a wide pressure range. It also appears that 
the rate of the gas phase chemical reaction is not a controlling 
factor at high pressure. 
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Comparison Between Hot Gas Ignition 
and Limit Flame Temperatures 


M. VANPEE? and 
H. G. WOLFHARD® 


Bureau of Mines, 
U.S. Department of the Interior 
Pittsburgh, Pa. 


Hot gas ignition temperatures of a variety of fuels have been determined and compared with limit 


flame temperatures of the same fuels. 


The ignition temperatures were taken under two condi- 


tions: A continuous jet of hot nitrogen flowing into a cold fuel-air mixture, and a continuous jet 


of hot air flowing into cold fuel. 


Limiting flame temperatures were determined by producing a 


diffusion flame in an atmosphere of air diluted with nitrogen. The nitrogen concentration re- 
quired to extinguish the flame was noted and the limiting flame temperature calculated assuming 
stoichiometric combustion. It was found that, with the exception of a few combustibles, a corre- 
lation holds between hot gas ignition temperature and limit flame temperature. The effects of 
various additives, such as bromine, chlorine and methyl bromide were also studied and corroborated 


this conclusion. 


N A SYSTEM containing a combustible gas mixture, 
ignition may occur through contact with hot burned 
gases. One such condition may occur when a laminar stream 
of a combustible gas mixture flows adjacent and parallel to a 
laminar stream of hot burned gases, as treated theoretically 
by Marble and others (1,2).4 Experimental ignition data 
are available for ignition of a virtually quiescent cold mix- 
ture by continuous hot inert jets (3). It is the purpose of 
our research to establish the parameters influencing ignition 
and thus to evaluate the extent to which present theoretical 
treatments are justified. 

In a previous paper (3) it was found that the ignition 
temperatures of a variety of fuels (for ignition by hot inert 
continuous jets) arranged themselves into an order which 
was very different from that known to exist for spontaneous 
ignition or spark ignition. However, there was a close 
relation to hot wire ignition temperatures. Hot gas ignition 
temperatures are much lower than the adiabatic flame 
temperatures of stoichiometric mixtures; they are close, 
however, to flame temperatures for limit mixtures. Thus, 
there is a possibility that the reactions governing limit flames 
are the same as those leading to hot gas ignition. If this is 
true, different fuels should show the same order of temperature 
in both cases. 

A close relation between lean or rich limit flame tempera- 
tures and hot gas ignition temperatures was not anticipated, 
as the fuel-air ratios are very different. A better relation was 
to be expected for the temperature of stoichiometric mixtures 
which had been diluted with inert gas until a limit of flame 
propagation was reached. For convenience, a diffusion 
flame was adopted for comparison, and its limit temperature 
by dilution was determined. The influence of characteristic 
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inhibitors on flame propagation or spontaneous ignition was 
also explored. 


Experimental Setup and Procedure 


Hot Gas Ignition 


The experimental procedure was the same as that described 
earlier (3). Air or nitrogen was heated in a ceramic furnace; 
it formed a continuous jet of hot gas flowing into the outer 
explosive mixture at room temperature. Two new platinum- 
wound ceramic furnaces which were less subject to cracking 
and aging (see (3)) were used, and results were more re- 
producible. Furnace A was 2-ft long and had an inner 
diameter of 10 mm; furnace B was 1-ft long and 4 mm in 
inner diameter. The gas flow for furnace A was 100 ec/sec, 
as this gave the flattest possible temperature distribution of 
the jet near the furnace exit. Furnace B was used for the 
work on additives, with gas flows at 35 cc/sec. 

As in our previous work, no corrections were made for the 
emissivity of the platinum-platinum rhodium thermocouple 
(wire, 0.005-in. diam) owing to the uncertainties of the 
relevant heat transfer numbers. The temperatures recorded 
were measured directly above the furnace exit in the center of 
the jet. The outer chamber containing the room-temperature 
explosive mixture was about 50-cm long and 9 cm in diameter. 
A slow flow of mixture (200 cc/sec) purged the vessel suf- 
ficiently to prevent dilution. Liquid fuels were introduced by 
saturating air with the fuel and admitting further air to get 
the desired concentration. 


Limiting Flame Temperatures 


Limiting flame temperatures were determined by producing 
a diffusion flame around a sintered stainless steel hemisphere 
(4). The fuel entered the flame from within the hemisphere, 
and air flowed slowly upstream against the hemisphere. 
This gave rise to a very stable flame. The air was diluted 
progressively with nitrogen until extinction occurred. This 
resembled the bursting of a soap bubble, as the flame extin- 
guished first at the stagnation point. The oxygen index of 
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the air at the extinction point was noted (oxygen index = 
oxygen in air divided by total amount of oxygen and nitrogen), 
and the limiting flame temperature calculated, assuming 
stoichiometric combustion. The burner system employed 
could be used for both gaseous and liquid fuels without 
modification. 

All gases used in this study were taken from cylinders 
(Matheson Co. cep grade). Benzene, methyl alcohol and 
ethyl alcohol were also of cp grade. Iron carbonyl] and 
tetraethyl lead were obtained from the Ethyl Corp. 


Experimental Results 


Ignition Temperatures 


Ignition temperatures were measured for the 10-mm jet 
under two conditions: (a) 100 cc/sec hot air going into pure 
fuel. This leads to an ignition zone 7 to 10 cm above the 
furnace exit, and a conventional diffusion flame strikes back 
and settles on the furnace exit; (b) 100 cc/sec hot nitrogen 
going into a fuel-air mixture. A luminous column forms in 
the center of the jet for temperatures below ignition tempera- 
tures. The column grows in length and width as the tem- 
perature of the jet is increased. At the ignition temperature 
the column branches out, and a propagating flame moves 
throughout the explosive mixture. All ignition temperatures 
refer to the temperature in the center of the jet directly above 
the furnace exit. 

Tables 1 and 2 give the ignition temperatures for cases (a) 
and (b) together with older values (3) for comparison. The 
values of Table 2 refer in both cases to the ignition tempera- 
tures of the most favorable mixture strength, which is always 
close to stoichiometric (3). The values with the 10-mm jet 
are as much as 150 C less than the previous values with smaller 
jets. 


Effect of Additives on Hot Gas Ignition Temperature 


The 4-mm-diameter furnace was used in these experiments; 
the igniting jet was either nitrogen or carbon dioxide of 
35 cc/sec total flow. 


Ethane-Air Mixtures 


Lead tetraethyl and iron carbonyl, when mixed with the 
outer explosive mixture and with carbon dioxide used as hot 
igniting jet, had no influence on the hot gas ignition tempera- 
ture within +5 C. In the case of iron carbonyl, a bright 
yellow luminous column formed before ignition, indicating the 
effective introduction of the additive. 

The results for chlorine and methyl bromide are shown in 
Fig. 1 for a hot nitrogen jet. Chlorine in the outer mixture 
raises hot gas ignition temperatures only slightly. For 


CH3 Br (outer) 

x 

1,200-- 
Cl2(jet) 

3 

1,000 

* ) 1 2 3 4 5 


ADDITIVE IN JET OR OUTER MIXTURE, PERCENT 


Fig. 1 Influence of additives on ignition temperature of ethane- 
air mixtures 
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Table 1 Ignition temperatures for hot air jet into fuel, C 
Jet diam, 10 mm Jet diam, 4 mm | 
Fuel 100 ce/sec air 35 cc/sec air(3) | 
methane 1040 1190 
ethane 840 945 
propane 885 1010 
n-butane 910 1025 
ethylene 765 875 
propylene 930 1060 
butylene 960 1070 
acetylene 655 755 
carbon monoxide 
(ep) 785 not tested 
CO (ep+ 5 per 
cent He) 660 not tested 
CO (commercial) not tested 765 
hydrogen 640 670 


Table 2 Ignition temperatures for a jet of 
hot nitrogen into fuel-air mixtures, C 
Fuel in Nz jet 10 mm Nz jet 7 mm 
fuel-air mixture 100 ce/sec 60 ce/sec (3) 
methane 1325 > 1200 
ethane 915 1015 | 
propane 990 1110 
n-butane 990 1095 
ethylene 900 985 
propylene 1060 1150 
iso-butylene 1070 1155 
CO (ep) 870 not tested 
CO (+ 10 per 
cent He) 720 not tested 
CO (commercial) not tested 795 
| hydrogen not tested <755 
benzene 1018 not tested 
| methyl! alcohol 946 not tested 
ethyl alcohol 964 not tested 
| diethy] ether 956 not tested 


comparison, Fig. 1 also shows the influence of chlorine on the 
ignition temperature when added to the hot jet. Itis apparent 
that this additive is more effective when present in the hot 
jet where the critical reactions leading to ignition occur. 
Methy! bromide has a strong inhibiting effect. Two per cent 
of this compound in the outer mixture raises the ignition 
temperature by nearly 150 C. Nitric oxide when added to 
the outer mixture in small concentrations lowers the ignition 
temperature slightly; in larger concentrations a small 
inhibiting effect becomes noticeable. 


Carbon Monoxide-Air Mixtures 


Both chlorine and methy!] bromide had a strong inhibiting 
effect on the ignition of carbon monoxide-air mixtures. In 
both cases a concentration of 4 per cent of the additive in 
the outer mixture raised the ignition temperature by about 
400 C. 


Limit Flame Temperatures 


The limiting oxygen indexes as measured with the hemispher- 
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Table 3 Limiting oxygen indexes for various fuels 
Limiting Limiting 
oxygen oxygen Limit 
index, index, temperature, 
Fuel this paper from (4) C 

methane 0.142 0.139 1518 
ethane 0.119 0.118 1363 
propane 0.129 0.127 1422 
n-butane 0.128 re 1434 
ethylene Te 0.105 1335 
propylene 0.127 ae 1489 
iso-butylene 0.130 Soil 1494 
CO (+10 per 

cent He) 0.065 1041 
CO (ep) eo 0.076 1177 
hydrogen 0.051 0.054 811 
benzene ee 0.133 1540 
methy] alcohol 0.111 1258 
ethyl] alcohol 0.126 1405 
diethyl ether 0.124 1386 


ical diffusion flame burner are given in Table 3. This table 
also contains values for relevant fuels from an earlier pub- 
heation (4). These values are in close agreement with those 
of the present investigation. Limit adiabatic flame tem- 
peratures are listed in the third column. 


Discussion 


A comparison between Tables 2 and 3 shows that a relation 
exists between hot gas ignition temperatures and limit 
temperatures of stoichiometric flames. In both cases re- 
actions occur under near-stoichiometric conditions in a large 
dilution with nitrogen. The implication of this relationship 
is that the same reactions are responsible for both flame 
propagation under limiting conditions and hot gas ignition 
under critical conditions. Fig. 2 shows a plot of hot gas 
ignition temperature against limit flame temperature. Al- 
though as a general rule it can be seen that high limit tem- 
perature always goes with high ignition temperature, excep- 
tions do exist. Methane clearly has an exceptionally high 
hot gas ignition temperature. Methyl alcohol and benzene 
also deviate from the proposed relationship. 

Accordingly, the correlation between hot gas ignition and 
limit temperature of stoichiometric flames must be considered 
a relatively loose one, and activation energies as calculated 
from the limit data are not directly applicable to hot gas 
ignition. Limit temperatures calculated from lean limit 
data show no improvement when plotted against hot gas 
ignition temperatures. The correlation is, in fact, much 
worse. 

Although the correlation between ignition temperature and 
limit temperature is only a loose one, it applies to observations 
made other than in direct ignition experiments. The influence 
of additives on ignition temperature can be understood from 
this correlation. Methyl bromide added to fuel-air mixtures 
narrows the region of flame propagation and increases the 
limit temperatures for all fuel-air mixtures (see for example 
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HOT- GAS IGNITION TEMPERATURE, °C. 


CO+10% H2 


800 900 1,000 1,100 1,200 1,300 1,400 = 1,500 
LIMIT TEMPERATURE, °C. 
Fig. 2 Correlation between hot gas ignition temperature and 
limit temperature of flames 


(5)). Even smaller details of flame inhibitions reflect them- 
selves in hot gas ignition. Chlorine is much less effective 
in inhibiting hydrocarbon flames than methyl bromide (see 
(6)) and the same tendency can be seen for the ignition 
temperature in Fig. 1. Carbon monoxide flames, on the other 
hand. are strongly inhibited by chloride (6), and again it was 
found that both methy] bromide and chlorine increase the igni- 
tion temperature of carbon monoxide drastically and by about 
the same amount. The effects of iron carbonyl and lead 
tetraethy! on limit temperature are unknown but are presum- 
ably negligible, as no influence on burning velocity could be 
found (7). Similarly, these two compounds have no in- 
fluence on the hot gas ignition temperature, whereas both 
compounds are known to influence spontaneous ignition 
temperatures greatly. 
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Technical Notes 


Hypersonic Viscous Shock Layer 
SINCLAIRE M. SCALA! 
General Electric Co., Philadelphia, Pa. 


DECADE ago Tsien (1)? (as well as others) and, more re- 

cently, Adams and Probstein (2) have attempted to de- 
fine the different regimes of gaseous interactions during high 
altitude flight. In this note some results are presented which 
are pertinent to the flight of hypersonic lifting vehicles com- 
posed of axially symmetric and two-dimensional elements; 
see Fig. 1. 


Analysis 


It may readily be shown that in the continuum regime the 
nominal boundary layer thickness at a stagnation point is 
directly proportional to the square root of the nose radius, 
i.e. 


8/V Rp = f (Va, Alt, Tw) (1) 


while the shock layer thickness or shock detachment distance 
at a stagnation point is directly proportional to the nose 
radius, i.e. 


6:/Re = f (Va, Alt) [2] 


It follows that a suitable parameter useful in defining the flow 
regime where incipient merging of the shock layer and 
boundary layer occurs is 


= f (Vo, Alt, Tw) [3] 


Upon restricting the flight Mach number and surface tem- 
perature at high altitude to the limits 


16<M < 24 
1000 < Ty < 3000 R 


and employing values of the nominal boundary layer thick- 
ness at the stagnation point (3,4) based on recent exact 
numerical solutions of the hypersonic boundary layer equa- 
tions and the results of Li and Geiger (5) for the shock de- 
tachment distance, one obtains the results shown in Fig. 2. 
That is, for the regime of interest defined by Equation [4], all 
the data fall on curves, all of which lie within the bands de- 
picted in Fig. 2. It is seen that at a given altitude the value 
of the merged layer parameter is larger for axially symmetric 
bodies than for two-dimensional bodies. This follows be- 
cause, although the boundary layer thickness 62p is of the order 
of +/2 times larger than 3p, the shock layer thickness 6,., is 
considerably larger than 6;,, for a given value of the nose ra- 
dius, consequently 


6s 2D bs 3D 
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TWO- DIMENSIONAL FLOW 


AXIALLY - SYMMETRIC FLOW 


Fig. 1 Schematic of merged hypersonic viscous layer for two- 
dimensional and axially symmetric flow 


For convenience, the mean free path was calculated up- 
stream of the shock, at the outer edge of the boundary layer, 
and at the wall. By definition 
1 
rnNo? 


and hence, if one takes ¢ = 3.6A (constant), then in fps units 
there follows 


[6 


9.196 x 10-7 [7] 

p(Ib/{t®) 
The initial decrease in the mean free path, as a free stream 
particle crosses the shock, is due to the shock compression, 
whereas, in the boundary layer, which is assumed to be at 
constant pressure, the mean free path drops directly with the 
temperature (for constant a). Thus. in Fig. 3 are shown Ao, 


Epitor’s Note: The Technical Notes and Technical Comments sections of ARS JouRNAL are open to short manuscripts describing 
new developments or offering comments on papers previously published. Such manuscripts are published without editorial review, 
usually within two months of the date of receipt. Requirements as to style are the same as for regular contributions (see masthead page). 
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Fig. 2 Merged layer parameter at the hypersonic stagnation 
point 


\. as functions of Mach number, and A,» as a function of 
Mach number for a surface temperature 7’, = 1000 R. 

In these calculations, the free stream conditions were based 
on a model of the Earth’s atmosphere proposed by ARDC (6). 
Normal shock tables were averaged based on work by Logan 
(7), Feldman (8) and Huber (9). 


Discussion 

The merged layer parameter shown in Fig. 2 immediately in- 
dicates the limits of applicability of the standard aerodynamic 
technique of separating the flow field into a shock wave, in- 
fluscid flow and boundary layer flow for the high altitude hy- 
personic flight regime. 

If preferred, one can also combine the data of Figs. 2 and 3 
to form a critical Knudsen number. Thus, if one sets the 
ratio 6/4, equal to unity, the data of Fig. 2 may be interpreted 
as a relationship between Rp and the altitude at which 
merging occurs. Then a critical Knudsen number may be 
defined which is based on the free stream mean free path \.. or 
on either A, or Ap. 

For example, if one selects 


Kno = X\o/Rp [8] 
then the data taken from Figs. 2 and 3 yield the limits 
6 X 10-? < < 10- 
2X 10-? < < 3 X 107 


Similarly, based on A, one obtains 
2X < Kttesn < 5 X 107? 
7X < $1.7 X 10-7, 


1959 


>|— 


M=16 
195 20 
al 24 
(Ty, =1000° R) 


lo” 


io8 | | | | 
0 100 200 300 400 


ALTITUDE « 107? FT. 


Fig. 3 Mean free path, upstream, local and at the surface of the 
body 


Other definitions can also be used, of course. 

In general, a small nose radius promotes the formation of a 
viscous shock layer. However, in view of inequality [5], it is 
possible for different portions of composite bodies to be opera- 
tive simultaneously in different flow regimes with catastrophic 
consequences. Fortunately, the trend indicates that with 
proper care two-dimensional lifting surfaces can be thinner 
than the center body and still operate in the same flight re- 
gime. 

It is understood, of course, that when 6/6; = 0(1), ordinary 
continuum concepts must be replaced by a new continuum 
model, such as perhaps the one recently considered by Prob- 
stein and Kemp (10), and independently by Hoshizaki (11). 
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Nomenclature 

Alt = altitude 

Kn = 2/Rz, Knudsen number, based on nose radius 
M = Mach number 


moles per unit volume 
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N = Avogadro number 

p = pressure 

Res = nose radius of body 

T = temperature 

Vo = flight speed 

6, = shock layer thickness 

6 = boundary layer thickness 
o = collision diameter 

= mean free path 
Subscripts 

2D = two-dimensional flow 
3D = three-dimensional flow 
e = edgeof boundary layer 
w = wall 

co = free stream, upstream of shock 
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Hypervelocity Vehicles at Large Angles 
of Inclination 


W. H. T. LOH! 


Chance Vought Aircraft, Inc., Dallas, Texas 


ESULTS on hypervelocity vehicles, either from boost 
glide or orbital re-entry, at small angles of inclination 
and constant lift-drag ratio were obtained recently by Allen 
and Eggers?; they are: 
1 Flight range 


ro\(L gro Ve 
2 Total aerodynamic heat input to the vehicle from the be- 


ginning to the end of the flight 


Q Cr'S 
(1/2)mV CpA 


[2] 
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tan 6 In 


3 The time rate and maximum time rate of average heat 
input per unit area 


2 
1 Cr mq [3] 


dt 2(L/D)CpA 


gro 
1 Cr'mg 
4 
( dt 3/3 (L/D)CpA Vore (4! 
J max 3 


4 The time rate and maximum time rate of local stagna- 
tion region heat input per unit area 


2 1/2 
6.8 xX 10- (: V2 [6 
dt CpA(L/D)o gro 
10-6 (gro) (7] 
CpA(L/D)o 


9 
Vidi. = [8] 


At large angles of inclination, the differential equation of 
motion is 


dV? [dé 2 29 ( L.. ) 
— 24 -- 6)=0 


When the angle of inclination is large but constant, closed 
form approximate solutions may be obtained as follows’: 
1 Flight range 


™ — (L/D) tan 6) — V2 
2 Total aerodynamic heat input to the vehicle from the 
beginning to the end of the flight 


— (L/D) tan 6] — _ 4) 
— (L/D) \v,) ** [10] 


3 The time rate and maximum time rate of average heat 
input per unit area 


2 
dHw _ 1 r'mg y( 
dt 2 (L/D)CpA gro 
| 
for V, > V(1/3) gre 
1 Cr’ mg cos 6 
12 
( dt 3/3 (L/D)CvA Vare 
V /dHay [13] 
dt V3 


4 The time rate and maximum time rate of local stagna- 


’ Other assumptions are the same as those in the Allen-Eggers 
paper. However, theoretically speaking, constant angle of in- 
clination can only be maintained by continuous programming of 
lift-drag ratio with velocity according to 


(4) (2c0t — (V2/gro)) 
D} \ Bro J (V*/gro)in(V/Vs)? 


But for regions of our interest, it may be treated approximately 
here as a constant. This is because, numerically, the maximum 
variation of (L/D) is only about 20 per cent for (V2/gro) from 0.8 
to 0.3 at (Vs?/gro) = 0.9; and it reduces to about 15 per cent for 
(V2/gro) from 0.6 to 0.2 at (Vs?/gro) = 0.7. Therefore, within the 
said regions of interest, (L/D) is approximately treated as con- 
stant in this derivation. Note also it is a function of 6. 

‘When V; is smaller than the quantity, the maximum occurs 
at the start of the unpowered flight. 
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tion region heat input per unit area 
2\ 1/2 
_ 68 x 10 (: v2 [14] 
dt (L/D)ACpo gro 
for Vy > W(2/3)gr0 


dH, 2 cos 0 


"(> ore [16] 
max 


5 Total time of flight 


on (: tan 0) 


1 L =1 
~ 2 cos 0 (: D 0) 


\Vhen the angle of inclination is small 


tan 


quations [9-16] reduce, respectively, to Equations [1-8]. 


Nomenclature 

= reference area for lift and drag expressions 
Cp = drag coefficient 

= equivalent skin friction coefficient = 


1 V 


skin friction coefficient 
drag, lb 
acceleration due to force of gravity, fps? 


x 


convective heat transferred per unit area 

L lift, Ib 

m mass of vehicle 

r distance of glider measured from center of Earth, ft. 
When altitude of the glider is small in comparison with 
Earth radius, r may be taken approximately as ro 

ro = radius of Earth, ft 

Q = convective heat transferred 

R = range, ft 

s = distance along flight path, ft 

S = surface area 

t = time, sec 

V = velocity, fps 

p = air density 


Vv 


EARTH'S 
SURFACE 


Fig. 1 Large but constant inclination glide trajectory 
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0 = angle of inclination or angle of flight path to local Earth 
horizontal, deg 
o = nose or leading edge radius of body or wing 
B = constant in density-altitude relation: p = poe~®¥. Here 
po = reference density, and y = altitude 
Subscripts 
e = condition at end of unpowered glide 
Fi = condition at end of powered boost or condition at be- 
ginning of unpowered glide 
= local conditions 
= stagnation conditions 
av = average values 
max = maximum values 


gro D tan 6 Vy 


[17] 


Stagnation Point Heat Transfer in a 
Subsonic Jet of Arc-Heated Air’ 


R. R. JOHN? and W. L. BADE’ 


Aveo Research and Advanced Development Div., 
Wilmington, Mass. 


Heat transfer from a free subsonic jet of arc-heated air 
to the stagnation region of a plane surface has been deter- 
mined experimentally. The enthalpy, pressure and ve- 
locity gradient of the flow at the edge of the boundary 
layer were estimated on the basis of mass flow and energy 
balance measurements. Heat flux values predicted by a 
correlation formula for stagnation point heat transfer in 
dissociated air agree approximately with the experimental 
values. 


Description of Experiments 


HE RECENT development of are jet and are tunnel 

facilities (1)4 has provided a means for continuously heat- 
ing a gas stream to the enthalpy and temperature levels 
associated with hypersonic flight. The present note reports 
preliminary measurements of heat transfer from a free sub- 
sonic jet of arc-heated dissociated air to the stagnation region 
of a plane surface. 

The arc jet unit used in the experiments is shown sche- 
matically in Fig. 1. The gas to be heated is injected tan- 
gentially into the chamber. The are discharge is ignited by 
exploding a tungsten bridge wire which initially connects the 
carbon cathode and water-cooled copper anode. The in- 
jected gas flows through the region occupied by the are dis- 
charge, momentarily becomes part of the discharge, and then 
leaves the system through a plenum chamber and a subsonic 
exit nozzle. The entire arc chamber, including the anode, 
plenum chamber, exit nozzle and cathode holder, is water- 
cooled. 

For the measurements reported below, the arc jet unit was 
connected in series with a direct current diesel generator 
(rated output 1500 amp at 250 v) and a 0.12-ohm ballast re- 
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sistor. The internal diameter of the exit nozzle was 0.50 in. 
The jet was exhausted into the open laboratory. 

The steady-state heat flux from the plasma jet to the stag- 
e nation region of a plane surface was determined using the 
Sie constant-flow water calorimeter shown schematically in Fig. 
ia 2. Measurements were made of the rate of heat transfer 
through the central circular area, which was aligned co- 
axially with the are jet nozzle. The water-cooled annular 
guard ring surrounding the calorimeter had the dual purpose 
of preventing heat transfer to the sides of the calorimeter and 
of establishing an “infinite plate” geometry, for which ve- 
locity distribution information was available in the liter- 
ature. Stagnation region heat transfer from the plasma jet 
was measured as a function of air mass flow; the results are 
shown in Fig. 3. It is noted that the observed heat flux is 
almost independent of air mass flow. 
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Fig. 1 Ejiectric arc heater 
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Fig. 2. Constant-flow water calorimeter 
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Estimation of Jet Properties 


In order to compare these results with existing information 
about heat transfer in dissociated air, it is necessary to con- 
sider the pertinent properties of the plasma jet. These are, 
as indicated in Equation [4a], the stagnation point veloc- 
ity gradient and the thermodynamic state ofthe air at 
the outer edge of the boundary layer. These quantities were 
estimated for the plasma jet in the manner described below. 

The mean stagnation enthalpy of the air in the jet was 
calculated from the relation 


h, = Pair/M = nP/M = 0.948 X 10-3 nVI/M 


The power absorbed by the air was estimated from an ar 
energy balance, in which the power to the arc was deter- 
mined from measurements of are voltage and current, and the 
power loss to the cooling water from the weight flow and 
temperature rise of the water. The air mass flow was 
measured using a flowmeter. 

The mean velocity of the jet was estimated from the 
measured air mass flow and the mean stagnation enthalpy, 
assuming the subsonic exit flow to be quasi one-dimensiona! 
and in thermal equilibrium, by simultaneous iterative solution 
of the equations of continuity and energy 


M = (2/4) D?p,U [2a] 

h, = ho + U?/2SQ [2b] 

and the equation of state (2) for equilibrium dissociated air 
po = f(ho, Po) [2c] 


The stagnation pressure was estimated from the jet velocity 
by assuming an isentropic compression (2) of the air flowing 
from the nozzle exit to the stagnation point. The stagnation 
pressure and the stagnation enthalpy [1] determine the ther- 
modynamic state of the gas at the edge of the boundary layer. 

Finally, on the basis of information available in the liter- 
ature (3,4) for the incompressible flow of a free jet against a 
large plane surface, the stagnation point velocity gradient was 
approximated by the expression 


B ~U/D [3] 


The variation of estimated jet properties with air mass flow 
is shown in Table 1. 


Table 1 Variation of estimated jet properties with air 
mass flow 

M hs Ps B 
lb/see Btu/lb atm sec™! R 
0.004 8800 1.02 31,000 11,800 
0.008 6500 1.08 49,000 10,800 
0.012 5100 1.18 62,000 9800 
0.016 4200 1.29 72,000 8600 
0.020 3700 1.42 81,000 7700 
0.024 3400 1.56 89, 000 7400 


Comparison With Theory 


Fay and Riddell (5) have established from numerical 
boundary layer solutions that stagnation point heat transfer 
can be represented by the correlation formula 


q = 0.763 X 
[1 + — 1) (hp/h,)] (hy — hy) 


This formula agrees with heat transfer rates measured in 
shock tubes (6) if the following nominal transport property 
values are used 


[4a] 


Np, = 0.71 Ny = 1.4 [4b] 
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where 
Mo = 1.16 X 10~* lb/ft-sec 
= 492R 
S = 204R 


Equations [4] have been applied to the plasma jet situation, 
using the estimated jet properties shown in Table 1. The 
curve of predicted heat flux is shown superimposed on the 
experimental] results in Fig. 3. The correlation formula [4] 
is in reasonably good agreement with the measured values, 
and suecessfully accounts for the empirical observation that 
heat flux is nearly independent of air flow. 

Work is currently under way on the measurement of heat 
flux at other arc power levels, jet diameters and exit pressures. 
These results will be published shortly. 


\cknowledgment 


The authors wish to acknowledge many helpful discussions 
with Dr. Malin, Dr. Polestra and Dr. Barriault of Avco, 
RAD. R. N. Schweiger designed the original Aveo, RAD 
clectric are air heater, using a water-cooled copper anode first 
designed and used by T. Brogan of Avco Research Labora- 
tory. 


Nomenclature 


nozzle diameter, ft 

gravitational constant, 32.2 ft/sec? 

enthalpy, Btu/lb 

dissociation enthalpy per unit mass o: air, Btu/Ib 


Ap = 

/ = are current, amp 

= Joule’s constant, 778 ft-lb/Btu 

M = air mass flow, lb/sec 

Vie = Lewis number for atom-molecule binary diffusion in air 

Vp, = Prandtl number for air 

p = pressure, lb/ft? 

P = power to arc, Btu/sec 

/ = stagnation point heat flux, Btu/ft? sec 

T = absolute temperature, R 

= jet velocity, fps 

V = arc voltage, v 

8 = stagnation point velocity gradient, sec™ 

n = are efficiency = 1 — (power loss to cooling water)+ 
(power to arc) 

un = viscosity of air, lb/ft-sec 

p = density of air, lb/ft® 

Subscripts 

0 = static conditions at nozzle exit 

s = stagnation conditions 

w = wall conditions 
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Long Period Oscillations During 
Atmospheric Entry 
GEORGE S. CAMPBELL! 

Hughes Aircraft Co., Culver City, Calif. 


An approximate theory is presented for the oscillatory 
characteristics of the long period motion during atmos- 
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pheric entry at constant C,; and Cp. Comparison with 
numerical results for the complete equations of motion 
indicates that the approximate theory provides an accurate 
description of the time-varying amplitude decay and fre- 
quency during entry. 


HE PROBLEM of atmospheric entry has been discussed 

in a number of papers, such as Chapman’s recent study of 
the problem (1).2 Numerical solutions of the equations of 
motion have shown that the use of lift during entry may cause 
a skipping motion. This motion has a long period but is 
quite different from the classical phugoid motion. 

It is the purpose of this note to present simple approxima- 
tions describing the oscillatory characteristics of the long 
period motion. The analysis presented is intended to pro- 
vide an understanding of the oscillatory behavior rather than 
a detailed computation of the trajectory path. It is hoped 
that the results obtained will be useful in the design of control 
systems for re-entry vehicles. 


Average Trajectory 

Before oscillatory characteristics can be described, it is 
necessary to have a simple description of the “‘average’’ tra- 
jectory during descent. Siinger’s gliding trajectory appears 
adequate for this purpose, but for the sake of clarity the 
derivation is repeated. 

The full equations of motion for descent at constant angle 


of attack are 


mV = —CpqgS + W sin y 
—mVy — (mV?/r) cos y = — W cos 


where y, the flight-path angle, is positive downward and is 
related to the altitude y 


sin-y = —9/V 
An exponential density-altitude relation is used (1) 
p = 
and the dynamic pressure is 
q = pV?/2 


The basic glide trajectory is calculated by omitting the 
terms (mVy) and (W sin y) from Equations [1] and by 
assuming smal] flight-path angles. Under these approxima- 
tions the above equations are satisfied by 


Vive = Vor tanh T (2] 
B vYave = log (Xo sinh? T) [3] 
= = J esch? T [4] 
Po No 
dave = (W/C1S) sech? T [5] 
in which 
Co fg 
T C, t) 
Xo 
2m 


The constant of integration ¢; must be determined so that the 
proper solution of the singular perturbation problem is chosen. 
For a trajectory such as that shown in Fig. 1, this parameter 
was determined by matching Equation [3] with the machine- 
calculated time for 30 miles altitude. 


2 Numbers in parentheses indicate References at end of paper. 
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Fig. 1 Altitude and velocity history for entry from 60 miles; 
yo = 0; Vo = 25753 fps; W/S = 18 1b/ft?; Cx = 0.54; Co = 
0.52 
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Fig. 2 Amplitude of long period oscillations during entry from 
60 miles; yo = 0; Vo = 25753 fps; W/S = 18 lb/ft?; Ci = 
0.54; Cp = 0.52 
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Fig. 3 Frequency of long period oscillations during entry from 
60 miles; yo = 0; Vo = 25753 fps; W/S = 18 lb/ft?; Ci = 
0.54; Cp = 0.52 
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Oscillatory Motion 


Superimposed on the average trajectory are altitude and 
velocity oscillations of the type shown by Chapman and in 
Fig. 1. In order to describe these oscillations, Equations [1 | 
are linearized about the average trajectory 


h< Yavg 
w < Vane 


= Yar th 
V = Vag t+ w 


When these substitutions are made in Equations [1], the rela- 
tions for w and hare obtained 


w qd y\w h qd 
( W/CpS y+ 
h Vh q yw q 
9 #9 yt ( W/CS* gr 
(6) 


The subscript avg has been omitted from q, V and y for pur- 
poses of clarity. The exponential term e~®" was replaced by 
(1 — 6A), so that the analysis, to be strictly valid, only ap- 
plies when Bh< 1. 

In computing frequency of the oscillations, only two terms 
of the lift equation are required 


Bh = 0 7a] 
L 
Substitution of average dynamic pressure, Equation [5], gives 
the equation 
d*h 
ar: + (K? sech* T)hkh=0 [7b] 
in which 
K = VBr (C1/Co) 
Since 1/8 r = 30 for the Earth’s atmosphere (1), K is large 
for (C./Cp) of order unity. An approximate solution to 
Equation [7] may therefore be obtained by the W B K method 
(2) 


Bh = SsechTaT 


The frequency of the altitude oscillations, obtained from h/h, 
is 


w = WBgsech T [Sa] 


or 


Q. q 
8b 
w= VBg W/C,S [8b] 

The second relation serves as a basis for simplifying Equa- 
tions [6], since w/w and h/h are of order w. If it is assumed 
that qave/(W/CdS), Vave/~/gr and (Cz/Cp) are of order 
unity, Equations [6] become, for small avg 


CpS 
+ + gph = 2 4 [9a 
w [9b] 
m 


Equation [7a] is substituted into Equation [9b], and the re- 
sulting relation is integrated to give 


w= 
The altitude oscillations are therefore given by 


Cp CS 
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d*h dh 
— 2d T + (K sech T)?h = 0 


AX = 3esch 27 + tanh T 


The damping term may be eliminated by the transformation 


h = HeS 


and the resulting equation solved approximately by the 
\'V B K method (2) to give 


B h = (sinh T)*/2 ettK Sf sechTdT [10] 


Equation [10] is the basic result describing the character 
ff the long period oscillations of a re-entry vehicle. The 
implitude envelope depends on density according to 


Bh ~ (Po/p)*/* 


so that the oscillations die out at the lower altitudes. Fre- 
juency of the oscillations is given by Equation [8a]. 


Numerical Results 


In order to check the validity of Equation [10], results 
from this equation have been compared with numerical solu- 
tions of the full equations of motion, Equations [1]. One 
of the trajectories analyzed is shown in Fig. 1. The average 
trajectory is described quite accurately by Equations [2 and 
3], witht, = 2965 sec. 

The amplitude decay and frequency of the altitude oscilla- 
tions of Fig. 1 are presented in Figs. 2 and 3, respectively. 
The approximate theory is compared with IBM results, 
and good agreement is shown. Similar agreement has been 
obtained for other lift-drag ratios, initial altitudes and flight- 
path angles. It is noted that the approximate equations, 
which were derived for small amplitudes (8 h< 1), give rea- 
sonable results even for relatively large-amplitude oscillations. 
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Viscous Interaction Experiments at Low 
Reynolds Numbers’ 


S. A. SCHAAF,? F. C. HURLBUT,* L. TALBOT‘ and 
J. AROESTY® 


University of California, Berkeley, Calif. 


HE DISTURBANCE in the surface pressure due to the 
displacement effect of a laminar boundary layer has been 
the subject of a considerable number of investigations during 
the past few years. The reader is referred to Bertram (1),® 
and the literature quoted therein, for a recent summary of 
these investigations. 
The effect is most pronounced, and has been most exten- 
sively investigated, for the case of a very sharp, leading edge 
flat plate at zero angle of attack in supersonic flow. For this 
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HYPERSONIC PARAMETER, X=MaVC/Re, 
Fig. 1 Induced pressure on flat plate (sharp leading edge) 


case, adiabatic laminar boundary layer theory (2) indicates 
that the surface pressure p should be roughly proportional 
to the “hypersonic interaction parameter” 


X = M*V/C/Re, 


where 
M = free stream Mach number 
Re, = Reynolds number (based on distance from the leading 
edge) 
Chapman-Rubesin constant 


According to this result, the pressure becomes infinite as the 
leading edge is approached. Of course, this is not expected 
actually to take place, since the boundary layer assumptions 
are not valid in this region. It is presumed that some ap- 
propriate modification in the analysis could be carried out, 
and that—providing the gas does not slip along the surface— 
such an analysis would result in a predicted value of the sur- 
face pressure at the leading edge of the order of that existing 
behind a normal shock. If, on the other hand, the gas did 
slip along the surface, according to the well-known slip ve- 
locity boundary condition of kinetic theory, a substantially 
lower value for the predicted surface pressure would be ex- 
pected. 

The indicated analyses are quite formidable, and neither 
has yet been successfully carried out. In the present note, 
data are presented which seem to be more in accordance with 
the latter expectation, i.e., a surface pressure near the leading 
edge substantially below that existing behind a normal shock. 

The new data are presented in Fig. 1, along with previous 
data by Bertram (1) and Kendall (3). The new data were ob- 
tained in the University of California low density wind 
tunnel. Two sets of data were obtained. In one of these, the 
Reynolds number Re, based on the leading edge thickness of 
0.01 in. was in the range 10 < Re, < 85; in the second, the 
leading edge thickness was less than 0.001 in., so that Re; < 8. 
No effect of leading edge was detectable over these ranges, so 
that the results were presumably equivalent to the ideal 
“sharp” leading edge case. The flat plates were at recovery 
temperature; i.e., the data are for the adiabatic wall case. 
Two Mach number ranges were used, one at M ~ 3.7 to 4.1 
and the other at M ~ 5.5 to 5.8. The higher of these is quite 
close to that investigated by Kendall (3), and the agreement 
between his data and the present experimental results is satis- 
factory in the region of overlap. The present data depart 
quite abruptly from the general trend exhibited by both the 
theoretical analyses and the data of Bertram at slightly higher 
values of X. For the M ~ 3.7 to 4.1 case, the departure is at 
about X = 5. For the M ~ 5.5 to 5.8 case, the departure is at 


about X = 7.5. 
Since the disturbance in the pressure is due to the displace- 
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ment effect of the boundary layer, it seems reasonable to sup- 
pose that any slip flow effect would be associated with the 
parameter \../5,, where A» is the mean free path at the wall, 
and 6, is the boundary layer displacement thickness. A 
simple calculation, for air at reasonably high Mach numbers 
and for the adiabatic boundary layer, leads to the result that 
Aw/5s = 0.27X/M. Hence both sets of data depart from the 
no-slip trends for a value of \»/54 = 0.35, a not unreasonable 
value. Conversely, the Bertram data, if this criterion should 
prove applicable, might be expected to depart for a value of 
X = 12. 

The limiting values of the surface pressure as the leading 
edge is approached appear to be about 3.5 po for the case M ~ 
3.7 to 4.1, and 5.5 po for the case M ~ 5.5 to 5.8. These are 
substantially below the values of 20 po» and 42 po at M = 3.9 
and 5.7, respectively, which would obtain at the stagnation 
point behind a normal shock. Thus, both the manner in 
which the data depart from the boundary layer trend and the 
limiting value approached at the leading edge seem to suggest 
the existence of slip flow phenomena. 
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Theoretical Study of Recombination 
Kinetics of Atomic Oxygen’ 


S. T. DEMETRIADES? and M. FARBER 
Aerojet-General Corp., Azusa, Calif. 


Estimates are presented of the rate constants for atomic 
oxygen recombination due to three-body collisions and 
other homogeneous and heterogeneous processes. 


HIS note presents a summary of a theoretical study 

(1)4 of the reaction kinetics of atomic oxygen. Estimates 
are presented of the mechanism and magnitude of the recom- 
bination rates for atomic oxygen in the homogeneous phase 
as well as of the recombination rates in a heterogeneous system 
(i.e., on a catalytic wall). A rate constant for the oxygen 
recombination mechanism O + O + M — O, + M of the 
order of 10" to 10" cc?/mole*-sec, depending on the tempera- 
ture, has been derived theoretically. These estimates are 
based on the results of collision rate theory, absolute reac- 
tion rate theory and the interpretation of experimental data 
available in the literature on other three-body recombination 
reactions and surface recombination reactions. On the basis 
of these estimates, Demetriades (2) has made a theoretical 
study of an upper atmosphere atomic oxygen powerplant. 


Homogeneous Recombination Mechanisms 


There are two important homogeneous mechanisms which 
contribute to the production of oxygen molecules in the 
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presence of oxygen atoms, ozone and inert bodies (A, N,, 
etc.). These are the triple-body collision mechanism 


k 
0+0+M=0.4+M (1) 


and the well-known ozone mechanism 


ke 
0; + O + 20; (2 
and 
ks 
M+0,;=M+0,+0 [3 


The total rate of production of oxygen molecules due to 
the combined mechanisms, in moles/ce-see, is 


2 k.(O) + k3(O3) (M) [4 
where (QO), (Oz) and (O;) are the concentrations in moles/cc 
of atomic oxygen, molecular oxygen and ozone, respectively, 
and (M) is the concentration of third bodies (i.e., usually the 
total concentration in moles/ec). The equilibrium con- 
stants K,,and K, are given by 


(Oz) or -12 

Ka (078 10 [(25,700/7) 1.390] [5] 

a= (Oz) (O) = 109[— (5,620/T) + 5.016 —logio T) [6] 
(Os) 


where 7 isin deg K. The rate constant k; was recently given 
by Kaufman (3,4). The rate k, was studied by Garvin (5) 
and also by Benson and Axworthy (6). Their experimental 
results yield 


ke = 1.07 X 1013 (7°/500) 10—1950/7 cc/mole-sec [7] 


The last two terms of the right-hand side of Equation [4] 
represent the rate of production of oxygen due to the ozone 
mechanism. Thus 


d(Oz) (O) (Oz) 


If the concentration of ozone obeys the steady-state approxi- 
mation and 


k-3(M) >> 


Equation [5] reduces to 


d(O.) ks (Os)? (0;)2 
k,— — = 9 
If ke(Os) >> ks (M) (O2), Equation [5] reduces to 
dO 
(02) (M) [10] 
dt ozone 


On the other hand, if the concentration of atomic oxygen is 
adjusted so that it obeys the steady-state approximation, and 
k3(M) > k(O), Equation [5] becomes 


HO)! Bhs 
dt ozone Kes (Oy) 
If k.(O) >> k3(M), Equation [5] reduces to 
dO. 
(00 () [12] 
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Finally, if the reaction given by Equation [2] is rate-con- 
trolling 


| 


& 2ke(O) (Os) (13) 
dt ozone 
and if 1 — (O)/Ka(O3)] 0 
d(Oz) 
dt 


It is clear that 4, and K-.s suffice to determine the rate of 
production of oxygen molecules in the important cases 
de-cribed by Equations [6, 8, 10 and 11]. However, the 
vaue of k; is required to determine the total rate of pro- 
du:tion of O, The authors of (3 to 6) presented a rate 
constant given by 


ky = 3 X (300/7')!-5 
1000K < T < 4000K 


cc?/mole*-sec 


Ou the basis of the same work, the temperature limits can be 
extended as follows 


5.6 X 10°(7'/500) 10-0987) < ky < BX 
(300/7')!-5 


< 1000K 


T > 4000K ky = 3 X 10" (300/T)!.5 


[15a 

This rate constant was derived in two different ways. 
First, the order of magnitude of k; in the temperature range 
of interest was established by considering the number of 
three-body collisions involving two oxygen atoms and one 
oxygen molecule. The efficiency for these collisions was 
then assumed to be the same as that for the well-known reac- 
tion given by Equation [2], since the activated complex 
consists of four oxygen atoms and the end products of two 
oxygen molecules in both cases. That efficiency is of the 
order of 0.05, and thus / is of the order of 5 X 10" in the 
temperature range from 1000 to 2000K. Second, the absolute 
reaction rate theory of Eyring, Gershinowitz and Sun (7) for 
the three-body recombination process was modified in such a 
manner that by employing the experimental rate constants 
for well-known reactions, such as H + H + H.— H, + Hand 
1+1+M-—I,+ M, an estimate of the rate constant k, 
was derived. At 2000K, the “scaling’’ on the basis of hydro- 
gen recombination yielded k, = 1.1 X 10", and the “scaling”’ 
on the basis of iodine recombination yielded k; = 1.9 X 10". 
Finally, the temperature dependence of k; was assumed to 
be the same as for the iodine recombination rate constant. 
The choice of 300K as the base temperature in the quantity 
300/7 was arbitrary. At temperatures lower than about 
1300K, the approximations introduced in “scaling” on the 
basis of the Eyring theory are such that the value of k; given 
by Equation [15] should be slightly higher than the correct 
value. On the other hand, the collision theory approxima- 
tion with the efficiency of the reaction O; + O — 20, should 
provide a definite lower limit. 

This estimate is for O. as the third body. For Ne, A or 
other inert molecules acting as third bodies, k, is slightly 
smaller. The rate constant for recombination in high tem- 
perature air may be as much as 10 times as great as given by 
Hquation [15] because of the uncertainty in the relative effi- 
ciencies of nitrogen and oxygen atoms as third bodies. Fur- 
ther, oxides of nitrogen may be formed and act as chain 
carriers for oxygen recombination and possibly increase the 
rate, although no reasonable chain mechanism has been 
found by these authors. Thus, as already stated by one of 
these authors (2), the rate constant k, for atomic oxygen as a 
third body is 


ky = 3 cc?/mole*-sec {16} 
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with essentially the same temperature limits. 
Some values of k; as a function of temperature from Equa- 
tion [12] are presented in Table 1. 


Table 1 Present estimate of rate constant ‘, for oxygen | 


recombination by three-body collision mechanism (for | 
O, as a third body) 
k,, cc?/mole?-sec | 
300 <3 X 
1000 4.9 10% | 
1500 2.7 x 10% | 
2000 1.7 X10" 
2500 1.2 x 10% | 
3000 9.5 104 
3500 7.6 10" 
| 4000 6.7 X 10 | 
5000 4.5 10" | 
8000 2.2 X 10% 


The integration of Equation [4] has been carried out in 
closed form for the two limiting cases when the triple-body 
mechanism or the ozone mechanism predominates, and [4] 
has also been integrated graphically for various intermediate 
cases (1). At the higher temperatures and for equilibrium 
initial concentration of Os, only triple-body recombination is 
significant. Tables were also presented (1) of the equilibrium 
composition of the O + O2 + O 3 system in the range of tem- 
peratures from 1000 to 5000K and at 107%, 5 x 107° and 
10-* atm. 


Heterogeneous Recombination Mechanisms 


The concentration of atomic oxygen determines which 
heterogeneous mechanism is fastest. The first possible mech- 
anism is 


O + O + wall — O2 + wall [17] 


i.e., when the wall acts as third body. Then, since the wall 
is very efficient as a third body 


d(Oz) wen 
= 1.5 X 105O)?(T)! mole/sec-em? 
dt | wall+ 2(0) 
The second possible mechanism is 
O + (O adsorbed on wall) ~ O2 [19] 


This is not a 100 per cent efficient process. The rate of pro- 
duction is 


dt | O + wall—O 

With an efficiency f = 0.001 (for glass, see (8)) this becomes 

d(Oz) 
t 


= 9 10*%O)f(T)} mole/sec-cm? [20] 


= 0.9(O) (T)!? mole/see-cm? [21] 
The efficiency f approaches unity for the noble metals and 
metallic oxides (especially in the presence of traces of OH), 
and may be low as 10~ for “poisoned” glass surfaces or glass 
surfaces treated with phosphoric acid. Both heterogeneous 
mechanisms occur simultaneously, but that given by Equa- 
tion [19] is faster when (QO) is small. 


Combined Production Rate 


For most practical purposes, the term (O2)/Ka(O)? in 
Equation [4] is small compared to 1, and, at least at high 
temperatures, the term due to the ozone mechanism is small 
compared to the triple-body mechanism. Then, at low pres- 
sures, the combined rate of production of oxygen becomes 
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| 

4 
= 
= 
' 


approximately 


To) 


dl > k,(O)*{M) + 9 X 10*%O) mole/sec-ce [22] 


where A is the catalytically active surface area in cm? of a 
control volume V in cc. When the control volume in cc is 
large compared to its surface area in cm?, wall recombination 
may be neglected. 


Conclusions 


With the rates presented in this paper, which are not so high 
as previously thought (9), either flow or static experiments 
can be designed to measure the specific rate constant /,. 
In these experiments, the visualization of the recombination 
reaction can be enhanced by injection of traces of substances 
whose reaction rates are known (e.g., NO, NOs, ete.) or in- 
clusion of N, to produce NO and the yellow-green afterglow 
caused by the reaction of oxygen atoms with nitric oxide. 
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Estimation of Performance of Fuels 


With Hydrogen Peroxide Oxidizer 


JOHN S. GORDON! 
Thiokol Chemical Corp., Denville, N. J. 


An empirical short method is described for determining 
rocket performance of fuels containing carbon, hydrogen, 
nitrogen and oxygen, with 90 per cent hydrogen peroxide 
oxidizer. The method is similar in approach to existing 
short methods for other oxidizers. The results in I;, and 
C* are uncertain by 0.3 per cent. 


HORT-CUT methods for estimating the performance 

parameters of various fuels containing carbon, hydro- 
gen, nitrogen and oxygen with nitric acid, oxygen and am- 
monium nitrate oxidants were developed by Johnston.? 
Data for a variety of fuels with 90 per cent H,O, have become 
available in the course of propellant investigations at Reaction 
Motors, Inc. These data are summarized in Table 1. Data 
from other sources are also included, but are older and be- 
lieved less reliable. In general, the semi-empirical approach 
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and nomenclature used by Johnston are retained in this an- 
alysis. Only the chemical formula and heat of formation of 
the fuel are required to make the performance calculation. 
Considerable time can be saved, compared to rigorous meth- 
ods, and the short-method results can be used as a starting 
point for further rigorous calculations. Johnston describes 
an approximate method of developing « performance curve 
(at several O/F ratios) from a single peak J, value; this type 
of extrapolation was not believed useful enough to warrant 
inclusion in this study. NH,ClO, oxidizer has been studied 
by the author.® 


Derivation of Method 


Analysis of the data in Table 1 shows that peak J,» occurs 
at a mixture ratio corresponding to the undissociated exhaust 
products No, (x/3) and (2x/3) CO, from a fuel with: 
the general formula C;H,O,N». Ninety per cent has thc 
composition H,O,.s%3 with a molecular weight of 31.237 an: 
a standard heat of formation of —49.012 Keal/mole, froin 
gaseous O2 and Hy at 25 C. Development of an expression 
for heat release per unit mass of propellant follows. 


Reaction Products 
CrH,O:Nw + a > 


22 2a+y w 


3 CO + H.O + Ns 
Oxygen Balance 
4 2 
3 2 
0.82634 = — z 
3 2 
a = 1.6132 + 0.605y — 1.212 


Generalized Reaction Products 


CrH,O.Nw + (1.6132 + 0.605y — 1.21z) H2Oj.5263 > 
x 
CoO + 3 CO. + (1.6132 + 1.105y — 1.212) + 5 N: 
Combustion Energy Release 


Quay 


») 
(26.42) + (94.05) + (1.6132 + 1.105y — 1.212) X 


(57.8) — (1.6132 + 0.605y — 1.212) (49.01) + Qr 
where Q, is the heat of formation of fuel. 


17.6la + 31.352 + 93.232 + 63.87y — 
69.942 — 79.052 — 29.65y + 59.302 + Qr 


Qay = 63.142 + 34.22y — 10.642 + Qe 
or 
U = 63.14 (C) + 34.22 (H) — 10.64 (0) + Qr {1} 


Mass of Oxidizer Per Mole of Fuel 
31.237 (1.6132 + 0.605y — 1.212) 
18.90 [(8/3) (C) + (H) — 2(0)] 


Energy Per Gram of ai 


3] 
Q Keal/em) [ 
where M is the molecular weight of fuel. 
The reader will note the similarity of this algebraic develop- 
ment to that given by Johnston for HNO; oxidizer. Values 
of Q were calculated for the eighteen fuels. Rigorous J,, 


3 Gordon, J. S., Bulletin of the 14th Meeting, JANAF Solid 
Propellant Group, vol. 4, 1958, p. 221. 
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Calculated frozen performance with 90 per cent H.O, (600/14.7 pressure ratio) 


H/C iss Ce Men ka ke Source 
atom 243.6 5159 18.97 1.203 1.250 RMI 
ratio 
4 253.5 5372 20.20 1.197 1.222 RMI | 
0.75 251.1 5312 22.36 1.195 1.223 RMI 
0.75 254.3 5388 22.08 1.198 1.225 RMI 
0 254 22.88 RMI 
0.6 253.1 5349 22.34 1.197 1.224 RMI 
1.2 252.3 5329 21:91 1.192 1.220 RMI 
1.4 254.6 5384 21.40 1.195 1.224 RMI 
2.4 244.9 5176 22.15 1.188 1.217 RMI 
1.97 247.5 5251 21 .55 1.187 1.217 RMI 


Table 1 

O/F Or Optimum 
(stoich) Fuel Keal/mole O/F To 
3:20 NH; —17.17 3.2 ar 2360 
5.03 —6.36 5.097 2652 
8.17 CeHie.se —67.145 6.9 0.73 2586 
2.36 12.05 1.6¢ 2560 
3.54 CH;0H —57.02 2.85 ers 2555 
3.97 —93.13 3.5 2668 
CH;NO, —21.28 0.50 2890 
—59.74 6.0 2545 
N.H,H2O — 57.95 1.6 2261 
5.03 C2HsN:; +12.72 4.0 0.80 2738 
5.52 C,H;N +71.39 4.0 0.725 2956 
5.52 +88.95 4.0 0.725 3013 
3.98  CyNe +147.14 2.3 0.58 3129 
5.64 C;H3N 105.82 3.9 0.69 3000 
5.52 C;H.O 30.63 4.1 0.7 2920 
6.29 C;H:N 73.52 4.7 0.75 2915 
3.32 CsH»N-O; —91.8 2.9 0.87 2770 
8.06 Co.66His —57.22 6.9 0.86 2770 

“ Accuracy questioned. 


values were plotted against Q as shown in Fig. 1, and a good 
straight-line relationship resulted with an average deviation 
of +1.3 see from the rigorous J,, 


Isp = 179.5 + 54.6 Q 


The 4-sec deviation for CyN: indicated that a correction 
factor for H/C ratio might improve the results. A good hy- 
perbolic relation between correction factor (X) and (H/C) 
was found empirically for the RMI-caleulated data only 


X) = -1.95(—) + 2.1 
(X + 5 


An arbitrary maximum value of 6.0 for X is dictated by the 
results for CyNe. The second approximation for /,, is thus 


= 171.6 + 54.6 Q — 1.95 [4] 


Here the average deviation has been reduced only 11 per 
cent, to +1.16 see, but the internal agreement for the RMI- 
calculated data is excellent, +0.33 sec. Noting that most of 
the fuels from sources outside RMI have negative heats of 
formation, correction factors based on Q,» and on Q? were 
attempted. Plots of A (long J,, — second approximation 
[,,) as functions of Qp and of Q? showed only a fair correlation 
with Qpr in the form 


A = —0.4 + 0.009 Q» 

The third approximation for J,,, is thus 
Iep = 171.2 + 54.6Q — 1.95 (C/H) + 0.009 Qe 

The average deviation for all data points has been reduced 
another 11 per cent, to +1.03 sec, but the internal consist- 
ency of the RMI-calculated data is poorer, with an average 
deviation of +0.47 see. The added complexity with the 
0.009 Qr correction term evidently is not justified, and this 
term has no physical significance anyway. 
Chamber Temperature 


; is proportional to the total energy release from a uni 
‘y proportional to the total energy rel from nit 
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mass of propellants divided by the mean specific heat of the 
reaction products. The mean molar specific heats of the 
primary combustion products are taken as: 


ts) 
T 3000 deg 


Species 
CO 8.148 
CO; 12.992 
Na 8.087 
H.0 10.856 


250 
° 


240 
1.1 1.08 1.26 1.34 1.42 
Fig. 1 Theoretical frozen J,, vs. Q = U/(M + X). First 


approximation: /,, = 179.5 + 54.6 Q 
531 


| 
P| The mean specific heat of products from one mole of fuel, at : 
245 - 


the mixture ratio for peak J,, is derived as follows 
CzH,O.Nw + (1.6132 + 0.605y — 1.212) 


= co + CO: + (1.613 + 1.105 y — 1.212) + 


2 
8.143 (= ) + 12.922 (5) + 10.856 (1.613 + 1.105y — 
1.212) + (5) 8.087 


(5.429 + 4.307) (C) + 17.511 + 11.996 (H) — 
13.136 (O) + 4.044 (N) 


P = 9.74 (C) + 17.51 + 12 (H) — 13.14 (O) + 4.04 (N) 


Values of P were then calculated for each fuel. Rigorous 
values of 7. were plotted against the heat release/average 
heat capacity parameter, U/P. A straight-line fit was ob- 
tained 

T. = 199 (U/P) + 2065 (deg K) [6] 


with an average deviation of +76 K for all 18 points and 
+46 K for the RMI-calculated points. No improvement in 
accuracy could be obtained by attempting a correction factor 
based on C atom ratio, or the atom-ratio parameter 2-0/ 
(4C + H). A straight-line plot of T, vs. Q fits the data 
with an accuracy of +120 K and is not very useful. The 
poor accuracy of this method for estimating 7’, is disappointing 
in view of the +15 K agreement obtained by Johnston with 
HNO; systems. 


Molecular Weight 


The ten sets of RMI-calculated data in Table 1 lead to a 
linear correlation between VM and the estimated flame tem- 
perature 


M = 6.45 + 0.00529 7, (average deviation +0.43) [7] 
The available data do not justify a more complex relationship. 


The estimated flame temperature from Equation [6] corre- 
lates better than the 7’, from exact calculation. 


Chamber Specific Heat Ratio (Frozen Composition) 

K, is essentially constant at 1.195 over the range of RMI- 
calculated data (average deviation +0.004). C,/C, at T, 
(for 600/14.7 expansion) is about 0.029 greater than at T,. 
Stoichiometric Mixture Ratio 

By inspection of Equation [2], this is seen to be 

— 18.9 [4(C) + (H) — 2(0)] 
M 
Mixture Ratio for Peak /,, 
This correlates roughly with the H/C atom ratio in the 


[8] 


fuel 


2 = 0.67 + 0.04 (2) < 0.92 9] 
Cc 


Characteristic Velocity 
Within the accuracy of the rigorously calculated data n 


Table 1, C* may be expressed as a linear function of the he.t 
release per unit mass 


C* = 3952 + 1.046Q (average deviation +16 fps)  [1')] 


An expression in the form C* = a + bQ — cQ* could have bec n 
developed from the limited data available, but this was not 
believed justified, since the best equation for J,» has no terin 
in Q*. The long-method C* values are based on k, rather 
than the effective y for frozen expansion (7./T, 
or the average k between chamber an 1 
exhaust k = 1/2 (ke + ke). 


Extrapolation to 100 Per Cent Hydrogen Peroxide Oxidizer 


Only four fuels were found in the literature with whic' 
both 90 and 100 per cent H,O, had been given theoretic:! 
treatment. In general, it seems reasonable to add 1 se 
I,, and 30 K flame temperature for each additional 1 per 
cent concentration, Optimum and stoichiometric mixtur:: 
ratios are 10 per cent less than for 90 per cent peroxide. 


Summary 

= 63.14 (C) + 34.22 (H) — 10.64 (O) + Qr 
X = 18.9 [(8/3) (C) + H — 2(0)] 

Q =U/(M +X) 


I.p = 171.6 + 54.6 Q — 1.95 (C/H) 
(frozen, 600/14.7 pressure ratio) 


9.74 (C) + 17.51 + 12(H) — 13.14 (O) + 4.04 (N) 
T. = 199 U/P + 2065 (deg K) 


| 


ll 


6.45 + 0.00529 T, 


_ 18.9(4-C + H — 2-0) [0.67 + 0.04 (H/C)) 
M 


= 


C* = 3952 + 1.046 Q 
“he expected uncertainties are 
1 sec 
T. + 70K 
M +04 
rp + 0.2 


C* + 16 fps 
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DIMAZINE 


unsym-Dimethylhydrazine, UOMH 


THE STORABLE FUEL 


Performance data on DIMAZINE-powered 
rocket engines in military missiles and satel- 
lite launching vehicles shows DIMAZINE 
to be a highly reliable rocket fuel. 


DIMAZINE provides fast, dependable hy- 
pergolic starts followed by smooth, stable 
combustion and easier shutdowns. Depend- 
able instant readiness is assured for years by 
its outstanding stability during storage in 
missiles. It also has high performance, high 


Putting 


thermal stability, low freezing point, low 
shock sensitivity, minimum susceptibility to 
contamination and high compatibility with 
almost all metals and appropriate sealing 
materials. 

These manifold advantages combine to 
make DIMAZINE the outstanding storable 
fuel. We will be pleased to work with you 
in evaluating DIMAZINE and to supply 
detailed data on its properties and handling. 


Work 


FOOD MACHINERY AND CHEMICAL CORPORATION 
Westvaco Chlor-Alkali Division 


General Sales Offices: 
161 E. 42nd STREET, NEW YORK 17 
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Book Reviews 


Combustion and Propulsion, Third 
AGARD Colloquium, edited by M. W. 
Thring, O. Lutz, J. Fabri and A. H. 
Lefebvre, Pergamon Press, New York, 
1958, 610 pp. $20.00. 

Reviewed by A. K. 
University of California 


This impressive volume contains the 
proceedings of the unforgettable Third 
Combustion Colloquium of the Advisory 
Group for Aeronautical Research and De- 
velopment of the North Atlantic Treaty 
Organization, which was held at Palermo, 
Sicily, March 17-21, 1958, and some se- 
lected papers of the 12th Meeting of the 
AGARD Combustion and _ Propulsion 
Panel, which took place the preceding 
year at Washington, D. C. The pre- 
dominating area in this worthy collection 
of papers is the gasdynamic and engineer- 
ing aspects of the subject matter, while 
the purely chemical considerations are 
significantly de-emphasized. 

The proceedings of the Palermo Collo- 
quium consist of five parts. The first— 
concerned with general features of pro- 
pulsion systems—consists of Levine’s dis- 
cussion of problems involved in the de- 
velopment of large liquid rocket engines, 
and of a comprehensive survey, given by 
Lombard, of the overall performance 
characteristics of propulsion by air- 
breathing engines, including some interest- 
ing remarks on the future of transport 
aircraft. 

The second part deals with interaction 
processes and includes the paper of Childs 
on the effect on turbojet combustors and 
afterburners of other engine components; 
an interesting discourse of Grivzs on “ram- 
rockets,” that is, ramjets incorporating 
rocket thrust chambers as the primary 
combustion system; the report by Napo- 
litano, Libby and Ferri on their recent 
work in the field of jet mixing, and the 
paper of Rudinger on the interaction be- 
tween shock waves and flames. 

The third part considers the problem of 
noise—primarily, of course, in the interest 
of its suppression. It is made up of a 
description by Sanders and North of 
methods used for reducing noise of jet 
engines, some thoughts of Richards on 
the design of nozzles for similar purpose, 
and the elucidation of the financial and 
physiological aspects of the problem sup- 
plied by Ducrot, Riehn, Roumilhac and 
Souvras. 

The fourth part is devoted to combus- 
tion problems. It contains a controver- 
sial paper by Spalding, a concise review 
by Gerstein summarizing the present 
knowledge of chemical processes in flames, 
a discussion by Fouré on the fundamental 
features of the development of combustion 
chambers for aircraft propulsion, and a 
report by Schultz, Green and Penner on 
the study of the decomposition mech- 
anism and related aspects of solid com- 
posite propellants. 

The fifth part is titled Aerophysical 
Chemistry; it contains two informative 
surveys: One by Meyerott on the radia- 
tion heat transfer to hypersonic vehicles, 


Ali Bulent Cambel, Northwestern University, Associate Editor 


and the other by Lees on related problems 
involving convective heat transfer, mass 
addition and chemical reactions. This is 
supplemented by a discussion by Laporte 
on the subject of high temperature shock 
waves, some interesting remarks by Schlii- 
ter on magnetohydrodynamics demon- 
strating the advantages of splitting pres- 
sure forces into components parallel and 
perpendicular to the direction of the mag- 
netic field, and the report by Liepmann, 
submitted in the discussion, on magnetic 
effects in high speed boundary layer flow. 

Part VI is made up of four selected 
papers from the 12th AGARD Combus- 
tion Panel, namely, an excellent review of 
flame stabilization by means of gaseous 
jets presented by Cambel, a report on a 
successful application of the crossflow 
jet principle for flameholding by Bertin 
and Salmon, the discussion of flame stabili- 
zation in boundary layers by Tau-Yi 
Toong, and Lutz’s presentation of a 
graphical method for the analysis of gas- 
dynamic mixing. 

Ali in all, the volume forms a worthy 
addition to the distinguished series of 
AGARD publications, forming a valuable 
reference in the field of combustion and 
jet propulsion. 


Some Problems in Chemical Kinetics and 
Reactivity, Volume I, by N. N. Semenov 
(translated from the Russian by Michael 
Boudart), Princeton University Press, 
Princeton, N. J., 1958, xii + 239 pp. 
$4.50. 

Reviewed by S. 8. PENNER 
California Institute of Technology 


Semenov prepared in 1954 an extensive 
survey on chemical kinetics in chain reac- 
tions as a prelude to a symposium held in 
Moscow in 1955. The present book com- 
prises the first of an expanded two volume 
edition based on the 1954 document. 

The distinguished scientist presents in 
this book a readable account of his inter- 
pretation of recent progress on chain reac- 
tions. To western readers, the discussion 
should be of interest not only as a signifi- 
cant contribution from one of the leading 
kineticists of our time, but also as a valu- 
able source of reference to relevant Rus- 
sian publications, many of which have not 
been studied widely by scientists working 
outside of the Soviet Union. 

The five chapters in Volume I include 
treatments on the following topics: Reac- 
tivity of monoradicals (94 pp.), competi- 
tion between monoradical reactions (50 
pp.), reactions of diradicals (22 pp.), dis- 
sociation of molecules and recombination 
of radicals (27 pp.), initiation of chain 
reactions by ions of variable valence (14 
pp.), and wall initiation and termination of 
chain reactions (29 pp.). The treatment 
follows the approach of the classical 
chemical kineticist and contains abundant 
information on the reaction rates of free 
radicals formed from hydrocarbons, as 
well as a brief compilation of kinetics data 
for reactions between simpler chemical 
species. For people who are primarily 


interested in propulsion applications, the 
most valuable part of Volume I may very 
well be an extensive discussion and tabul:- 
tion of bond energies and of reaction rates 
given in chapter I. Particular mention 
should also be made of the treatment of 
heterogeneous reactions given in chapter 
VI, since this discussion covers a field in 
which Russian scientists have made m- 
cently some important basic coatribution-. 

Among the principal attractions of the 
book are the clarity and simplicity of the 
physical descriptions, the insight whic 
the author provides for the mechanism of 
reaction processes, and the panoramic 
coverage on the rates of different classes of 
chemical transformations. The scientifi: 
community owes a debt of gratitude to the 
translator Michael Boudart who has 
rendered a singularly readable account of 
Semenov’s work. 


Aircraft Engines of the World, 1958/59, 
by Paul H. Wilkinson, Paul H. Wilkin- 
son Publications, Washington, D. C., 
1958, 320 pp. $15.00. 

Reviewed by P. Roy CHoupHuRY 
University of Southern California 


“Aircraft Engines of the World” gives « 
good coverage of all the current gas tur- 
bine, ramjet and rocket engines. It also 
lists the major reciprocating engines cur- 
rently being used. A section on nuclear 
aircraft engines acquaints the reader with 
the latest activities in this area of propul- 
sion. Available accessories and equip- 
ment, and fuel and lubrication are also 
included in this volume. 

The author first describes different en- 
gines and their characteristics, and then 
lists them in tabular form for easy refer- 
ence. Russian engines, which are of ma- 
jor interest these days, are also listed. 

In short, the book is an excellent com- 
pilation of figures, tables and engine char- 
acteristics of all the known aircraft pro- 
pulsion systems of the world. 


Books Received 


Liquids and Gases, by Alexaader Efron 
John F. Rider, New York, 1958, 117 pp. 
$1.20. 

Mechanics, by Alexander Efron, John 
F. Rider, New York, 1958, 112 pp. $1.50 

Nuclear Energy, by Alexander Efron, 
John F. Rider, New York, 1958, 72 pp 
$1.25. 

Basic Pulses, by I. Gottlieb, John F. 
Rider, New York, 1958, 176 pp. $3.50. 

Fundamentals of Nuclear Energy and 
Power Reactors, by Henry Jacobowitz, 
John F. Rider, New York, 1959, 144 pp. 
$2.95. 

Fundamentals of Transistors, by Leon- 
ard Krugman, John F. Rider, New York, 
1958, 176 pp. $3.50. 

Gas Tubes, edited by A. Schure, John 
F. Rider, New York, 1958, 80 pp. $1.50. 

Impedance Matching, edited by A. 
Schure, John F. Rider, New York, 1958, 
128 pp. $2.90. 
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JOIN REPUBLIC IN AN 


INTEGRATED ATTACK 
ON PROBLEM AREAS OF 


SPACE EXPLORATION 


It’s the fervent conviction of engineers 
and scientists at Republic Aviation that 
the courageous “Space Time-Table” 
above is entirely feasible — given a tra- 
dition-free, integrated approach to the 


Electronics 


Systems Engineering 
Information Theory 


Inertial Guidance & Navigation 
Digital Computer Development 


problems. Such an approach is evident 
at Republic Aviation. Here, groups of 
specialists from many disciplines are 
working in close collaboration to solve 
problems across the entire spectrum of 
space technologies, which limit today’s 
interplanetary and upper atmosphere 
flight capabilities. 

Expanded by $35,000,000 last: year, 
Republic’s integrated Research and 
Development program has already 
produced signal advances in space guid- 
ance concepts; in new propulsion sys- 


Thermo, Aerodynamics 


Theoretical Gasdynamics 
Hyper-Velocity Studies 
Astronautics Precision Trajectories 


A°TIME-TABLE” 
FOR SPACE 
CONQUEST 


BY 19638 -— INSTRUMENTED 
PLANETARY SOFT LANDING 


BY 1968 SPACE STATION 
FOR STAGING TO MOON AND PLANETS 


BY19'70-'75-—Moon BASE 


These predictions were made by wef. 
Alexander Kartveli, Vice-President for 
Research & Development at Republic fe 
Aviation, and one of the most opti- Fr 
mistic of the 56 leading space experts 
of the world who were consulted by 
the U.S. House of Representatives 
Committee on Astronautics & Space 
Exploration for its report: ‘‘The Next 
10 Years in Space, 1959-1969." 


tems (plasma, nuclear); in radiation 
physics; in new materials and process- 
ing techniques; in unique hypersonic 
configurations; and in prototype devel- 
opment of hardware (as an example: 
hydraulic systems that operate reliably 
up to 1000°F). 

Professional men—who can meet new 
challenges with enthusiasm and dedica- 
tion — are urged to look into openings 
with our R&D groups, working in an 
atmosphere of exhilarating intellectual 
adventure. 


Plasma Propulsion 

Plasma Physics 

Gaseous Electronics 

Hypersonics and Shock Phenomena 
Hydromagnetics 


Telemetry-SSB Technique 

Doppler Radar - Countermeasures 
Radome & Antenna Design 
Microwave Circuitry & Components 
Receiver & Transmitter Design 
Airborne Navigational Systems 
Jamming & Anti-Jamming 
Miniaturization — Transistorization 
Ranging Systems 

Propagation Studies 

Ground Support Equipment 


Airplane/ Missile Performance 

Air Load and Aeroelasticity 

Stability and Controls 

Flutter & Vibration 

Vehicle Dynamics & System Designs 
High Altitude Atmosphere Physics 
Re-entry Heat Transfer 
Hydromagnetics 

Ground Support Equipment 


Physical Chemistry 
Combustion and Detonation 
Instrumentation 

High Power Pulse Electronics 


Nuclear Propulsion 
and Radiation enomena 


Nuclear Weapons Effects 

Radiation Environment in Space 

Nuclear Power & Propulsion Applications 
Nuclear Radiation Laboratories 


Send resume in complete confidence to: Mr. George R. Hickman, Engineering Employment Manager, Dept. 10G 


AVIATION 


Farmingdale, Long Island, New York 
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Information retrieval is a major area of study at IBM. Current investi- 
gations may lead toward such benefits as the instant accessibility to 
knowledge in the Library of Congress —or toward a system which can 
translate any of the earth’s languages into English in real time. 


Problems in information retrieval have defined entirely new concepts 
for the design of storage, input-output and ‘‘memory”’ units —achiev- 
ing far greater capacities than any known today. These facilities will 
provide for the handling of the tremendous amount of updated infor- 
mation needed by business, science and government. With extremely 
rapid accessibility to vast amounts of information electronically stored, 
industrial and research efforts can be materially expedited. IBM needs 
engineers and scientists with the vision and the ability to pave the 
way to tomorrow. 


You will enjoy unusual professional freedom and the support of a 
wealth of systems know-how. Comprehensive education programs are 
available plus the assistance of specialists of many disciplines. Work- 
ing independently or as a member of a small team, your individual 
contributions are quickly recognized and rewarded. This is a unique 
opportunity for a career with a company that has an outstanding 
growth record. 
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TRNAL 


THESE AREAS 


QUALIFICATIONS 


SOME TYPICAL 
ASSIGNMENTS 


JuLy 1959 


CAREERS AVAILABLE IN 


Applied mathematics Operations research 
& statistics Optics 
Circuit design & Programming 
research Real-time engineering 
Component engineering semiconductors 
Human factors solid state development 
Inertial guidance Systems analysis 
Information theory & design 
Logic Transistor device design 


B.S., M.S., or Ph.D. in Electrical or Mechanical 
Engineering, Physics or Mathematics—and proven 
ability to assume a high degree of technical re- 
sponsibility in your sphere of interest. 


704 PROGRAMMER ANALYST to study data flow 
diagrams and write differential equations of a cir- 
cuit diagram; to investigate analog and digital 
real-time control systems, using high-speed elec- 
tronic digital and/or analog computers. Must be 
familiar with variable length alphabetic data, 
transforms, numerical analysis. 


COMPUTER OR SYSTEMS ENGINEER, MATHEMATICIAN 
OR PHYSICIST to design advanced computer, and 
work on development of new information retrieval 
program. Must have strong interest in transistor 
circuit design or in logical or systems applications 
of solid state circuitry. 


MATHEMATICIAN to do programming of information 
retrieval research and plan construction of ad- 
vanced systems. Will play an active part in auto- 
matic programming techniques, numerical analy- 
ses, criteria selection, probability and game theory. 


SENIOR ENGINEER, MATHEMATICIAN OR PHYSICIST 
interested in systems; experienced in operations 
research, communications, missiles or radar. 


For details, write, outlining background and 
interests, to: 

Mr. R. E. Rodgers, Dept. 572G 

IBM Corporation 

590 Madison Ave. 

New York 22, N. Y. 


INTERNATIONAL BUSINESS MACHINES CORPORATION 


New Patents 


George F. McLaughlin, Contributor 


Radiant energy reflector (2,866,971). Kk. 
S. Kellher, Alexandria, Va. 

A semispheric: al metallic reflector placed 
at one side of a dielectric sphere, having a 
variable index of refraction. 

a deploying apparatus (2,869,426). 

J. Wilkie, Superior Township, Wash- 
County, Mich., assignor to Aveo 
Mfg. Corp. 

Apparatus for discharging missile 
from an aircraft. Actuation of an explo- 
sive means within a chamber forcefully 
ejects the missile from spaced retaining 
hooks attached to the aircraft structure. 
Temperature responsive rocket nozzle 
oe T. M. Long, Somerville, 


Rocket body with throat portion of vari- 
able area, Fluid in the chamber, having a 
pressure which is a function of tempera- 
ture, adjusts the throat area. A collap- 
sible tube surrounding the chamber actu- 
ated by combustible gases, seals the cham- 
ber against entry of the gases. 

Variable ejector for iris nozzles (2,870,- 
000). ©. R. Brown, Mission, Kan., as- 
signor to the U.S. Navy. 

Sleeve mounted on a casing for longitu- 
dinal movement, and « movable closure 
overlying a variable area nozzle. A ring 
connected to a bar actuates the closure. 
Afterburner exhaust nozzle variable area 
(2,870,602). R. G. Glenn, Merriam, 
Kan., assignor to the U.S. Navy. 

Leaves forming an annular member, 
hinged to permit inward tilting. Rollers 
on tapered surfaces force the leaves to tilt 
toward the center of the nozzle to reduce 
its area. 

Flow controlling mechanism for reaction- 
type motors (2,870,603). T. M. Long, 
Somerville, N 

Valve and body each provided with a 

bore serving as a passage for the flow of 
fluid oxidizer and fuel. Movement of one 
valve seat into and out of sealing engage- 
ment limits the flow, restricting movement 
of another valve seat. 
Flame stabilizer for high velocity gas 
stream (2,870,604). Kk. W. Conrad, N. 
Royalton, Ohio, assignor to the U. S. 
Navy. 

V-shaped gutter extending transversely 
of the gas stream and having coarse and 
fine screens upstream of the leading edge 
to reduce the local velocity of the stream 
laterally outward from the trailing edges. 
Projectile launcher (2,870,678). P. H. 
Girouard, G. L. Schuyler and M. H. 
Baller, W: ‘ashington, D. C., assignors to the 
Navy. 


Epitor’s Nore: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies of 
patents may be obtained from the Com- 
missioner of Patents, Washington 25, 
D. C., at a cost of 25 cents each; design 
patents, 10 cents 
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Missile guide tube and dummy gun 
mounted for independent movement in 


elevation. Means for bringing the guide 
tube into loading position when empty, 
and into firing position when loaded. 
Missiles are fired in the desired elevation 
angle by the position of the dummy gun. 
Acceleration amplifier (2,870,675). W. 
W. Salisbury, Lafayette, Calif., assignor to 
Zenith Radio Corp. 

Cylindrical coils placed end-to-end to 
form tubular bores. Phase-displaced cur- 
rents and a motor armature comprise a 
conductive projectile movement within a 
motor stator. A localized field traversing 
the stator supports and moves the projec- 
tile at increasing acceleration. 

Missile retaining device (2,870,677). F. 
Herrmann, W. Englewood, N. J., assignor 
to the U.S. Army. 

Mechanism for releasably retaining a 
missile in launching position within a tube. 
A pair of plungers in the housing engage 
bores, urged radially outward by springs. 
A cap screw limits inward movement of the 
plungers. 

Air intake control (2,870,956). H. B. C. 
Dhonau and S. E. Nelson, Indianapolis, 
Ind., assignors to General Motors Corp. 

Blockage segmental flaps pivotally sup- 
ported on links and movable from a posi- 
tion parallel to the wall of a gas compressor 
air intake to a position oblique to the wall 
when segments are moved axially. 
Controls for aircraft (2,870,978). A. A. 
Griffith and A. W. Thorpe, Farnborough, 
England, assignors to Rolls-Royce Ltd., 
and the British Minister of Supply. 

Pair of auxiliary jet nozzles from a gas 
turbine engine spaced apart on an aircraft 
wing to produce mutually opposed turning 
moments on the aircraft. The ratio in 
which the supply is divided between the 


nozzles is varied to alter the effective 
turning moment on the aircraft. 
Method for transmitting intelligence (2,- 
871,463). R. W. Beckwith, N. Syracuse, 
N. Y., assignor to General Electric Co. 
Means responsive to signals of either of 
two frequencies to actuate a respective 
control device. Signals having a fre- 
quency varying regularly between the two 
actuating frequencies prevent actuation of 
the control. 
Sustainer exhaust gas deflector (2,- 
871,658). W.J. Keck, Sacramento, Calif., 
assignor to the U. 8. Navy. 


Safety device for preventing complete 

ignition of a pair of tandem rockets due to 
ignition of either rocket. A conical de- 
flector deflects gases from the fore rocket 
to vent holes, sealing off the forward inter- 
ior part to provide an isolating chamber 
for the fore rocket. 
Jet diverting equipment (2,871,656). W. 
E. P. Johnson, Hempstead, London, Eng- 
land, assignor to Power Jets (Research & 
Development) Ltd. 

Pair of pivotably mounted diversion 
members at the discharge end of the jet 
nozzle. A thermocouple is mounted with 
its hot junction in the air intake; another 
is mounted at a position where heating of 
the aircraft structure by the diverted jet 
ean occur. Linkages control the move- 
ment of the members. 

Explosively released nut (2,871,750). 
N. C. Parrish, Redondo Beach, Calif., 
assignor to Northrop Aircraft, Inc. 

Breakaway attachment and fastening 
element. An explosive charge in a cavity 
may be detonated when desired by means 
of a pair of electrical leads terminating in 
the cavity. 

Devices for measuring the temperature of 
a flame (2,871,759). A. G. L. Moutet, 


reventer (2,- 
C. L. Johnson, Encino, Calif., 
assignor to Lockheed Aircraft Corp. 


Turbine engine blow-out 
870,684). 


_System to prevent interruption of tur- 
bine engine operation by the induction of 
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Rocket Firing 
Timer 


oc 


high temperature gases released upon firing 
large caliber armament. Fuel flow is auto- 
matically reduced at the same rate as the 
air mass flow is reduced, thereby main- 
taining fuel-air ratio to assure stable com- 
bustion. 


Villaine par Massy, France, assignor to 
ONERA Society. 

Measurement of the radiant energy 
emitted by a given area of flame alone, and 
comparison with a source of radiation of 
known and constant energetic brightness. 


Dropping mechanism (2,871,597). P. 
Yost, Hugo, Minn., assignor to Gener:l 
Mills, Ine. 

Load releasing device for balloon. 
Load units spaced around a circular mem- 
ber mounted on a support, and severing 
means for releasing individual load units. 


Rocket (2,872,131). M. E. Lattman, 
Zurich, Switzerland, assignor to Contraves 
AG. 

Rocket guided by a high-frequency sig- 
nal rav. Receiving and transforming dc- 
vices for the signals produce two electric: | 
voltages, the variable instantaneous values 
of which determine the distance of the CC 
of the rocket from the axis of the ray. 


Device for decelerating the sinking speed 
of a dropped load (2,872,138). R. Vogi, 
Pacific Palisades, assignor to Aerophysics 
Development Corp. 

Rocket secured to a support above « 
load, its exhaust directed downward. 
Plates deflect portions of the load later- 
ally. Angles of plates may be varied to 
cause corresponding variation in the mag- 
nitude of the vertical thrust component. 


Rocket and launching mechanism (2,- 
871,762). E. Schmued, Los Angeles, 


Calif., assignor to Northrop Aircraft, Inc 


Means for controlling a finless rocket to 
have maximum stability as it leaves the 
launcher in a spinning condition without 
reducing the propulsive power, or causing 
loss of range. Spin is imparted in the 
launching tube by a source of power apart 
from the rocket. 


Oxidized chromium-alumina metal ce- 
ramic protective tube (2,872,724.) L. A. 
Conant, Tonawanda, N. Y., assignor to 
Union Carbide Corp. 

Sintered ceramic having a combined 
chromium and alumina content constitut- 
ing 62.8 per cent (by weight). To in- 
crease resistance to penetration by molten 
metals, a metal wie having an electrical 
resistance of 100 megohms is coated on 
the surfaces. 


Pulse jet engine with acceleration chamber 
(2,872,780). P. Schmidt, Munich, Ger- 
many. 

Chamber having an outlet end subdi- 
vided by walls aligned with the flow direc- 
tion, and covered with material absorbing 
sound and shock waves. Air admitted at 
the front end passes solely through the 
combustion and acceleration chambers. 


Liquid fuel control for jet engines (2,872,- 

784). O. N. Laurence, Dorridge, Eng- 

a assignor to Joseph Lucas (Industries) 
std. 

Means for preventing a temporary 
excess supply of fuel to an engine during 
acceleration. A liquid operated servo- 
mechanism controls the rate of fuel supply. 
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Technical Literature Digest 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Jet Propulsion Engines 


Modified Tubular Combustors as High- 
temperature Gas Generators, by Robert 
Friedman and Eugene V. Zettle, NACA 
Res. Mem. E55H25, Oct. 1955, 17 pp. 
(Declassified from Confidential by author- 
ity of NASA Pub. Announcement 1, p. 15, 
11 14/58.) 

Experimental Determination of Gas 
Motion Accompanying Screeching Com- 
bustion in a 6-inch Simulated Afterburner, 
by Perry L. Blackshear, Warren D. Rayle 
and Leonard K. Tower, NACA Res. Mem. 
£53128, Dec. 1953, 63 pp. (Declassified 
from Confidential by authority of NASA 
Pb. Announcement 1, p. 12, 11/14/58.) 

Some Screeching-combustion Charac- 
teristics of a Transpiration-cooled After- 
burner Having a Porous Wall of Wire 
Cloth, by William K. Koffel, James L. 
Harp Jr. and Lively Bryant, NACA Res. 
Mem. ¥:54H27, Nov. 1954, 12 pp. (De- 
classified from Confidential by authority 
otf NASA Pub. Announcement 1, p. 13, 
11/14/58.) 

Effects of Variations in Combustion- 
chamber Configuration on Ignition Delay 
in a 50-pound-thrust Rocket, by Dezso J. 
Ladanyi, NACA Res. Mem. 22, Oct. 
1956, 53 pp. (Declassified from Con- 
fidential by authority of NASA Pub. 
Announcement 1, p. 16, 11/14/58.) 

An Analysis of Ramjet Engines Using 
Supersonic Combustion, by Richard J. 
Weber and John 8. MacKay, VACA TN 
4386, Sept. 1958, 49 pp. 

Application of Well Stirred Reactor 
Theory to the Prediction of Combustor 
Performance, by Allan Hardy Bonnell, 
Proj. Squid, Tech. Rep., MIT-18-T (MIT, 
D.Sc. Thesis), 1958, 238 pp. 

Atmosphere Breathing Engines in Astro- 
nautics, by 8S. W. Greenwood and D. 8. 
Carton, Cranfield Coll. of Aeron., Note 
88, Oct. 1958, 17 pp. 

Hydrostatic Gas Bearings, by John H. 
Laub, Calif. Inst. Techn., Jet Prop. Lab., 
Progr. Rep. 20-353, Oct. 1958, 16 pp. 

A Rapid Method for Calculating the 
“Off Design’’ Performance of Compres- 
sors and Turbines, by J. H. Horlock, 
Aeron. Quart., vol. 9, Nov. 1958, pp. 346— 
360. 

Measurements of Rocket Exhaust-gas 
Temperatures, by I. Warshawsky, /nstr. 
Soc. of America J., vol. 5, Nov. 1958, pp. 
91-97. 

Hydro-rockets for Underwater Pro- 
pulsion, by G. Partel, Fusées et Recherche 
Aéron, vol. 3, no. 2. July 1958, pp. 59-64. 
(In French.) 


Heat Transfer and 
Fluid Flow 


Approximate Formula for the Thermal 
Conductivity of Gas Mixtures, by F. A 
Mason and §S. C. Saxena, Phys. Fluids, 
vol. 1, no. 5, Sept.-Oct. 1958, pp. 361-369. 

Properties of Gases at Very High Tem- 
peratures, by I. Amdur and FE. A. Mason, 


Note: Contributions from Pro- 
- sssors E. R. G. Eckert, J. P. Hartnett, T. 
Irvine Jr. and P. J. Schneider of ‘the 
fae Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged 
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Phys. Fluids, vol. 1, no. 5, Sept.-Oct. 1958, 
pp. 370-383. 

Criteria for Thermodynamic Equilib- 
rium in Gas Flow, by Morton Rudin, 
Phys. Fluids, vol. 1, no. 5, Sept.-Oct. 
1958, pp. 384-392. 

Plasma Heating of Supersonic Air- 
streams, by Raymond L. Chuan, Phys. 
Fluids, vol. 1, no. 5, Sept.-Oct. 1958, pp. 
452. 


Some Examples of Turbulent Diffusion 
in Liquid Metals, by A. Paoletti and M. 
Vincentini, Phys. Fluids, vol. 1, no. 5, 
Sept.-Oct. 1958, pp. 453-454. 

Energy Content and Ionization Level in 


an Argon Gas Jet Heated by a High In- | 


tensity Arc, by Gordon L. Cann and 
Adriano C. Ducati, J. Fluid Mech., vol. 4, 
Sept. 1958, pp. 529-541. 

Vaporization into a Hypersonic Laminar 
Boundary Layer, by Sinclaire M. Scala, J. 
Aero/Space Sci., vol. 25, Oct. 1958, pp. 
655-656. 

Magnetogasdynamics of Hypersonic 
Couette Flow, by Z. O. Bleviss, J. Aero/- 
Space Sci., vol. 25, Oct. 1958, pp. 601 615. 

Transport Phenomena in Completely 
Ionized Gas Considering Electron-Elec- 
tron Scattering, by M.S. Sodha and Y. P. 
Varshni, Phys. Rev., vol. 111, no. 5, Sept. 
1, 1958, pp. 1203-1205. 

A Theory of the Cylindrical Ejector 
Supersonic Propelling ot by H. 
i J. B. Holliday and 8. F. Smith, 

J. Roy. Aeron. Soc., vol. 62, Oct. 1958, pp. 
746-751. 

Momentum and Mass Transfer by 
Eddy Diffusion in a Wetted-wall Channel, 
by A. M. Dhanak, A/ChE J., vol. 4, June 
1958, pp. 190-196. 

Flow Against a Vertical Plate with Large 
Suction, by J. R. Foote, J. Aeron. Scz., 
vol. 25, July 1958, pp. 462-463. 

Preliminary Experimental Investigation 
of Transpiration Cooling for an After- 
burner with a Sintered, Porous Stainless 
Steel Combustion Chamber Wall, by W. 
K. Koffel, NACA RM no. E53D08, June 
1953, 47 pp. (Declassified by authority 
of NACA Res. Abstr. no. 128, 5/2/58, p. 


28.) 
Analysis of a Transpiration-cooled 
Hemisphere-cylinder, by C. J. Scott, J. 


Aeron. Sci., vol. 25, June 1958, p. 397. 

Exploratory Tests of Transpiration 
Cooling of a Porous 8° Cone at M = 2.05 
Using Nitrogen Gas, Helium Gas and 
Water as the Coolants, by L. T. Chauvin 
and H.S. Carter, NACA RM no. L55C29, 
June 1955, 22 pp. (Declassified by au- 
thority of NACA Res. Abstr. no. 126, 
5/2/58, p. 35.) 

Transition from Laminar to Turbulent 
Shear Flow—A Review of Some Recent 
Advances in its Understanding, by M. W. 
Morkovin, Trans. ASME, vol. 80, July 
1958, pp. 1121-1128. 

Mixing and Chemical Reaction in the 
Laminar Wake of a Flat Plate, by S. I. 
Cheng and A. A. Kovitz, J. Fluid Mech., 
vol. 4, May 1958, pp. 64-70. 

Partition of Energy in Hydromagnetic 
Turbulence, by S. Chandrasekhar, Annals 
of Phys., vol. 2, Dec. 1957, pp. 615-626. 

Measurements of Total Hemispherical 
Emissivity of Various Oxidized Metals at 
High Temperature, by W. R. Wade, 
NACA TN no. 4206, March 1958, 42 pp. 


MENGINEERS—SCIENTISTS @ 


Work on 


TOTAL 
SOLUTIONS 


to Mayor 
Defense Droblems 


at General Electric’s 
Defense Systems 
Department 


You'll find greater opportunities in 
Systems Engineering in an organi- 
zation whose charter has a breadth 
and scope seldom met in industry: 


..t0 find total solutions to specific 
large scale defense problems requiring 
the integration of diverse fields of 
knowledge and equipments. These 
problems will be of sufficient magni- 
tude and duration to justify the 
allocation by General Electric of con- 
siderable numbers of highly qualified 
scientists and engineers to contribute 
Systems program management and sys- 
tems engineering support—on pro- 
grams such as: 


ATLAS e DYNA-SOAR e SENTRY 
...and other highly classified systems 
which cannot be listed here. 


Immediate Opportunities for: 


Systems Program Engineers 
Systerns Management Engineers 
Guidance Equation Engineers 
Data Processing Engineers 
Electronic Systems Management Engrs. 
Operation Analysts 
Systems Logistics Engineers 
Engineering Psychologists 


Forward your confidential resume at an early 
date. Whereas growth potential here is evi- 
dent—both for DSD and the engineers who 
join us—the positions filled during these 
early months will carry significant ‘‘ground- 
floor’ benefits. 

Write fully to Mr. E. A. Smith 


Dept. 7-M 


Derense Systems Department 
A Department of the Defense Electronics Division 


GENERAL ELECTRIC 


300 South Geddes Street 
Syracuse, New York 
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Composition and Thermodynamic Prop- 
erties of Air in Chemical Equilibrium, by 
W. E. Moeckel and K. C. Weston, NACA 
7'N no. 4265, April 1958, 39 pp. 

Approximations for the Thermodynamic 
and of High Tem- 
perature Air, by C. F. Hansen, NACA TN 
no. 4150, March 1958, 67 pp. 

Boundary Layer Along Annular Walls in 
a Swirling Flow, by H. Yeb, Trans. ASME, 
vol. 80, May 1958, pp. 767-776. 

Some Further Results on the Benard 
Problem, by W. H. Reid and D. L. Harris, 
Phys. Fluids, vol. 1, March-April 1958, pp. 
102-110. 

The Laminar Free-convection Heat 
Transfer from the Outer Surface of a 
Vertical Circular Cylinder, by K. Millsaps 
and K. Pohlhausen, J. Aeron. Sci., vol. 
25, June 1958, pp. 357-360. 

Experimental Study of Thermal Con- 
vection in Rotating Liquid, by R. Hide, 
Roy. Soc., London, Phil. Trans. vol. A 
250, July 31, 1958, pp. 441-478. 

Flow of a Viscous Liquid on a Rotating 
Disk, by A. KE. Emslie, T. F. Bonner and 
L. G. Peck, J. Appl. Phys., vol. 29, May 
1958, pp. 858-862. 

A New Method for the Production of 
Plasma for Ion Sources, by A. Ziegler, 
Zeit. Angew. Phys., vol. 10, April 1958, pp. 
185-186. (In German.) 

A Study of Temperature Measurement 
in a Solar Furnace, by B. B. Brenden, H. 
W. Newkirk Jr. and S. H. Woodcock, 
Solar Energy, vol. 2, Jan. 1958, pp. 13-17. 

On Heat Transfer from the Condensa- 
tion of Steam in the Presence of Inert 
Gases, by J. Madjeski, Chemie-I ngenieur- 
Technik, vol. 29, Dec. 1957, pp. 801-813. 
(In German. ) 

Heat Transfer in Film Condensation, 
by E. Baer and J. M. McKelvey, AIChE 
J., vol. 4, June 1958, pp. 218-222. 

Forced Convection for Gas Flow Normal 
to Circular Cylinder; Influence of Large 
Temperature Difference, by W. M. Kays, 
Stanford Univ., Dept. Mech. Engng. TR 
no. 35, June 1, 1958, 26 pp. 

Transitional Correction to the Drag of a 
Sphere in Free Molecule Flow by R. M. 
L. Baker Jr. and A. F. Charwat, Phys. 
Fluids, vol. 1, March-April 1958, pp. 73- 
81. 

Calculation of the Thermal Conductivity 
of Porous Media, by W. Woodside, Cana- 
dian J. Phys., vol. 36, July 1958, pp. 815- 
823. 

On One-dimensional Temperature Dis- 
tribution in Two-layered Slabs with 
Contact Resistance at the Plane of Con- 
tact, by P. Seide, J. Aeron. Sci., vol. 25, 
Aug. 1958, pp. 523-524. 

Heat Transfer to Water in Turbulent 
Flow in Internally Heated Annuli, by R. 
P. Stein and W. Begell, AJChE J., vol. 4, 
June 1958, pp. 27-31. 

Steady Laminar Heat Transfer in Cir- 
cular Tube with Prescribed Wall Heat 
Flux, by R. Siegel, E. M. Sparrow and 
T. M. Hallman, Appl. Sci. Res., Section 
A, vol. 7, no. 5, 1958, pp. 386-392. 

Heat Transfer in Swirling Laminar 
Pipe Flow, by R. Siegel and M. Perlmut- 
ter, J. Appl. Mech., vol. 25, June 1958, 
pp. 295-297. 

Correlation of Turbulent Heat Transfer 
in a Tube for the Dissociating System 
2NO,, by R. S. Brokaw, VACA 
RM no. E57K19a, March 1958, 17 pp. 

On Thermometer with Recovery Factor 
r Greater Than 1, by J. A. Rietdijk and A. 
Valstar, Appl. Sci. Res., Section A, vol. 7, 
no. 4, 1958, pp. 251-355. 
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On Heat Transfer in Slip Flow, by S. 
H. Maslen, J. Aeron. Sci., vol. 25, June 
1958, pp. 400-401. 

Note on the Stabilizing Effect of Centrif- 
ugal Forces on the Laminar Boundary 
Layer over Convex Surfaces, by L. Lees, 
J. Aeron. Sci., vol. 25, June 1958, pp. 407- 
408. 


On Hypersonic Stagnation-point Flow 
with a Magnetic Field, hy N. H. Kemp, J. 
Aeron. Sci., vol. 25, June 1958, pp. 405- 
407. 

Prandtl-Meyer Expansion of Chemically 
Reacting Gases in Local Chemical and 
Thermodynamic Equilibrium, by S. P. 
Heims, NACA TN no. 4229, April 1958, 
54 pp. 

Evaporation, Heat Transfer, and Veloc- 
ity Distribution in Two-dimensional 
and Rotationally Symmetrical Laminar 
Boundary Layer Flow, by N. Fréssling, 
NACA Tech. Mem. no. 1432, Feb. 1958, 37 
pp. 

Periodic Boundary Layer Experiments 
in Liquid Helium, by R. J. Donnelly and 
A. C. H. Hallett, Annals of Phys., vol. 3, 
March 1958, pp. 320-345. 

Analysis of Turbulent Flow and Heat 
Transfer in a Flat Plate at High Mach 
Numbers with Variable Fluid Properties, 
bv R. G. Diessler and A. I. Loeffler Jr., 
NACA TN no. 4262, April 1958, 61 pp. 

Summary of Experimental Heat Trans- 
fer Measurements in Turbulent Flow for a 
Mach Number Range from 0.87 to 5.05, 
by M. J. Brevoort and B. D. Arabian, 
NACA TN no. 4248, May 1958, 43 pp. 

Effect of Distributed Granular Type 
Roughness on Boundary Layer Transition 
at Supersonic Speeds with and without 
Surface Cooling, by A. L. Braslow NACA 
RM no. L58A17, March 10, 1958, 22 pp. 

The Effect of Gas Properties on the 
Heat Transfer in Stagnation Flows, by I. 
KE. Beckwith, J. Aeron. Sci., vol. 25, Aug. 
1958, pp. 533-534. 

Fundamental Solutions for Heat Trans- 
fer from Nonisothermal Flat Plates, by 
Donald C. Baxter and W. C. Reynolds, 
J. Aeron. Sct., vol. 25, June 1958, pp. 403- 
404. 

New Method of Flow Visualization for 
Low-density Wind Tunnels, by R. A. 
Evans, J. Appl. Phys., vol. 28, Sept. 1957, 
pp. 1005-1010. 

An Approximate Analytical Method for 
Studying Entry into Planetary Atmos- 
pheres, by D. R. Chapman, NACA TN 
no. 4276, May 1958, 101 pp. 

Two-dimensional Effects in Gaseous 
Detonation Waves, by James A. Fay and 
Gerhard Opel, J. Chem. Phys., vol. 29, 
Oct. 1958, pp. 955-956. 


Combustion, Fuels and 
Propellants 


Effect of Pressure and Duct Geemetry 
on Bluff-body Flame Stabilization, by 
Andrew E. Potter Jr. and Edgar L. Wong, 
NACA TN 4381, Sept. 1958, 31 pp. 

Identification of Ions in Flame and 
Measurement of Their Concentration, by 
A. van Tiggelen, J. Deckers, J. Poncelet 
and R. Berendsen, Univ. of Louvain, Lab. 
Analyt. Chem., Tech. 
(Scientific) Note 1, March 1958, 43 pp. 

Preliminary Study of the Light Scatter- 
ing Technique for Determination of Size 
Distributions in Burning Sprays, by D. A. 
Dobbins, Princeton Univ., Dept. Aeron. 
Enang., Rep. 430 (AFOSR TN 58-822; 
ASTIA AD 202906), July 1958, 54 pp. 

Some Effects of High-energy Fuels on 


Aircraft Performance, by Abraham Hyatt, 
Aero/Space Engng., vol. 17, Nov. 1958, 
pp. 45-50. 

Thermal Conductivity cf Liquid Ozone, 
by Thomas E. Waterman, Donald P. 
Kirsh and Robert I. Brabets, J. Chen. 
Phys., vol. 29, Oct. 1958, pp. 905-908. 

Present Status of Detonation Theory, 
by W. W. Wood and J. G. Kirkwood, //. 
Chem. Phys., vol. 29, Oct. 1958, pp. 957- 
958. 

Precision Flash X-Ray Determination 
of Density Ratio in Gaseous Detonations, 
by Russell E. Duff, Herbert T. Knight 
and John P. Rink, Phys. Fluids, vol. 1, no. 
5, Sept.-Oct. 1958, pp. 393-398. 

Transition from Slow Burning to Deto- 
nation in Gaseous Explosives, by Fred: - 
rick J. Martin, Phys. Fluids, vol. 1, no. 5, 
Sept.-Oct. 1958, pp. 399-407. 

Effect of Pressure on Ion Formation in 
Propane Air Flames, by I. R. King, Pro,. 
Squid, Tech. Rep. EXP-7-P (ASTIA Al) 
159256), May 1958, 5 pp. (Available 
only on microcards. ) 

A Study of the Kinetics of the Hydrogen 
Oxygen Reaction in a New Flow Reactor, 
by Rudolph J. Swigart, Princeton Univ., 
Dept. Aeron. Engng., Rep. 432, Aug. 1958, 
58 pp., 23 figs. (MSE Thesis. ) 

Carbon Disulfide Flames Influence on 
Moisture, by G. J. Gibbs, I. R. King, and 
H. F. Calcote, Experiment Inc., CF 2264 
(TM-642), May 1954, 27 pp. 

Effect of Fuel Drop Size and Injector 
Configuration on Screaming in a 200- 
pound-thrust Rocket Engine Using Liquid 
Oxygen and Heptane, by Charles E. Feiler, 
NACA Res. Mem. F58A20a, June 1958, 
27 pp. (Declassified from Confidentia! 
by authority of NACA Res. Abstr. 130, p. 
25, 9/30/58. ) 

Effect of Fuels on Screaming in 200- 
pound-thrust Liquid-oxygen-fuel Rocket 
Engine, by Isaac Pass and Adelbert O. 
Tischler, NACA Res. Mem. E56C10, 
June 1956, 25 pp. (Declassified from 
Confidential by authority of NACA Res. 
Abstr. 130, p. 25, 9/30/58. ) 

Application of Well Stirred Reactor 
Theory to the Prediction of Combustor 
Performance, by Allen Hardy Bonnell, 
Proj. Squid, Tech. Rep. MIT-18-T (AS- 
TIA AD 200821), Aug. 1958, 58 pp. 
(Available only on microcards. ) 

The Effects of Several Composition 
Factors on the Burning Rate of an Am- 
monium Perchlorate Solid Propellant, by 
Harold J. Taback, Princeton Univ., Dept. 
Aeron. Engng., Rep. 429, Sept. 1958, 68 
pp., 34 figs. (MSE Thesis.) 

Aluminum Borohydride-Hydrocarbon 
Mixtures as a Source of Ignition for a 
Turbojet Combustor, by Hampton H. 
Foter, Edward A. Fletcher and David M. 
Straight, NACA Res. Mem. F54K12, Feb. 
1955, 24 pp. (Declassified from Confi- 
dential by authority of NACA Res. Abstr. 
130, p. 24, 9/30/58. ) 

Introductory Study to Aerothermcchem- 
istry, by L. G. Napolitano, L’ Aerotecnica, 
vol. 38, no. 2, April 1958, pp. 96-108. (In 
Italian. ) 


Kinetics of the Thermal Reactions of 
Ethylene, by George Dahlgren Jr. and 
John E. Douglas, J. Amer. Chem. Soc., 
vol. 80, Oct. 5, 1958, pp. 5108-5110. 

Numerical Calculation of the Tempera- 
ture Rise in Ignition by Hot Wires, by 
H. L. Selberg and C. H. Johansson, Arkiv 


fiir Fysik, vol. 13, no. 5, 1958, pp. 423-427. 


(In English. } 

Hydroxyl Concentrations in Rich Hy- 
drogen-air Flames Held on Porous Bur- 
ners, by W. E. Kaskan, Combustion and 
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Another important advance in man’s 
knowledge of outer space was provided 
by Pioneer III. This, like many others of a 
continuing series of space probes, was 
designed and launched by Jet Propulsion 
Laboratory for the National Aeronautics 
and Space Administration. JPL is admin- 
istered by the California Institute of 
Technology for NASA. 
and During its flight of 38 hours, Pioneer III 
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PIONEERING IN SPACE 


CALIFORNIA 


was tracked by JPL tracking stations 
for 25 hours, the maximum time it was 
above the horizon for these stations. 
The primary scientific experiment was 
the measurement of the radiation en- 
vironment at distances far from the Earth 
and telemetering data of fundamental 
scientific value was recorded for 22 hours. 
Analysis of this data revealed, at 10,000 
miles from the Earth, the existence of a 


PASADENA, CALIFORNIA 


Bur- OPPORTUNITIES NOW OPEN APPLIED MATHEMATICIANS . ENGINEERING PHYSICISTS . COMPUTER ANALYSTS 
AERONAUTICAL ENGINEERS . RESEARCH ANALYSTS . 


DESIGN ENGINEERS 


RESEARCH 


belt of high radiation intensity greater 
than that observed by the Explorer 
satellites. 

This discovery is of vital importance as 
it poses new problems affecting the dis- 
patch of future vehicles into space. The 
study and solution of such problems 
compose a large part of the research and 
development programs now in extensive 
operation at the Laboratory. 
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Flame, vol. 2, Sept. 1958, pp. 229-243. 

The Spontaneous Ignition of Alkyl Sil- 
anes, by S. T. Griffiths and R. R. Wilson, 
Combustion and Flame, vol. 2, Sept. 1958, 
pp. 244-252. 

The Flashback of Laminar and Turbu- 
lent Burner Flames at Reduced Pressure, 
by B. Fine, Combustion and Flame, vol. 2, 
Sept. 1958, pp. 253-266. 

Formation of Polycyclic Hydrocarbons 
in Diffusion Flames, b R. Arthur, B. 
T. Commins, J. A. 8. Gilbert, A. J. Lind- 
sey and D. H. Napier, Combustion and 
Flame, vol. 2, Sept. 1958, pp. 267-272. 

Use of the Constant-volume Bomb 
Technique for Measuring Burning Veloc- 
ity, by R. C. Eachenbach and J. T. Ag- 
new, Combustion and Flame, vol. 2, Sept. 
1958, pp. 273-285. 

The Concentration of Hydroxyl and of 
Oxygen Atoms in Gases from Lean 
Hydrogen-air Flames, by W. E. Kaskan, 
Combustion and Flame, vol. 2, Sept. 1958, 
pp. 286-304. 

A Note on Spalding’s Centroid Formula, 
by J. Menkes, Combustion and Flame, vol. 
2, Sept. 1958, pp. 328-330. 

Note on Monocromatic Radiation Py- 
rometry for Measuring Flame and Exhaust 
Gas Temperatures, by R. H. Tourin and 
M. Grossman, Combustion and Flame, vol. 
2, Sept. 1958, pp. 330-332. 

Rocket Working Fluids. Hydrogen and 
Helium, by R. L. Potter, G. E. Guffner 
and T. F. Reinhardt, Ind. & Engng. 
Chem., vol. 50, Oct. 1958, pp. 1557-1560. 

Screw Extrusion Theory with Applica- 
tion to Double-base Propellant, by M. L. 
Jackson, F. J. Lavacot and H. R. Rich- 
ards, Ind. & Engng. Chem., vol. 50, Oct. 
1958, pp. 1569-1576. 

Ignition Limits of the Gaseous System 
Isopropyl Alcohol-Oxygen-Nitrogen, by 
Harry Sello, Ind. & Engng. Chem., vol. 50, 
Oct. 1958, pp. 1561—1564. 

Generalized Thermodynamic Proper- 
ties of Diatomic and Triatomic Gases, by 
B. E. Saltzman, Ind. & Engng. Chem., vol. 
50, Oct. 1958, pp. 1593-1598. 

Gaseous Detonations, XI, Double 
Waves, by Mark Cher and G. B. Kistia- 
kowsky, J. Chem. Phys., vol. 29, Sept. 
1958, pp. 506-511. 

The Development of a New Erosive 
Burning Law, by G. Robillard and J. M. 
Lenoir, Calif. Inst. Tech., Jet. Prop. Lab., 
Publication 91, April, 1957, 27 pp. 

Screaming Tendency of the Gaseous- 
hydrogen-liquid-oxygen Propellant Com- 
bination, by Louis Baker Jr. and Fred W. 
Steffen, NACA Res. Mem. E58E09, Sept. 
1958, 24 pp. 

Phase Behavior in the Nitrogen Am- 
monia System, by H. H. Reamer and B. 
H. Sage, Calif. Inst. Tech., Jet Prop. Lab., 
Progr. Rep. 20-358, July 1958, 6 pp. 

Temperature Gradient Tests of Inert 
Propellant Grains, by James F. Wiedman, 
Redstone Arsenal, Ordn. Missile Labs., 
Rep. 3M78N 10, Aug. 1958, 4 pp, 23 figs. 

Determination of Temperature Gra- 
dients in Solid Propellant Grains by Use of 
Scale Models, by Myron L. Cohen, Red- 
stone Arsenal, Ordn. Missile Labs., Rep. 
3R5F, July 1958, 67 pp. 

The Thermal Decomposition of Liquid 
Nitric Acid, II, by H. F. Cordes, N. R. Fet- 
ter and J. A. Happe, J. Amer. Chem. Soc. ; 
vol. 80, Sept. 20, 1958, pp. 4802-4805. 

Free Radicals for High Energy Fuels, 
by D. E. Carr, Fox and D. 
Guth, Aviation ‘Age, vol. 30, Oct. 1958, pp. 
22-23, 191. 

Homogeneous Catalysis for Homoge- 
neous Reactors, Catalysis of the Reaction 
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Between Hydrogen and Oxygen, by H. F. 
McDuffie, E. L. Compere, H. H. Stone, 
L. F. Woo and C. H. Secoy, J. Phys. 
Chem., vol. 62, Sept. 1958, pp. "1030-1035. 

An ‘Analytical Study of Turbulent and 
Molecular Mixing in Rocket Combustion, 
by David A. Bittker, NACA TN 4821, 
Sept. 1958, 22 pp. 

Mass Analysis of Flames and Flue 
Gases, by E. A. Bunt, Univ. of Witwaters- 
rand, Dept. Mech. Engng., Rep. 2758, 
June 1958, 43 pp. 

Propellant Vaporization as a Criterion 
for Rocket-engine Design; Calculations 
of Chamber Length to Vaporize Various 
Propellants, by Richard J. Priem, NACA 
TN 3883, Sept. 1958, 36 pp. 

Study of Explosive Sensitivity of Cobalt 
Amine Complexes, by Taylor B. Joyner, 
NAVORD Rep. 5639, Oct. 1957, 13 pp. 

Oxidation of 1,1-Dialkylhydrazines- 
chromatographic Separation and Ultra- 
violet Analysis of Some Tetraalkyltetra- 
zenes, by Everett M. Bens and William 
R. McBridge, NAVORD Rep. 5879 
(NOT'S 2016), June 1958, 34 pp. 

Molecular Structure and Infrared Spec- 
trum of Nitrotrichloromethane, by Alex- 
ander H. Castelli, Picatinny Arsenal, 
Feltman Res. and Engng. Labs., Tech. Rep. 
2537, Aug. 1958, 13 pp. 

A High Rate Tension Tester, by Ste- 
phen Strella, Mitchel Chmura, Benjamin 
Holman and Harry Sigler, Picatinny 
Arsenal, Feltman Res. and Engng. Labs., 
Tech. Rep. 2487, Aug. 1958, 12 pp. 

Rate of Reaction of Gaseous Fluorine 
with Water Vapor at 35°C, by Vernon A. 
Slabey and Edward A. Fletcher, NACA 
TN 4374, Sept. 1958, 16 pp. 

Study of Combustion by Cinespectros- 
copy (Plenary Conference), by R. G. W. 
Norrish, Inst. Francais du Petrole, Revue 
et Ann. Combustibles Liquides, vol. 8, 
April 1958, pp. 313-332. (In French.) 

Inhibition in the Gaseous Phase of 
Combustion of Carbon Oxides, I: The 
Phenomenon of Flames Long Retarded 
and Its Application to the Determination 
of the Kinetic Constants of Combustion 
Reactions, by H. James; II: The Charac- 
ter of the Phenomenon of Inhibition by 
Ammonia, by H. James and P. Lafitte, 
Inst. Francais du Petrole, Revue et Ann. 
Combustibles Liquides, vol. 8, April, pp. 
338-358. (In French.) 

A Method cf Extracting Ions in Flames, 
by J. Deckers and A. Van Tiggden, Jnst. 
Francais du Petrole, Revue et Ann. Com- 
bustibles Liquides, vol. 8, April 1958, pp. 
349-358. (In French.) 

Study of Oxidation of Hydrocarbons in 
Extremely Rigorous Experimental Con- 
ditions, by Z. G. Szabo and D. Gal, Inst. 
Francais du Petrole, Revue et Ann. Com- 
bustibles Liquides, vol. 8, April 1958, pp. 
414-427. (In French.) 

Thermal Aspects of the Combustion at 
Low Temperature of n-Pentane-Oxygen 
Mixtures, by R. Ben Aim and M. Luc- 
quin, Inst. Francais du Petrole, Revue et 
Ann. Combustibles Liquides, vol. 8, April 
1958, pp. 438-442. (In French.) 

Flames of Mixtures of Hydrazine with 
Oxygen, Nitrous Oxide and Nitric Oxide, 
by P. Gray and J. C. Lee, Inst. Francais 
du Petrole, Revue et Ann. Combustibles 
Liquides, vol. 8, April 1958, pp. 443-450. 
(In French. ) 

Flame Propagation in Mixtures of 
Hydrazine, by S. deJaegere and A. Van 
Tiggelen, Inst. Francais du Petrole, Revue 
et Ann. Combustibles Liquides, vol. 8, 
April 1958, pp. 451-456. 

Analysis of the Decomposition Products 
of Hydrazine, by G. Pannethier, H. Guene- 


baut and R. de Hartoulari, Inst. Francais 


du Petrole, Revue et Ann. Combustibles 
Liquides, vol. 8, April 1958, pp. 457-465. 
(In French. ) 

Burning of Gaseous Combustibles with 
Nitric Oxide, by G. K. Adams and G. W. 
Stocks, Inst. Francais du Petrole, Revue et 
Ann. Combustibles Li utdes, vol. 8, April 
1958, pp. 483-490. ( n French.) 

Burning-rate Studies, Part 7: Onset of 
Turbulent Combustion of Liquids Con- 
tained in Small Tubes, by A. yet le 
Whittaker, T. M. Donovan and H. W)!- 
liams, J. Phys. Chem., vol. 62, no. 8, Aug. 
1958, pp. 908-911. 

The Theoretical Study of the Florine 
Molecule, by J. Eve, Proc. Roy Soc., vol. 
246A, no. 1247, Aug. 26, 1958, pp. 582- 
591. 

Problems of Continuous Combustion 
in Aviation Engines, by A. Capetti, L - 
Aerotecnica, vol. 38, no. 1, Feb. 1958, p). 
3-7. (In Italian.) 

Stabilization of Flames in Airbreathinz 
Jet Engines, by F. Filippi, L’ Aerotecnic:, 
vol. 38, no. 1, Feb. 1958, p. 19-28. (In 
Italian.) 

Liquid Propellants for Rockets, by |. 
Macioce, L’ Aerotecnica, vol. 38, no. }, 
Feb. 1958, p. 44-47. (In Italian.) 

The Indication of the Specific Fuel Con- 
sumption of Turbojets, by P. Formen- 
tini, L’ Aerotecnica, vol. 38, no. 1, Feb. 
1958, pp. 39-43. (In Italian.) 

Theoretical Investigation on the Fluid 
Dynamics of ‘‘Wave Engines’’ Operating 
according to the Principle of the ‘‘Com- 
prex,’? by Dino Dini, L’ Aerotecnica, vol. 
38, no. 1, Feb. 1958, pp. 8-18. (In Italian. 


Aerodynamics of Jet 
Propelled Missiles 


Drag Transformation and Reduction for 
Bodies of Revolution, by Roger P. John- 
son, Rand Corp., Res. Mem. RM-2107 
(ASTIA AD 156008), Aug. 1957, 73 pp. 

Life of Slender Nose Shapes According 
to Newtonian Theory, by J. D. Cole, Rand 
Corp., Paper P-1270, Feb. 1958, 18 pp. 

On the Problem of a Streamlined Pro- 
file in a Near Sonic Flow, by A. F. Kryu- 
chin, Rand Corp., Translation T-74 (trans- 
lation from Acad. Scei., Inst. 
Mech., Appl. Math. and Mech., vol. 18, 
1954), Aug. 1957, 25 pp. 

The High Speed Flow of Gas around 
Blunt Bodies, by Hyman Serbin, /tand 
Corp., Paper P-1172, Sept. 1957, 31 pp. 


Materials of Construction 


New Plasma Tools Handle Ultra-hard 
Metals, by Peer Fossen, Missiles and 
Rockets, vol. 4, Dec. 22, 1958, pp. 26-27. 

Falcon Airframe Uses Die Casting, 
Magnesium Structure, by Irwin Stambler, 
Space/Aeron., vol. 30, Dec. 1958, pp. 18- 
19. 

Aeroelastic Problems Associated with 
High Speeds, and High Temperature, by 
E. G. Broadbent, J. Roy. Aeron. Soc., 
vol. 62, Dec. 1958, pp. 867-883. 

Methods of Protecting Reaction Motor 
Materials under Unstable Combustion 
Conditions, by B. Langenecker, Fusées et 
Récherche Aéron., vol. 3, no. 2, July 1958, 
pp. 77-85. (In French.) 

Thermal Stresses in Plates, by M. J. 
Forray, J. Aero/Space Sci., vol. 25, Nov. 
1958, pp. 716-717. 

Behavior of Transparent Plastics under 
Thermal Shock, by Kh. Iablokoff and M. 
Hediard, La Récherche Aéron., no. 66, 
Sept.-Oct. 1958, pp. 35-42. (In French.) 
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Space Technology Laboratories new corporate symbol represents a bright history in a stimulating age. STL has 
provided the over-all systems engineering and technical direction for the Air Force Ballistic Missile Program since it 
ipas assigned the highest national priority in 1954. Five years of accelerated effort produced epic advances in science 
und technology, and propelled the art of missilery through three distinct generations of progress. STL. contributed 
‘echnical leadership to the science/government/industry team which has built this solid, expandable foundation 
jor future advances in space, and is daily adding new strength to our national security. * In addition to its major 


management functions, STL also conducts advanced space probe experiments for the Air Force at the direction of 


such agencies as NASA and ARPA. * To those scientists and engineers with capabilities in propulsion, electronics. 
thermodynamics, aerodynamics, structures, astrophysics, computer technology, and other related fields and 


disciplines, STL now offers unique professional opportunities. Inquiries regarding staff positions at STL are invited. 


a new symbol 
for a new era of 
technology 
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ENGINEERS * SCIENTISTS 


Seek and find 
Daring New Goals 
and Rewards with 
General Electric’s 
Growing 


| 
FLIGHT PROPULSION 


DIVISION IN CINCINNATI | 


(*formerly Aircraft Gas Turbine Di- 
vision, recently renamed to match 
its expanding product goals.) 


Our engineers and scientists enjoy 
the daily challenge of new and stim- 
ulating technical problems in the 
largest number of advanced propul- 
sion system projects in our history. 
Men with imagination and ideas are 
r as Managers and Specialists 
in: 


...high Mach powerplants using 
high energy fuels (Jet Engine Dept.) 


-...commercial jet powerplants 
Commercial Engine Operation) 


+++.@xotic propulsion systems for 
rockets and space vehicles using 
plasma and ion power sources | 
Flight Propulsion Laboratory Dept.) 


..+.-developing production tech- 
niques for unusual materials and 
design configurations (Production 
Engine Dept.) 


Advancement opportunities 
are multiplied— 


nearly 50% of our professional 
people earn promotion every year. 
Our engineering staff has more than 
doubied since 1955. 


LOOK INTO NEW POSITIONS 

NOW OPEN IN OUR 4 DEPART- 
MENTS IN CINCINNATI for Engineers 
and Scientists experienced in: 


Rocket Motor Prelim. Design 
Liquid Rocket Applications 
Solid Rocket Casings 

Nozzle Design 

lon Acceleration 

Plasma Generation 

Advanced Space Applications 


(U. S. Citizenship Required) 


For further information on the oppor- 
tunities that will let you use your tal- 
ents better, send your resume to Dr. 
Mark Peters. Bldg. 100, Dept. 35-MG 


FLIGHT PROPULSION DIVISION 


GENERAL ELECTRIC 


Cincinnati, Ohio 


Instrumentation, Data 
Recording, Telemetering 


The Constant Current Hot Wire Ane- 
mometer in Subsonic Compressible Flows, 
by D. I. Paddison, Brown Univ., Div. 
Engng., Tech. Rep. WT-29 (AFOSR-TN- 
58-713; ASTIA AD 162248), 1958, 37 
pp. 

Chronograph Error in Drag Measure- 
ments, by G. D. Kahl and F. D. Bennett, 
Aberdeen Proving Ground, Ball Res. Labs., 
Mem. Rep. 1161, Aug. 1958, 15 pp. 


Static Calibration of Pressure and 
Force Gages (Joint Army Navy Air Force 
Solid Propellant Static Test Panel), Solid 
Propellant Info. Agency, Aug. 1958, 28 
pp. 

High Speed Framing Disc Camera, by 
C. H. Bagley, Stanford Res. Inst., Poulter 
Labs., Tech. Rep. 015-58, April 1958, 11 
pp. 

New Possibilities in the Design of 
Vibration Measuring Apparatus, by C. V. 
Basel, Archiv fiir Tech. Messen, no. 273, 
Oct. 1958, pp. 201-204. (In German.) 

The Application of a Magnus Effect 
Flow Meter as a Laboratory Teaching 
Device, by Wesley J. Trathen, Gen. Motors 
Engng. J., vol. 5, Oct.-Dec. 1958, pp. 12- 
15. 


Manometric Pressure Measurement, by 
R. J. Cartier, Instr. and Automation, vol. 
31, Dee. 1958, pp. 1980-1982. 

Electronic Density Probe, by William 
G. Burt Jr., Instr. and Automation, vol. 31, 
Dec. 1958, pp. 1983-1984. 

Proceedings, Conference on Radio 
Interference Reduction, 3d, Feb. 26-27, 
1957, Illinois Inst. Tech., Armour Res. 
Found. (sponsored by Army Signal 
Engng. Labs.), Feb. 1957, 503 pp. 

Pressure Probes Selected for Three 
Dimensional Flow Measurements, by 
D. W. Bryer, D. E. Walshe and H. C. Gar- 
ner, Gt. Brit., Aeron. Res. Council, Rep. 
and Mem. 3037 (formerly ARC Tech. 
Rep. 3037), 1958, 26 pp. 

A Survey of New Developments in 
Pressure Measuring Techniques in the 
NACA, by John Dimeff, NATO, AGARD, 
Rep. 166, March 1958, 25 pp. 

Tests on Deflection of Flat Oval Bour- 
don Tubes, by G. Kardos, ASME, Paper 
58-A-67, Nov. 1958, 6 pp 

Theory of High Pressure Bourdon 
Tubes, by W. Wuest, ASME, Paper 58- 
A-119, Nov. 1958, 12 pp. (Ing. Arch., 
vol. 24, 1956, pp. 92-110.) 

Pressure Pick Ups and Miniature 
Manometers Usable in Wind Tunnels and 
in Flight, by M. Bassiére, NATO, 
AGARD, Rep. 173, March 1958, 18 pp. 
(In French.) 

Electromagnetic Flowmeter Primary 
Elements, by V. P. Head, ASME, Paper 
58-A-126, Nov. 1958, 7 pp. 

Multi-variable Function Generator, by 
Thomas J. Kelley, Jnstr. and Automation, 
vol. 31, Nov. 1958, p. 1819. 

Who Uses Medium-Size Computers, 
by Phyllis Huggins, Jnstr. and Auto- 
mation, vol. 31, Nov. 1958, pp. 1827- 
1831. 

Mathematical Elements for Analog 
Computers, by Charles Lax and Phillip I. 
Rafield, Instr. and Automation, vol. 31, 
Nov. 1958, pp. 1832-1833. 

Pressure Transducers for Missile Test- 
and Control, by Yao Tzu Li, JSA J. 
(Instr. Soc. America), vol. 5, Nov. 1958, 
pp. 81-84. 

Evolution of Test-range Instrumenta- 
tion, by Dr. Joachim W. Muehlner, [SA 


J. (Instr. Soc. America), vol. 5, Novy. 
1958, pp. 85-90. 

Calibration of High-frequency-response 
Pressure Transducers, by Ralph Bower- 
sox, JSA J. (Instr. Soc. America), vol. 5, 
Nov. 1958, pp. 98-103. 

Magnetic Tape Storage of Telemetry 
Data, by Kenneth A. Hall, JSA J. (Instr. 
Soc. America), vol. 5, Nov. 1958, pp. 104- 
106. 

Optical Instrumentation—Obsolescence 
or Revival, by Alexis B. Dember, 
(Instr. Soc. America), vol. 5, Nov. 195s, 
pp. 107-113. 

Review of High-temperature Immersion 
Thermal Sensing Devices for In-flight 
Engine Control, by Victor D. Sander-, 
Rev. Sci. Instr., vol. 29, Nov. 1958, p). 
917-928. 

High Speed Multiple-spark Light 
Source, by Merle R. Wilson and Richard J. 
Hiemenz, Rev. Sci. Instr., vol. 29, Nov. 
1958, pp. 949-950. 

Microcomparator Illuminations Systeni 
for Improved Resolution of Rayleigh Dif- 
fusionmeter Diagrams, by J. L. Oncley 
M. Ludwig and T. E. Thompson, Rev. Sc: 
Instr., vol. 29, Nov. 1958, pp. 985-987. 

Convection Free Instrument for Meas- 
uring Infrared Radiation in the Atmos- 
phere, by Robert L. Aagard, Rev. Sc: 
Instr., vol. 29, Nov. 1958, pp. 1011-1015. 

Temperature Error Associated with 
Imbedded Thermocouples, by A. R 
Thomas, B. Schurin and J. C. Morris, 
Rev. Sci. Instr., vol. 29, Nov. 1958, pp. 
1045-1046. 

Pressure Calibrator for Field Use, by 
Earl T. Pearson, Rev. Sci. Instr., vol. 29, 
Nov. 1958, pp. 1052-1053. 

Tube Diameter and Thickness Gauge 
for High Temperatures, by Ian Finnie, 
Rev. Sci. Instr., vol. 29, Nov. 1958, p. 
1054. 


Guidance Systems and 
Components 


Inertial Guidance for Hypersonic and 
Orbital Vehicles, by George W. Johnson 
and R. B. MacDonald, Jnst. Aeron. Sci., 
Rep. 59-33, Jan. 1959, 23 pp. 

Guidance Requirements of Ballistic 
Satellites, by Murray J. Stateman, Jnsf. 
Aeron. Sci., Rep. 59-30, Jan. 1959, 16 
pp. 


Servomechanisms and 
Controls 


An Introduction to the Time Modulated 
Acceleration Switching Electrohydraulic 
Servomechanism, by Stephen A. Mur- 
taugh Jr.,. ASME, Paper 58-A-159, Nov. 
1958, 6 pp. 

A Weight Comparison of Several Atti- 
tude Controls for Satellites, by James J. 
Adams and Robert G. Chilton, NASA 
Memo 12-30-58L, Feb. 1959, 19 pp. 

Design of Automatic Flight Control 
Systems for Manned Space Vehicles, by 
Albert J. Kelly, Znst. Aeron. Sci., Rep. 
59-32, Jan. 1959, 14 pp. 

A High Resolution Servo, by W. A. 
Good, Wissenschaftlichen Gesellschaft fiir 
Luftfahrt, Jahrbuch, 1957, Braunschweig, 
Friedr. Vieweg & Sohn, 1958, pp. 368- 
370. 


Space Flight 
A Note on ‘‘A Nuclear Electric Pro- 
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Now from Du Pont: 


a complete line of photorecording papers and films 


Du Pont now makes photorecording papers and films for 
your various recording needs. Here's what our complete 
photorecording line includes: 

PHOTORECORDING PAPERS 
Lino-Writ 1. Orthochromatic paper designed for average test 
requirements at lower writing speeds. 
Lino-Writ 2. Moderately high-speed orthochromatic paper 
for use in intermediate speed recording requirements, 
Lino-Writ 3. For high-frequeney oscillographic traces with- 
out loss of detail. 
These three papers are available in two types: Type B. stand- 
ard-weight stock when greater opacity is required, and Type 
W, a thin, 100% all-rag stock. 
Lino-Writ 4. It’s the toughest, fastest. whitest paper you can 
use; 20% thinner than Type W. accommodates writing speeds 
as high as 5000 cps at maximum amplitude. Exceptional 
exposing and processing latitude. All-rag stock has unusual 
wet or dry strength. 


*CRonar is Du Pont's trademark for its polyester photographic films 


This advertisement was prepared exclusively by Phototypography 
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PHOTORECORDING FILMS 

CRONAR* Recording Film (High Speed) —relative speed 
(tungsten): 30. 

Lino-Flex 1 (Regular Speed) —relative speed (tungsten): 10. 
Both of these films are on Du Pont’s exclusive CRONAR 
polvester photographic film base. which otfers unexcelled 
strength, exceptional dimensional stability and flexibility. 
These fast-drying films have a fine matte surface which 
aceepts peneil and ink markings. 


PHOTORECORDING CHEMICALS 


Lino-Writ Rapid Processing Chemicals Kit. Du Pont’s 
high-speed, low-odor chemicals designed for rapid stabiliza- 
tion processing of papers and films. 


For more information on our complete photorecording line, 
write: E. I. du Pont de Nemours & Co. (Inc.). Photo Prod- 
ucts Department. Wilmington 98, Delaware. In Canada: 
Du Pont Company of Canada Limited, Toronto. 


REG. U.S, Pat. OFF. 
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ENGINEERS 

& SCIENTISTS 

\ 

\ Aeronutronic Systems, 
Inc., a dynamic new 
name in science and 
research, has immediate 

\ need for qualified people Y 
to staff senior positions 
at its new Research 
Center in Newport 
Beach, Southern Cali- 

\ fornia. 

The Space Technology Y 

Division of this rapidly 


expanding Ford Motor 
Company _ subsidiary 
has career openings in 
the following fields of 


interest: 


VEHICLE TECHNOLOGY 

Aerodynamic design 
and testing 

Rocket engine develop- 
ment 

Rocket nozzle and re- 
entry materials 

High temperature 
chemical kinetics 

Combustion and deto- 
nation theory 

Combustion thermody- 
namics 

High temperature 
structural plastics & 
ceramics 

Advanced structures 

Rocket vehicle systems 


MISSILE DEFENSE 

Supersonic 
aerodynamics 

Aerothermodynamics 

High temperature heat 
transfer 

Space physics 

Re-entry programs 


ASTRO SCIENCES 

Space electronics 

Guidance & control 

Communications 

Instrumentation 

Experimental physics 

Plasma and magneto- 
hydrodynamics 
studies 


Visit Aeronutronic’s ex- 
hibit booth 3822-24 at 
the WESCON show. 


Qualified applicants are invited to 
send resumes and inquiries to 
Mr. R. W. Speich, Aeronutronic 
Systems, Inc. 


AERONUTRONIC 


a subsidiary of 

FORD MOTOR COMPANY 
Bidg. 19 Box 451 

Newport Beach, Calif. ORiole 3-2520 


Newport Beach e Santa Ana 
Maywood, California 
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pulsion System,’’ by H. Preston-Thomas, 
J. Brit. Interplanet. Soc., vol. 16, Sept.- 
Oct. 1958, pp. 508-516. 

Multidirectional G-protection in Space 
Vehicles, by Harald J. vonBeckh, J. 
Brit. Interplanet. Soc., vol. 16, Sept.-Oct. 
1958, pp. 525-533. 

Observations at the Mullard Radio 
Astronomy Observatory (on Satellites), 
Cambridge, Proc. Roy. Soc., vol. 248A, Oct. 
28, 1958, pp. 3-9. 

Interferometer Measurements on the 
First Satellites, by A. N. Beresford, Proc. 
Roy. Soc., vol. 248A, Oct. 28, 1958, pp. 
10-15. 

Observations on the USSR Earth 
Satellite Radio Signals, by H. V. Grif- 
fiths, Proc. Roy. Soc., vol. 248A, Oct. 28, 
1958, pp. 16-23. 

Radar Observations of the Russian 
Earth Satellites and Carrier Rocket, by 
The Staff of the Jodrell Bank Exp. Sta., 
Proc. Roy. Soc., vol. 248A, Oct. 28, 1958, 
pp. 23-33. 


Flight Vehicle Design and 
Testing 


Mathematical Theory of the Optimum 
Trajectories of a Rocket, by Angelo Miele, 
Purdue Univ., School of Aeron. Engnq., 
Rep. A-58-10 (AFOSR-T R-58-154; ASTIA 
AD 206361), Nov. 1958, 35 pp. 

Dynamic Stability of Re-entry Vehicles, 
by Norris F. Dow, Astronavtics, vol. 4, 
March 1959, pp. 32-33, 66. 


Aerophysics, Astrophysics 


New Research on Old Gravitation, by 
R. H. Dicke, AEC, NYO-8144, Oct. 1958, 
11 pp. 

Dirac’s Cosmology and the Dating of 
Meteorites, by R. H. Dicke, AEC, NYO- 
8145, Nov. 1958, 4 pp. 

Bakerian Lecture: The Nature of the 
Cosmic Radio Sources, by M. Ryle, 
Proc. Roy. Soc., London. vol. A248, Nov. 
25, 1958, pp. 289-301. 

Discontinuities in Spherically Sym- 
metric Gravitation Fields and Shells of 
Radiation, by W. Israel, Proc. Roy. Soc., 
London, vol. A248, Nov. 25, 1958, pp. 404- 
416. 

On the Interplanetary Gas and Its Mag- 
netic Field, by L. Block, Arkiv f. Fysik, 
vol, 14, no. 2, 1958, pp. 179-193. 

Cosmic Ray Orbits in Interplanetary 
Magnetic Fields, by L. Block, Arkiv f 
Fysik, vol. 14, no. 2, 1958, pp. 161-178 

Interplanetary Magnetic Field and Its 
Control of Cosmic-ray Variations, by 
J. H. Piddington, Phys. Rev., vol. 112, 
Oct. 15, 1958, pp. 589-596. 

Structure of Particles of Linearized 
Gravitational Theory, by R. Sachs and 
P. G. Bergmann, Phys. Rev., vol. 112, 
Oct. 15, 1958, pp. 674-677. 

Project Prairie Grass, a Field Program 
in Diffusion, by Air Force Cambridge Re- 
search Center, Geophysics Res. Director- 
ate, Geophysical Res. Paper 59 (Air 
Force Cambridge Res. Center, Tech. Rep. 
58-235), vol. I, 280 pp. (ASTIA AD 
152272); vol. II, 209 pp. (ASTIA AD 
152573), July 1958. 

Contributions to Stratospheric Meteor- 
ology, by Air Force Cambridge Res. Cen- 
ter, Geophysics Res. Directorate, GRD 
Res. Note 1 (Air Force Cambridge Res. 
Center, TN 58-450; ASTIA AD 152626), 
Aug. 1958, 134 pp. 

Thomas-Fermi Equation of State for 
Dilute Gases, by Richard Latter, Rand 


Corp., Res. Mem. RM 1847-AEC (AST/A 
AD 123511), Jan. 1957, 20 pp. 

Spectroscopic Observations of Venus 
for Rotation Made at Mount Wilson in 
1956, by Robert S. Richardson, J. Brit. 
Interplanet. Soc., vol. 16, Sept.-Oct. 1958, 
pp. 517-524. 

Ionizing Radiation Associated with 
Solar Radio Noise Storm, by K. A. 
Anderson, Phys. Rev., vol. 112, no. 1, 
Oct. 1, 1958, pp. 333-334. 


Atomic Energy 


Cohesion in Plasma, by Melvin A. Cook 
and William S. McEwan, J. Appl. Phy:., 
vol. 29, Nov. 1958, pp. 1612-1613. 

Observations on the Plasma Produced 
in a Magnetic Mirror Geometry, by W. 33. 
Johnson, Phys. Rev. Letters, vol. 1, no. )), 
Nov. 1, 1958, pp. 333-334. 

The Diffusion of Magnetic Fields in a 
Cylindrical Conductor, by J. H. Adlan 
and R. J. Taylor, Gt. Brit., Atom. Ener: y 
Res. Estab., AERE T/M 160, 1958, '2 
pp. 
Tests of a Direct Cycle Nuclear Turbo- 
jet System, by D. R. Shoults, Cincinnati, 
Interavia, vol. 13, Dec. 1958, pp. 1316- 
1317. 

Magnetically Confined Plasmas, by 
A. C. Kolb, Phys. Rev., vol. 112, Oct. 15, 
1958, pp. 291-295. 

Vlasov Instability in Longitudinal 
Plasma Oscillations, by Peter L. Auer, 
Phys. Rev. Letters, vol. 1, no. 11, Dee. |, 
1958, pp. 411-412. 

Our Stake in Thermonuclear Power, by 
John M. Carroll, Electronics, vol. 31, 
Dec. 19, 1958, pp. 75-77. 

USSR May Have Tested Atom Rocket, 
by Alfred J. Zachringer, Missiles and 
Rockets, vol. 4, Dec. 22, 1958, pp. 13-14. 

Studies of Neutron Production in 
Linear Deuterium Pinches, by P. Ohlin, 
K. Siegbahn, T. Sundrstrom and _ §&. 
Svennerstedt, Nuclear Insir., vol. 3, no. 
4, Oct. 1958, pp. 237-244. 

The Mechanism of Neutron Emission in 
High Energy Fission, by D. M. Skyrme 
and G. N. Harding, Nuovo Cimento, vol. 
9, Sept. 16, 1958, pp. 1082-1092. (In 
English.) 

The Runaway Effect in a Fully Ionized 
Plasma, by FE. R. Harrison, Phil. Maqg., 
Ser. 8, vol. 3, Nov. 1958, pp. 1318-1327. 

Safety Hazards of Nuclear Propulsion, 
II, by Leo Seren, Space/Aeron., vol. 30, 
Dec. 1958, pp. 48-52. 

Experiments on the Ohmic Heating and 
Confinement of Plasma in a Stellarator, 
hy T. Coor, S. P. Cunningham, R. A. 
Ellis, M. A. Heald and A. Z. Kranz, Phys. 
Fluids, vol. 1, no. 5, Sept.-Oct. 1958, pp. 
411-420. 

‘“‘Runaway”’ Electrons and Cooperative 
Phenomena in B-1 Stellarator Discharges, 
by W. Bernstein, F. F. Chen, M. A. Heald 
and A. Z. Kranz, Phys. Fluids, vol. 1, no. 
5, Sept.-Oct. 1958, pp. 430-437. 

The Divertor, a Device for Reducing 
the Impurity Level in a Stellarator, by 
C. R. Burnett, D. J. Grove, R. W. Palla- 
dino, T. H. Stix and K. E. Wakefield, 
Phys. Fluids, vol. 1, no. 5, Sept.-Oct. 1958, 
pp. 488-445. 

Interferometric Measurement of Elec- 
tron Concentrations in Plasmas, by Ralph 
A. Alpher and Donald R. White, Phys. 
Fluids, vol. 1, no. 5, Sept.-Oct. 1958, pp. 
452-453. 

What Are the Safety Hazards of Nuclear 
Propulsion? by Leo Seren, Space/Aero- 
nautics (formerly Aviation Age), vol. 30, 
Nov 1958, pp. 18-22. 
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Russian Technical Articles* 


Interaction of a Shock Wave with a 
Boundary Layer by the Leading Edge of a 
Flat Plate at “y-% Supersonic Speeds with 
Radiation, by G. A. Kulonen, Leningrad 
Univer., Vestnik, vol. 13, no. 7, 1958, pp. 
172-188. (In Russian.) 


The Behavior of a Completely Ionized 
Plasma in a Strong Magnetic Field, by 
S. I. Braginskii, Soviet Phys.-J ETP, vol. 6, 
March 1958, pp. 494-498. 

Stationary Convective Flow of an Elas- 
tically Conducting Liquid between Parallel 
Plates in a Magnetic Field, by G. E. Ger- 
shuni and E. M. Zhukhovitskii, Soviet 
Phys.-JETP, vol. 34 (7), no. 3, Sept. 
1°58, pp. 461-464. 

Stability of the Stationary Convective 
Flow of an Electrically Conducting Liquid 
between Parallel Vertical Plates in a Mag- 
netic Field, by G. E. Gershuni and E. M. 
Ziukhovitskii, Soviet Phys.-JETP, vol. 
34 (7), no. 3, Sept. 1958, pp. 465-470. 

Radiation of Plasma in a Magnetic 
Field, by B. A. Trubnikov, Akad. Nauk 
SSSR, Doklady, vol. 118, Feb. 11, 1958, 
pp. 913-916. (In Russian. ) 

Unsteady Flow with Heat Transfer in a 
Viscous Incompressible Fluid between 
Two Rotating Disks when There is Injec- 
tion of the Fluid, by G. A. Tirskii, Akad. 
Nauk SSSR, Doklady, vol. 119, March 11, 
1958, pp. 226-228. (In Russian.) 

Stationary Convection in a Plane Liquid 
Layer Near the Critical Heat Transfer 
Point, by L. P. Gor’kov, Soviet Phys.- 
JETP, vol. 6, Feb. 1957, pp. 311-315. 

Study of Heat Transfer in a Spiral 
Canal, by G. Fastovsky and A. FE. Rovin- 
sky, Teploenergetica, no. 1, Jan. 1957, pp. 
39-42. (In Russian.) 

Some Problems of Magnetogasdynam- 
ics with Account of Finite Conductivity, 
by G. 8. Golitsyn and D. P. Staniukovich, 
Soviet Phys.-JETP, vol. 6, June 1958, pp. 
1090-1093. 

Investigations of the Burning of Homo- 
geneous Mixtures in Turbulent Streams 
by Recording Fluctuating Temperatures, 
by I. V. Kokushkin, Akad. Nauk SSSR, 
Izvestiia, Otdelenie Tekh. Nauk, no. 8, Aug. 
1958, pp. 1-11. (In Russian. ) 

On Drag Coefficients in the Motion of 
Burning Particles, by V. N. Basov and V. 
A. Popov, Akad. Nauk SSSR, [zvestiia, 
Oldelente Tekh. Nauk, no. 8, Aug. 1958, 
pp. 12-14. (In Russian.) 

The Phase Diagram of the Hydrogen- 
Deuterium System, by V. 8. Kogan, B. G. 
Lazarev and R. F. Bulatova, Soviet Phys.- 
JETP, vol. 34 (7), no. 1, July 1958, pp. 
165-166. 

The Radiation from an Electron Moving 
in a Magnetoactive Plasma, by V. Ia. 
Kidman, Soviet Phys.-JETP, vol. 34 (7), 
no. 1, July 1958, pp. 91-95. 

The Effect of Coulomb Correlations on 
the Spectrum of Electron-Plasma Oscil- 
lations, by P. 8. Zyrianov, Soviet Phys.- 
JETP, vol. 34 (7), July 1958, pp. 160-161. 

On the Theory of the Stability of Liquid 
Jets in an Electric Field, by G. A. Glonti, 
Zhurnal Eksper. i Teoret. Fiziki, vol. 34, 


*The Editors of Technical Literature 
Digest are making a systematic search for 
pertinent Russian articles, both in the orig- 
inal Russian and in translated form. For 
the balance of 1959, these will appear in a 
separate section with this heading, mainly 
‘o draw the readers’ attention to the new 

sting. After that, they will appear item 
‘y item under the proper subject headings. 
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no. 5, 1958, pp. 1329-1330. (In Russian.) 

Distribution of Momentum in a Con- 
tinuous Fluid Jet at a Supersonic Velocity, 
by A. A. Semerchan, L. F. Vereshchagin, 
F. M. Filler and N. N. Kuzik, Zhurnal 
Tekhn. Fiziki, vol. 28, Sept. 1958, pp. 
2062-2071. (In Russian.) 


Acceleration of Plasma in a Magnetic 
Field, by G. V. Gordeyev and A. L. 
Gubanov, Zhurnal Tekhn. Fiziki, vol. 28, 
Sept. 1958, pp. 2046-2054. (In Russian.) 

Weak Waves in a Compressible Fluid 
With Account of Emission, by V. A. Pro- 
kof’yev, Prikladnaya Matematika i Mekh- 
anika, vol. 21, Nov.-Dec. 1957, pp. 775- 
782. (In Russian.) 

Ion Gun Operating at 20 to 30 Kilo- 
volts, by L. N. Bukhovskaya, Pribory 7 
Tekhnika Eksperimenta, no. 2, March- 
April 1958, pp. 49-51. (In Russian.) 

Heat Transfer to Liquid Metals, by S. S. 
Kutateladze, EK. N. Borishanskiy and 
I. Novikov, Atomnaya Energiya, vol. 4, 
May 1958, pp. 422-436. (In Russian.) 

On Riemann Waves in Magnetohydro- 
dynamics, by A. G. Kulikovskiy, Doklady 
Akadamit Nauk SSSR, vol. 121, Aug. 21, 
1958, pp. 987-990. (In Russian.) 

Investigation of the Kinetics of the 
Combustion of Hydrogen With Oxygen 
Above the Lower Limit of Self-ignition, 
by L. V. Karmilova, A. B. Nalbandyan 
and N. N. Semenov, Zhurnal Fiz. Khimii, 
vol. 32, June 1958, pp. 1193-1204. (In 
Russian. ) 

On the Stability of Flame Propagation, 
by Ya. B. Zel’dovich and G. I. Baren- 
blatt, Prikladnaya Matematika i Mekh- 
anika, vol. 21, no. 6, 1957, pp. 856-859. 
(In Russian.) 

The Interaction Between ee and 
Methyl Hydroperoxide, by N. A. Kley- 
menov and A. B. Nalbandyan, ‘poblady 
Akademii Nauk SSSR, vol. 118, Jan. 
1958, pp. 125-127. (In Russian.) 

On the Combustion of and Heat Trans- 
fer by Dispersed Liquid Fuel in a Turbu- 
lent Air Flow, by G. N. Delyagin and 
B. V. Kantorovich, [nzhenerno-F izicheskiy 
Zhurnal, no. 3, March 1958, pp. 24-39. 
(In Russian.) 

The Reactions of Nitrogen in Explo- 
sions, by A. Ya. Apin, Yu. A. Lebedova 
and O. I. Nefedova. Zhurnal Fiz. Khimii, 
vol. 32, April 1958, pp. 819-823. (In 
Russian.) 

The Glowing of Air During a Strong 
Explosion and the Minimum Brightness 
Effect of a Fireball, by Yu. P. Rayzer, 
Zhurnal Eksper. i Teoret. Fiziki, vol. 34, 
no. 2, 1958, pp. 483-493. (In Russian.) 


Measurement of Ionization Intensity in 
(Acetylene) Flames, by V. S. Ressikhin 
and N. A. Nesterko, Zhurnal Fiz. Khimii, 
vol. 31, no. 12, 1957, pp. 2663-2667. (In 
Russian.) 

Pressure Impact on the Lower Concen- 
tration Limit of the Flame Disruption in a 
Flow, by M. A. Peshkin, Zhurnal Fiz. 
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